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PREFACE 


This publication is a collection of the presentations given at the NASA Computational Fluid 
Dynamics (CFD) Conference held at NASA Ames Research Center, Moffett Field, California, 
March 7-9, 1989. The objectives of the conference were to disseminate CFD research results to 
industry and university CFD researchers, to promote synergy among NASA CFD researchers, and to 
permit feedback from researchers outside NASA on issues pacing the discipline of CFD. The focus 
of the conference was on the application of CFD technology but also included fundamental 
activities. The conference was sponsored by the Aerodynamics Division, Office of Aeronautics and 
Space Technology (OAST), NASA Headquarters, Washington, DC 20546. 

The conference consisted of twelve sessions of papers representative of CFD research conducted 
within NASA and three non-NASA panel sessions. For each panel session, the panel membership 
consisted of industry and university CFD researchers. A summary of the comments made during the 
panel sessions have been included in this publication. 

The conference proceedings are published in two volumes. Volume 1 contains the papers presented 
in Sessions I- VI; Volume 2 contains those given in Sessions VII-XII. Each volume contains the 
same front matter, and each contains a list of attendees as an appendix. 
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PANEL SESSION SUMMARY 



NASA CFD Conference 
NASA Ames Research Center 
March 7-9, 1989 

Panel Sessions Summary 

The NASA CFD Conference was held at Ames Research Center on March 7-9, 
1989. To conclude each day’s presentations, a panel session with participation 
from the audience furnished a great deal of excellent feedback from the industry 
and academic communities. During the conference it was evident that the panel 
members proffered comments only after having spent considerable time in 
preparing them. 

The members of the panel sessions are listed below: 

March 7 P. Rubbert - Boeing Commercial Airplanes 

R. Melnik - Grumman Aerospace Corporation 
D. Whitfield - Mississippi State University 

March 8 I. Bhateley - General Dynamics - Fort Worth Division 

R. Agarwal - McDonnell Douglas Research Laboratories 
R. MacCormack - Stanford University 

March 9 V. Shankar - Rockwell International Science Center 
J. Carter - United Technologies Research Center 
A. Jameson - Princeton University 

The crucial comments from the three panel sessions have been combined and are 
summarized as follows: 

- NASA's CFD program is now too heavily focused on applications: program 
balance has swung from fundamentals (1970's) to applications (1980's) 

- Three critical "needs" emerged: 

(1) More algorithm research is needed; especially for Navier-Stokes 
solvers with unstructured grids 

(2) More research is required on geometric modelling; need rapid, 
accurate, and effective surface definition techniques 

(3) More research is needed on grid generation methods with the focus 
on speed, efficiency, and grid quality to reduce set up time and 
complexity 

- Developers of CFD need to understand the needs of the users; designers of 
aerospace vehicles have requirements that are different than the CFD 
researchers perceptions 

- Industry needs more reliable and cost effective CFD tools 
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Additional detail comments from the three panel sessions are listed below: 

- CFD has matured during the last decade and is being used to solve real 
problems; however, industry lacks confidence in Navier-Stokes solutions 

- Industry needs codes that have been validated to increase confidence in 
CFD technology 

- Improved communality between codes would increase usability; standards 
are needed 

- Improved data storage, networking, data transfer, and graphics required to 
assimilate information provided by CFD 

- Improved turbulence modeling for separated flows 

- Accurate prediction of drag for complete powered aerospace vehicles 

- Develop multidisciplinary CFD technology with optimization capability 

- NASA must maintain focus on technology development and high risk 
research 

- Technology transfer is not complete until design engineers are using CFD 
codes successfully 

- Industry needs NASA to improve CFD technology for codes simpler than 
Navier-Stokes solvers 

- NAS program has been extremely helpful to industry in transferring CFD 
technology 

- Industry needs to be more aggressive in their use of CFD 

- Improved understanding of CFD by design engineers required; cooperative 
programs or workshops were suggested to bring CFD researchers and 
designers together 

- Design cycle time needs to be reduced with CFD; codes must be cost 
effective, reliable, and useable, and robust to work at flight Reynolds 
Numbers 

- Improved coordination/reduced overlap of CFD applications between NASA 
centers 
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SESSION I 


CENTER OVERVIEWS 


Chairman: 

Randolph A. Graves, Jr. 

Director, Aerodynamics Division 
NASA Headquarters 



COMPUTATIONAL FLUID DYNAMICS PROGRAM 
AT NASA AMES RESEARCH CENTER 


Terry L. Holst 

Chief, Applied Computational Fluids Branch 
NASA Ames Research Center, Moffett Field, California 


ABSTRACT 

The Computational Fluid Dynamics (CFD) Program at NASA Ames Research Center is 
reviewed and discussed. The presentation is broken into several sections as follows: First, the 
technical emements of the CFD Program are generally listed and briefly discussed. These 
elements include algorithm research, research and pilot code development, scientific 
visualization, advanced surface representation, volume grid generation, and numerical 
optimization. Next, the discipline of CFD is briefly discussed and related to other areas of 
research at NASA Ames including Experimental Fluid Dynamics, Computer Science Research, 
Computational Chemistry, and Numerical Aerodynamic Simulation. These areas combine with 
CFD to form a larger area of research, which might collectively be called computational 
technology. The ultimate goal of computational technology research at NASA Ames is to 
increase the physical understanding of the world in which we live, solve problems of national 
importance, and increase the technical capabilities of the aerospace community. 

Next, the major programs at NASA Ames that either use CFD technology or perform research 
in CFD are listed and discussed. Briefly, this list includes turbulent/transition physics and 
modeling, high-speed real gas flows, interdisciplinary research, turbomachinery demonstration 
computations, complete aircraft aerodynamics, rotorcraft applications, powered lift flows, high 
alpha flows, multiple body aerodynamics, and incompressible flow applications. Some of the 
individual problems actively being worked in each of these areas is listed to help define the 
breadth or extent of CFD involvrnent in each of these major programs. 


State-of-the-art examples of various CFD applications are presented to highlight most of these 
areas. The main emphasis of this portion of the presentation is on examples which will not 
otherwise bd treated at this conference by the individual presentations. Thus, a good survey of 
CFD applications research at NASA Ames can be obtained by looking at this presentation in 
conjunction with the individual NASA Ames presentations made at this conference. 

Finally, this overview is concluded with a list of principal current limitations and expected 
future directions. Some of the future directions include algorithm research, turbulence/transition 
research, multidisciplinary research, graphics and workstation research and applications which 
will address more realistic simulations in the engineering world. 
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MAJOR PROGRAMS USING CFD 

NASA Ames Research Center 


TURBULENT/TRANSITION PHYSICS AND PHYSICAL MODELING 
HIGH-SPEED REAL GAS FLOWS 

• THERMO- AND CHEMICAL-NONEQUILIBRIUM 

• RADIATION 

• COMBUSTION 

• RAREFIED FLOW EFFECTS 

INTERDISCIPLINARY RESEARCH 
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• CFD + COMPUTATIONAL STRUCTURAL MECHANICS 

• CFD + ACTIVE CONTROLS 
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TURBOMACHINERY DEMONSTRATION COMPUTATIONS 
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• JET FREESTREAM MIXING 

• THRUST AUGMENTOR EJECTORS 

• STOVL DELTA WING IN GROUND EFFECT 


HIGH ALPHA 

• HARV APPLICATIONS (F-1 8) 

• OGIVE CYLINDER COMPUTATIONS 

• UNSTEADY FLOWS 



MAJOR PROGRAMS USING CFD (CONCLUDED) 

NASA Ames Research Center 


MULTIPLE BODY AERODYNAMICS 

• SPACE SHUTTLE (LAUNCH CONFIGURATION) 

• SRB/ET-ORBITER SEPARATION 

• AIRCRAFT STORE SEPARATION 

• SPACE SHUTTLE C/ SPACE SHUTTLE II 


INCOMPRESSIBLE NAVIER-STOKES 

• SSME APPLICATIONS 

• HYDRODYNAMICS 

• HIGH LIFT CONFIGURATIONS 

• ARTIFICIAL HEART BLOOP FLOV^ SIMULATION 



ADVANCED SIMULATION AND ANALYSIS 

PROJECT (ASAP) 

VAN DALSEM, VOGEL, LUH, SORENSON, ATWOOD 

OBJECTIVE 

• REDUCE THE “CLOCK TIME” 

REQUIRED TO OBTAIN THE SURFACE 
DEFINITION AND GRID ABOUT A 
COMPLEX CONFIGURATION BY AT 
LEAST AN ORDER OF MAGNITUDE 

APPROACH 

* • DEVELOP AN INTEGRATED, 

INTERACTIVE SURFACE DEFINITION 
AND GRID GENERATION CAPABILITY 
TAILORED TO THE CFD ENVIRONMENT 

FUTURE DIRECTIONS 

• EXPLORE APPLICATION OF Al TO ENHANCE NONEXPERT USER PERFORMANCE 

• INVESTIGATE: 

- GRID QUALITY MEASURES 

- SOLUTION-ADAPTIVE TECHNIQUES 

- NEW GRID GENERATION APPROACHES (STRUCTURED AND UNSTRUCTURED) 

PAYOFF 

• A POWERFUL, EASY-TO-USE TOOL THAT SIGNIFICANTLY REDUCES THE 
TURNAROUND TIME FOR CURRENT AND FUTURE CFD ANALYSES 



ADVANCED SIMULATION AND ANALYSIS 

PROJECT (ASAP) 

VOGEL, LUH, SORENSON, ATWOOD 



F-18 FOREBODY GRID BY 3DGRAPE: SELECTED AXIS-NORMAL SURFACES 



SURFACE GRID FOR GENERIC 2D ZONING COMPARISON: HUMAN 

HYPERSONIC AIRPLANE EXPERT vs EXPERT SYSTEM (EZGRID) 




PANEL METHOD APPLICATIONS 

ASHBY, IGUCHI, BROWN 


OBJECTIVES 


• DEVELOP CAPABILITY TO ANALYZE 
COMPLEX GEOMETRIES VERY QUICKLY 

-INCLUDES LEADING-EDGE SEPARATION, 
JET PLUMES, AND UNSTEADY EFFECTS 
(TIME STEPPING) 

= APPROACH 

• LOW-ORDER PANEL METHOD WITH TIME- 
STEPPED WAKES 

FUTURE DIRECTIONS 

• COUPLE WITH BOUNDARY LAYER CODE 
TO INCLUDE VISCOUS EFFECTS 

• USE FOR DESIGN AND ANALYSIS OF WIND 
TUNNEL MODELS AND TO DETERMINE 
WIND TUNNEL WALL INTERFERENCE 

PAYOFF 

• EFFICIENT, RELIABLE TOOL FOR USE IN 
LOW SPEED APPLICATIONS 





PANEL METHOD APPLICATIONS 

ASHBY, IGUCHI, BROWN 



80- BY 1 20-FOOT WIND TUNNEL INLET 



VELOCITY VARIATION ACROSS INLET 



E-7 STOVL MODEL 



CHORDWISE PRESSURE DISTRIBUTION 
ON E-7 WING PMARC COMPUTATION; 
2y/b = 0.6, a = 8° 


u> 


TWO-DIMENSIONAL COMPUTATIONS OF 
MULTI-STAGE COMPRESSOR FLOWS 


GUNDY-BURLET, RAI 


OBJECTIVE 

• DEVELOP CAPABILITY TO CALCULATE 
UNSTEADY VISCOUS FLOWS WITHIN 
MULTI-STAGE TURBOMACHINES 

CURRENT APPROACH 

• SOLVE THE TWO-DIMENSIONAL 
NAVIER-STOKES EQUATIONS USING A 
ZONAL METHOD TO SIMULATE ROTOR- 
STATOR INTERACTION 


FUTURE DIRECTIONS 

• EXTEND TO THREE-DIMENSIONS 

PAYOFF 

• BETTER UNDERSTANDING OF 
UNSTEADY FLUID DYNAMICS IN 
TURBOMACHINES 

• INCREASED RELIABILITY AND 
EFFICIENCY OF TURBOMACHINES 


INSTANTANEOUS TEMPERATURE CONTOURS 
Mqo = 0.07, Re = 100,000/in. 






EFFECT OF TANGENTIAL LEADING EDGE 
BLOWING ON VORTICAL FLOW 

YEH, TAVELLA, ROBERTS 


OBJECTIVE 

• TO INVESTIGATE THE ABILITY OF 
TANGENTIAL LEADING EDGE 
BLOWING TO CONTROL VORTICAL 
FLOW AT HIGH ALPHA 

APPROACH 

• UTILIZE DELTA WING GEOMETRY 

• SOLVE THIN-LAYER NAVIER-STOKES 
EQUATIONS USING MULTIPLE-ZONE GRID 
APPROACH TO ACCOMODATE JET-SLOT 
GEOMETRY 

« UTILIZE ALGEBRAIC TURBULENCE 
MODEL FOR SURFACE BL AND WALL JET 

FUTURE DIRECTIONS 

• EXTEND TO FULL AIRCRAFT CONFIGURATIONS 

• INVESTIGATE BLOWING CONTROL CONCEPTS 

PAYOFF 

• INCREASED UNDERSTANDING OF VORTICAL 
FLOW PHYSICS 

• NEW TOOL FOR STUDYING BLOWING 
CONTROL CONCEPTS 


COMPUTED “OIL FLOW’ 

ON DELTA WING SURFACE 
Moo = 0.3, Re = 1.3 x 106, a 

BLOWING 




EFFECT OF TANGENTIAL LEADING EDGE 
BLOWING ON VORTICAL FLOW 

YEH, TAVELLA, ROBERTS 
(Moo = 0.3, a = 40®, Re = 1.3 x 106, X/C = 0.36) 


NORMALIZED 

TOTAL 

PRESSURE 

CONTOURS 


SURFACE 

PRESSURE 

DISTRIBUTION 


o EXPERIMENT 
CALCULATED 


NO BLOWING 




BLOWING 




0 .2 .4 .6 .8 1 

Y/S 






EULER VALIDATION/PRESSURE INTEGRATION 

MELTON, ROBERTSON, MOYER 


OBJECTIVES 


M =08° EULER FL057 

6 ‘ A EXPERIMENTAL (INTEGRATED Cp’s + SKIN FRICTION) 
\~ □ EXPERIMENTAL (BALANCE) 


• CFD VALIDATION FOR FL057 
EULER CODE 

• ENHANCE WIND TUNNEL 
PRESSURE INTEGRATION FOR 
PREDICTING FORCES AND 
MOMENTS 

APPROACH 

• FL057 FINITE VOLUME 3D EULER 
CODE 

• COMPUTE DISCRETIZATION 
ERROR BY COMPARING CFD 
FORCE AND MOMENT 
INTEGRATION WITH CFD 
PRESSURES INTERPOLATED AND 




0 .050 .100 .02 

C D 


0 

C M 


INTEGRATED AT MODEL TAP LOCATIONS 



FUTURE DIRECTIONS 

• INVESTIGATE NEW METHODS FOR INTEGRATING CFD AND EXPERIMENTAL RESULTS 

PAYOFF 

• VALIDATION OF FLO 57 FOR DELTA CONFIGURATIONS 

• REDUCE INSTRUMENTATION CONSTRAINTS ON COMPLEX WIND TUNNEL MODELS 

• INCREASE ACCURACY OF FORCE AND MOMENT PREDICTIONS FROM WIND 
TUNNEL PRESSURE INTEGRATIONS 



EULER VALIDATION/PRESSURE INTEGRATION 

MELTON, ROBERTSON, MOYER 

M = 0.8 DRAG POLAR 


00 


UPPER SURFACE Cp’s 



.5 


.3 


o 


.2 


.1 


-.1 


EXPERIMENTAL BALANCE DATA 

□ EXPERIMENTAL PRESSURE 
INTEGRATION 

o EXP PRESSURE INTEGRATION 
CORRECTED FOR 
DISCRETIZATION ERROR 

a FINAL EXP PRESSURE 
INTEGRATION CORRECTED 
FOR DISCRET ERROR AND 
SKIN FRICTION 


'SKIN FRICTION CORRECTION 
DISCRETIZATION CORRECTION 


.01 


.02 .03 


.04 


.05 


.06 


.07 


-.2 


SPACE SHUTTLE LAUNCH CONFIGURATION 

STEGER, RIZK, OBAYASHI, MARTIN, CHIU, BUNING 


OBJECTIVE 

• DEVELOP CAPABILITY TO 
COMPUTE FLOW OVER 
INTEGRATED SPACE SHUTTLE 
IN ASCENT 

APPROACH 

• SOLVE 3D REYNOLDS- 
AVERAGED NAVIER-STOKES 
EQUATIONS 

• USE CHIMERA GRID APPROACH 



FUTURE DIRECTIONS 

• IMPROVE PLUME SIMULATION CAPABILITY AND CODE EFFICIENCY 

• VALIDATE UNSTEADY MODE AND STUDY FAST SEPARATION SIMULATIONS 

PAYOFF 

• PREDICTIVE TOOL FOR UNDERSTANDING AND REFINING AERODYNAMIC 
PERFORMANCE OF MULTIPLE BODY VEHICLES 


SPACE SHUTTLE ASCENT-MODE RESULTS 


to 

o 



SURFACE PRESSURES 
ON SHUTTLE C 



g 1 i i i i i i i I i i- 1 1 1 1 — 

'900 1300 1700 2100 900 1300 1700 2100 


X X 

Cp COMPARISONS ALONG 
ORBITER FUSELAGE DURING ASCENT 



a = -3 WIND TUNNEL 

COMPARISON OF SURFACE PRESSURES 


UNSTEADY MULTIPLE BODY AERODYNAMICS 




OBJECTIVE 

• TO DEVELOP A GENERAL CAPABILITY FOR TIME- 
ACCURATE SIMULATION OF 3-D MULTIPLE BODY 
VISCOUS FLOWS GIVEN ARBITRARY GRID 
COMBINATIONS, BODY SHAPES, AND RELATIVE 
MOTION BETWEEN GRID SYSTEMS 


CURRENT APPROACH 

• UNSTEADY CHIMERA COMPOSITE GRID 
TECHNIQUES 

• IMPLICIT TIME-ACCURATE SOLVER FOR THE THIN- 
LAYER NAVIER-STOKES EQUATIONS 

FUTURE DIRECTIONS 

• DEVELOP TRAJECTORY PREDICTION ROUTINES 

• IMPROVE EFFICIENCY AND ACCURACY OF BASIC 
ALGORITHMS 

• CODE VALIDATION STUDIES 


PAYOFF 

• A VALIDATED COMPUTATIONAL TOOL FOR 
ANALYZING COMPLEX AERODYNAMIC PROBLEMS 
INVOLVING MULTIPLE BODIES IN RELATIVE MOTION 




TIME-ACCURATE SIMULATION OF THE SPACE SHUTTLE SRB 

SEPARATION SEQUENCE 


PRESSURE CONTOURS 

4.5 

a = + 2 ° 

R = 6.95 x10 s 

e 


IO 

to 


ASSUMPTIONS: 

• SIMPLIFIED 
GEOMETRY 

• NO PLUMES 

• PRESCRIBED 
SRB 

TRAJECTORY 


(MEAKIN, SUHS) 



DISCRETIZATION 

GRIDS: 

ET 73X39X45 
SRB 53X37X21 
ORB 74X77X33 


TIME-STEP = 

1 .36 X 1 0 3 sec 

500 STEPS THROUGH BSM 

BURN 

BOOSTER 

SEPARATION 

MOTOR 

(BSM) 

BURN-TIME = 0.68 sec 








AIRCRAFT STORE SEPARATION 

MEAKIN, SUHS 

FREESTREAM CONDITIONS: M =1.05, a = +2°, Re = 2.4x10 6 

oo 7 7 

WING LOWER SURFACE C p DISTRIBUTION MACH CONTOURS ABOUT STORE 

TIME-ACCURATE 

STEADY-STATE STORE-SEPARATION SIMULATION 



TURBULENCE MODELING FOR HYPERSONIC FLOWS 

COAKLEY, HORSTMAN, KUSSOY, MARVIN 


OBJECTIVE 

• IMPROVE AND DEVELOP MODELS 
FOR HYPERSONIC FLOWS 


CURRENT APPROACH 

• PERFORM COMBINED 
COMPUTATIONAL AND 
EXPERIMENTAL STUDIES 


TEST MODEL GEOMETRIES USED FOR 
COMPARISONS OF COMPUTATION AND 
EXPERIMENT 



FUTURE DIRECTIONS 

• IMPROVE COMPUTATIONAL EFFICIENCY 

• DEVELOP SECOND ORDER CLOSURE MODELS 

• PERFORM EXPERIMENTS USING NEWLY DEVELOPED NON INTRUSIVE 
INSTRUMENTATION 


PAYOFF 

• ACCURATE COMPUTATIONS OF HEAT TRANSFER, SKIN FRICTION, AND 
COMPLEX FLOW STRUCTURES 



ORIGINAL PAGE IS 

OF POOR QUALITY 


K) 


PRESSURE CONTOURS FOR AN IMPINGING 

SHOCK WAVE FLOW 


Moo = 7.2 

BALDWIN-LOMAX MODEL 


e = 15° 



q-co MODEL c 



1.0 1.2 1.4 


PRESSURE CONTOURS FOR AN IMPINGING 

SHOCK WAVE FLOW 

Moo =7.2 e = 15° 


COMPUTATION q-co MODEL c 


EXPERIMENT 



HYPERSONIC APPLICATIONS 


OBJECTIVE 

• DEVELOP CAPABILITY TO COMPUTE REAL- 
GAS AEROTHERMODYNAMIC 
CHARACTERISTICS OF HYPERSONIC 
VEHICLES 

• USE CAPABILITY TO GUIDE VEHICLE DESIGNS 

APPROACH 

• SOLVE 3D REYNOLDS-AVERAGED NAVIER- 
STOKES AND PARABOLIZED NAVIER-STOKES 
EQUATIONS USING VARIOUS TRANSITION/ 

TURBULENCE MODELS WITH PERFECT GAS, 

EQUILIBRIUM REAL GAS AND 
NONEQUILIBRIUM REAL GAS MODELS 

FUTURE DIRECTIONS 

• IMPROVE TRANSITION/TURBULENCE MODELS, 

REAL GAS MODELS, AND COMPUTATIONAL 
EFFICIENCY 

• EXTEND APPLICATIONS TO MORE COMPLEX GEOMETRIES 

PAYOFF 

• PROVIDE DESIGN INFORMATION NOT POSSIBLE TO MEASURE IN GROUND-BASED 
TEST FACILITIES 

• ENABLE DEVELOPMENT OF AEROASSISTED VEHICLES AND AIRBREATHING 
HYPERSONIC AIRCRAFT 


NORMALIZED PRESSURE 


GENERIC HYPERSONIC 
AEROTHERMODYNAMIC RESULTS 

LAWRENCE 

= 11.4, a = 0°, Re L = 29 X 106 


00 


SYMMETRY PLANE 

PRESSURE CONTOURS 



WINDWARD SYMMETRY PLANE PRESSURES 


WINDWARD SYMMETRY PLANE HEAT TRANSFER 





HYPERSONIC EXHAUST PLUME/AFTERBODY INTERACTION 

EDWARDS 


NOZZLE/AFTERBODY MODEL 



EXHAUST GAS CONCENTRATION 
CONTOURS 

M oo = 6, NPR = 4, a = 0° 


o 



AFTERBODY PRESSURE COMPARISON 
BETWEEN 2ND AND 3RD NOZZLES 



GENERIC HYPERSONIC VEHICLE 


PRESSURE CONTOURS 
6, NPR = 100, a = 0° 



M 00 = 1.46, NPR = 6.27, a = 0° 
□ EXPERIMENT 



x, cm 

PRESSURE COMPARISON 
WINDWARD PLANE OF SYMMETRY 
DOWNSTREAM OF NOZZLE 



BLUNT 5° CONE WITH SHOCK GENERATORS 

Molvik,Strawa 


o 


Conditions: M ^ = 14.4 
Re L = 10 6 
Poo = 0.0932atm 
Poo = 298° K 
a = 6.35° 

Experimental Facility: Ames Ballistic Range 

Flow Solvers: Three-Dimensional Navier- Stokes 
with Finite Rate Chemistry 
(TUFF and STUFF) 

REAL GAS COMPUTATION 




REAL GAS COMPUTATION 



ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 


DIRECT PARTICLE SIMULATION OF 
HYPERSONIC FLOWS 


OBJECTIVE: 

• DEVELOP THE CAPABILITIES OF A NEW DISCRETE 
PARTICLE SIMULATION METHOD FOR RAREFIED 
HYPERSONIC FLOWS IN 3D WITH NON- 
EQUILIBRIUM CHEMISTRY. 

APPROACH: 

• FLUID IS MODELED AS A LARGE COLLECTION OF 
DISCRETE PARTICLES THAT INTERACT WITH EACH OTHER 
THROUGH COLLISIONS. 

• SIMPLIFIED PHYSICAL MODELS ARE USED ALLOWING 
ORDERS OF MAGNITUDE INCREASE IN COMPUTATIONAL 
EFFICIENCY WHILE ENHANCING STATISTICAL ACCURACY. 

FUTURE DIRECTIONS: 

• REALISTIC 3D GEOMETRIES WITH MORE GENERAL BOUNDARY CONDITIONS. 

• EXTENDED MOLECULAR MODELS TO ACCOUNT FOR ADDITIONAL INTERNAL DEGREES OF 
FREEDOM, CHEMISTRY AND WALL-PARTICLE INTERACTIONS. 

PAYOFF: 

• DIRECT PARTICLE SIMULATION IS APPLICABLE AT LOW DENSITIES AND HIGH MACH 
NUMBERS BEYOND THE REACH OF CONTINUUM METHODS. 

• ENABLES PARTICLE SIMULATIONS ON A MUCH LARGER SCALE THAN 
PREVIOUSLY POSSIBLE. 

• PROVIDES NEEDED INSIGHT IN THE DESIGN OF PROPOSED HYPERSONIC VEHICLES. 



PRINCIPAL CURRENT LIMITATIONS 


PHYSICAL MODELS/ALGORITHMS 

• BOUNDARY LAYER TRANSITION MODELS 

• TURBULENCE MODELS FOR SEPARATING AND REATTACHING FLOWS 

• TRANSITION/TURBULENCE MODELS FOR REAL GAS FLOWS AND FLOWS WITH 
COMBUSTION 

• REAL GAS FLOW VALIDATION DATA 

• FAST, USER FRIENDLY GEOMETRY DEFINITION/GRID GENERATION SOFTWARE 

• FAST, ACCURATE ALGORITHMS FOR COMPLETE SIMULATIONS 

• SCIENTIFIC VISUALIZATION SOFTWARE 

COMPUTER SYSTEMS 

• COMPUTATIONAL SPEED 

• NETWORK BANDWIDTHS 

• HIGH-SPEED LARGE-VOLUME MASS STORAGE 

• TOOLS FOR ANALYZING MASSIVE RESULT FILES 



FUTURE DIRECTIONS 


ALGORITHM RESEARCH 

. IMPROVED ALGORITHMS FOR COMPUTING REAL-GAS TURBULENT FLOWS 
. NEW ALGORITHMS TO EXPLOIT ADVANCED MULTIPLE PROCESSOR 

COMPUTER ARCHITECTURES 

• NEW GRID-GENERATION CONCEPTS FOR COMPLEX CONFIGURATIONS, 
MULTIPLE MOVING BODIES AND UNSTEADY FLOWS 

TURBULENCE RESEARCH 

• IMPROVED TURBULENCE MODELS FOR PERFECT-GAS SEPARATING AND 
REATTACHING FLOWS 

• NEW TURBULENCE MODELS FOR REAL-GAS FLOWS 

• METHODS FOR MANAGING TURBULENCE TO REDUCE DRAG, IMPROVE 
COMPONENT PERFORMANCE, MINIMIZE HEAT TRANSFER, AND CONTROL 
COMBUSTION PROCESSES 

MULTIDISCIPLINARY RESEARCH 

. NUMERICAL METHODS FOR SOLVING FULLY COUPLED COMBINATIONS OF 
EQUATIONS FOR AERODYNAMICS, GAS CHEMISTRY, STRUCTURES, 
CONTROLS, PROPULSION AND ELECTROMAGNETICS 



FUTURE DIRECTIONS 

(CONCLUDED) 


GRAPHICS AND WORKSTATION RESEARCH 

• IMPROVED USER EFFICIENCY THROUGH ADVANCES IN GRAPHICS AND 
WORKSTATION TECHNOLOGY 


APPLICATIONS CODES 

• AIRCRAFT MANEUVERING NEAR PERFORMANCE BOUNDARIES 

• POWERED LIFT AIRCRAFT OPERATING IN AND OUT OF GROUND EFFECT 

• HYPERSONIC VEHICLES INCLUDING INLET, ENGINE AND EXHAUST FLOWS 

• ROTORCRAFT IN HOVER, TRANSITION AND FORWARD FLIGHT 

• TURBOMACHINERY INCLUDING PUMPS, COMPRESSORS AND TURBINES 

• METHODS FOR NUMERICALLY OPTIMIZING DESIGNS 

• DESIGNER-FRIENDLY CODES WITH 'EXPERT SYSTEMS' ELEMENTS 



COMPUTATIONAL FLUID DYNAMICS RESEARCH AND APPLICATIONS 
AT NASA LANGLEY RESEARCH CENTER 

Jerry C. South, Jr. 

Head, Analytical Methods Branch 
Fluid Mechanics Division 


Research at Langley includes all aeronautics disciplines and selected space disciplines. Nearly all the 
aeronautics disciplines have a significant component of CFD research which complements wind-tunnel and flight 
experimental research; in space research, aerothermodynamics of planetary entry vehicles is heavily dependent 
on CFD as a research tool. Langley's CFD strategy contains four major thrusts: Focus efforts on critical CFD 
barriers; focus efforts on critical aerodynamics barriers; validate CFD codes; and transfer technology to United 
States users. The Langley presentations in this conference are representative of our strategy. They are grouped 
in six broad areas: 

1. Direct Simulation of Transition and Turbulence 

2. Hypervelocity Aerothermodynamics and Rarified Flows 

3. Hypersonic External and Internal (scramjet) Flows 

4. Unsteady Aerodynamics and Aeroelasticity Applications 

5; Grid Generation and Applications for Complex Configurations 

6. Supersonic and Transonic Wing Design Applications 

Three examples of important work not shown in this conference are described in the conclusion of the 
presentation. 

Or 
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LANGLEY CFD STRATEGY 


FOCUS EFFORT ON CRITICAL CFD BARRIERS 

• GRIDS & ALGORITHMS FOR COMPLEX CONFIGURATIONS 

• TRANSITION & TURBULENCE MODELS FOR RANS CODES 

• NEW ALGORITHMS FOR MASSIVELY-PARALLEL 
PROCESSORS 

• ENABLE ROUTINE APPLICATIONS 

FOCUS EFFORT ON CRITICAL AERODYNAMICS BARRIERS 

• DIRECT SIMULATION OF TRANSITION AND TURBULENCE 

• SIMULATION/PREDICTION OF HIGH-ALPHA FLOWS 

• SIMULATION/PREDICTION OF HYPERSONIC PROPULSION 


LANGLEY CFD STRATEGY 
(Concluded) 


U> 


• VALIDATE CFD CODES 

• CODE-ON-CODE 

• CODE-ON-EXPERIMENT (GROUND & FLIGHT) 

• TRANSFER TECHNOLOGY 

• RESEARCHER EXCHANGES (ARC*-»LaRC-. 

• NAS NETWORK DATA/CODE EXCHANGES 

• TRAINING APPLICATIONS RESEARCHERS 


LeRC 



CFD FIVE-YEAR PLAN 


THRUST 

YEAR 

GOAL 

89 90 91 92 93 

CFD 

DEVELOPMENT 


PREDICTIVE 
CAPABILITY FOR 
COMPLEX 3D 
VISCOUS FLOWS 
FOR ADVANCED 
A/C MISSILES, & 
PROP. SYS. 

IMPROVE SPEED/ACCURACY N-S STOKES 


COMPLEX GEOMETRY/GRIDS 


PARALLEL PROCESSOR ALGORITHMS 


TRANSITION & TURBULENCE MODELS 


HIGH-SPEED REACTING FLOWS 


FLOW PHYSICS 


HIGHLY DETAILED 
DATA BASE TO 
EXTRACT UNDER- 
STANDING OF 
PHYSICS OF 2D & 
3D FLOWS 

LOW SPEED-TRANSITION, TURB., & SEPARATION 


HIGH SPEED TRANSITION & TURB. 


TURBULENCE/CHEM. KINETICS INTERACTION 





3-BLOCK GRID TOPOLOGY 

Nearfield view F-18 ; 300,000 points total 





EFFECT OF PATCHING ALGORITHM 

Moo = .60 a = 20 ° R~c = . 8x 10 6 


Spatial-flux conservation 




F-18 SURFACE FLOW 




NAVIER-STOKES (CFL3D) 


ENTROPY AND VORTICITY EFFECTS IMPROVE ACCURACY OF 
UNSTEADY TRANSONIC SMALL-DISTURBANCE (TSD) THEORY 

• UNMODIFIED THEORY GIVES • ENTROPY AND VORTICITY CORRECTIONS 

INACCURATE SHOCK PREDICTION YIELD EULER-LIKE RESULTS 



I |V 

x/c x/c 



EFFECTS OF REDUCED FREQUENCY ON FIRST HARMONIC 
COMPONENTS OF UNSTEADY PRESSURES DUE 
TO AIRFOIL PITCHING 

• NACA 0012 airfoil at M~ = 0.755, ot 0 = 0.016°, on = 2.51° 





UNSTRUCTURED ADAPTIVE GRID FOR INVISCID SOLUTION 
FOR 3-ELEMENT AIRFOIL (MAVRIPLIS, ICASE) 




FINITE-VOLUME EULER SOLUTION FOR 3-ELEMENT AIRFOIL (MAVRIPLIS, ICASE) 





UNSTRUCTURED GRID FOR AIRFOIL WITH SLAT. HIGHLY-STRETCHED TRIANGLES TO RESOLVE LAMINAR BOUNDARY LAYER. 

(MAVRIPLIS, ICASE) 
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Computational Fluid Dynamics at the Lewis Research Center 

An Overview 

By 

Robert M. Stubbs 
NASA-Lewis Research Center 
Cleveland, Ohio 44135 

Lewis is a multidisciplinary Center with strong research and 
development programs in aeronautical and space propulsion, power, 
space communications, space experiments and materials. 
Computational fluid dynamics (CFD) is playing an important and 
growing role in most of these areas. This presentation describes 
how CFD is integrated into these programs and highlights elements 
of the CFD activities. Examples are presented of codes developed to 
predict flow fields in advanced propulsion systems and several of 
the code validation experiments are described. As will be evident in 
the several Lewis authored papers to be presented at this 
conference, the CFD effort at Lewis ranges from basic research on 
new and improved algorithms through code development to the 
application of these codes to specific engineering problems. 

Because of the substantial improvement in CFD’s predictive 
capability its use at Lewis is on a steep growth path, spreading 
rapidly into new areas which had not traditionally taken advantage 
of the techniques of numerical simulation. The presentation 
concludes with a discussion of multidisciplinary codes and the 
future direction of CFD at Lewis. 
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LEWIS COMPUTER RESOURCES 


CRAY X-MP/24 

SCIENTIFIC VAX CLUSTER 

2 AMDAHLS (VM/CMS AND MVS/XA) 

ALLIANT FX/8 PARALLEL PROCESSOR 

ADVANCED WORK STATIONS (SILICON GRAPHICS, SUN, 
APOLLO) 


T1 LINK TO NAS 




OBJECTIVE 

0 TO DEVELOP ADVANCED COMPUTATIONAL METHODS REQUIRED FOR THE SOLUTION 
OF INTERNAL FLUID MECHANICS AND STRUCTURAL PROBLEMS 

ROLE 


0 GENERATE NEW IDEAS/APPROACHES TO PROPULSION RESEARCH THROUGH 
INTERACTION WITH LeRC STAFF 

0 PROVIDE OPPORTUNITY TO ICOMP PERSONNEL TO DO ORIGINAL RESEARCH 
UTILIZING WORLD-CLASS COMPUTERS 



CFD'S ROLE IN THE PROBLEM SOLVING PROCESS 
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COMPUTATIONAL FLUID DYNAMICS AT LEWIS 

! s f -1 ? :• 

;• - 

. ' i * : 

WHERE DOES CFD PLAY A ROLE? ' 

• AERO PROPULSION 

t 

• SPACE PROPULSION 

• OTHER 

- SPACE POWER 

- MATERIALS PROCESSING 

: . \ f< , • * 

«!,.»■ . • ; 

- FLUIDS IN MICROGRAVITY 

f t 



Notional Aeronautics and 
Space Administration 

Lewis Research Center 


INTERNAL FLUID MECHANICS DIVISION 




CONCEPTUAL ADVANCED TURBOFAN ENGINE 
PHYSICAL PHENOMENA REQUIRING ANALYSES 




/ 


T 


3D VISCOUS FLOWS 
SECONDARY FLOWS 
i nws 

FLOW SEPARATION 
SHOCK/BOUNDARY LAYER 


30 INVISCIDI 
FLOWS 
SHOCKS 






UNSTEADY LOADING 
WAKE EFFECTS 
MULTISTAGE EFFECTS 


SPRAY SIZE/DISTRIBUTIONI 
IGNITION/EXTINCTION 
TURBULENT BURNING 
CHEMISTRY 
JET INTERACTIONS 
RADIATION 


HEAT TRANSFER 
FILM COOLING 




RVC3D (ROTOR VISCOUS CODE 3-m 

BY DR. R. V. CHIMA 

NAVIER-STOKES ANALYSIS FOR 3D FLOWS IN TURBOMACHINERY 

STACKED C-TYPE GRIDS FOR AXIAL OR CENTRIFUGAL MACHINES 

THIN-LAYER NAVIER-STOKES FORMULATION RETAINS HUB-TO-TIP 
AND BLADE-TO-BLADE VISCOUS TERMS > 

; ' * ■ r‘ ' 1 ’ 

BALDWIN-LOMAX TURBULENCE MODEL 

EXPLICIT 4-STAGE RUNGE-KUTTA TIME-MARCHING SCHEME WITH 
VARIABLE TIME STEP AND IMPLICIT RESIDUAL SMOOTHING 

HIGHLY VECTORIZED FOR CRAY X-MP 



RPLUS 3D 


BY DRS. S.T. YU, J. S. SHUEN & P. TSAI 


3D NAVIER-STOKES CODE FOR CHEMICALLY REACTING FLOWS 

FINITE RATE CHEMISTRY MODEL 

O 9 SPECIES INVOLVING O, H, AND N 
O 18 REACTION STEPS 

TIME-MARCHING LU SCHEME 


FAST AND ROBUST 



PROTEUS 


IMPLICIT 2D AND 3D NAVIER-STOKES CODE 

READILY MODIFIED, WELL DOCUMENTED GENERAL SOLVER 

FULL, COMPRESSIBLE NAVIER-STOKES WITH ENERGY EQUATION 

EULER AND THIN-LAYER OPTIONS 

LINEARIZED BLOCK IMPLICIT (LBI) SOLVER 

IMPLICIT BOUNDARY CONDITIONS 

FIRST OR SECOND ORDER IN TIME 

ALGEBRAIC (BALDWIN-LOMAX) AND 2-EQUATIONTURBULENCE 
MODELS CURRENTLY BEING IMPLEMENTED 



MULTISTAGE AVERAGE PASSAGE CODE 

BY DR. JOHN ADAMCZYK 

MULTISTAGE SIMULATIONS FOR ARBITRARY GEOMETRIES 
3D VISCOUS 

EXPLOITS MACROTASKING 

FIRST SIMULATION OF A 3D MULTISTAGE TURBINE 

COMPATIBLE WITH THE DESIGN ENVIRONMENT 

CODE/METHODOLOGY CURRENTLY BEING USED IN DESIGN OF 
TURBOMACHINERY 



Notional Aeronautics and 
Space Administration 

Lewis Research Center 


INTERNAL FLUID MECHANICS DIVISION NASA 


PUTTING IT ALL TOGETHER 


on 


INLETS. DUCTS. 
AND NOZZLES 





CHEMICAL 

TURBOMACHINERY 


REACTING 

FLOWS 


COMPUTATIONAL 
AND EXPERIMENTAL 
TECHNOLOGY 


INTEGRATED MULTIDISCIPLINARY ANALYSIS AND TEST 










LEWIS RESEARCH CENTER 


NUMERICAL SIMULATION OF NONLINEAR DEVELOPMENT OF INSTABILITY WAVES 
TIME DEPENDENT VISCOUS INCOMPRESSIBLE NAVIER-STOKES EQUATIONS 
PROGRESS TOWARD THE DEVELOPMENT OF AN AIRFOIL ICING ANALYSIS CAPABILITY 

THE BREAKUP OF TRAILING-LINE VORTICES 

THREE DIMENSIONAL SIMULATION OF SUPERSONIC REACTING FLOWS WITH FINITE RATE 
CHEMISTRY 

SIMULATION OF TURBOMACHINERY FLOWS 

AUTOMATED DESIGN OF CONTROLLED DIFFUSION BLADES 

NUMERICAL ANALYSIS OF FLOW THROUGH OSCILLATING CASCADE SECTIONS 

NUMERICAL ANALYSIS OF THREE-DIMENSIONAL VISCOUS INTERNAL FLOWS 

A NUMERICAL STUDY OF THE HOT GAS ENVIRONMENT AROUND A STOVL AIRCRAFT 
IN GROUND PROXIMITY 

CFD ANALYSIS FOR HIGH SPEED INLETS 

FLUX SPLITTING ALGORITHMS FOR TWO-DIMENSIONAL REAL GAS FLOWS 


R. MANKBADI 
J. GOODRICH 

M. POTAPCZUK 

C. Bl DWELL 

B. BERKOWITZ 

D. JACQMIN 

S. YU 

J. SHUEN 
P. TSAI 

J. ADAMCZYK 

J. SANZ 

D. HUFF 

R. CHIMA 
J.YOKOTA 

T. VAN OVERBEKE 
J. HOLDEMAN 

T. BENSON 

J. SHUEN 
M. LIOU 



FUTURE DIRECTIONS 


SEVERAL FACTORS CONTRIBUTING TO A CFD EXPLOSION 
AT LEWIS 

- IMPROVED PREDICTIVE CAPABILITY 

- INCREASE IN MACHINE CAPACITY AND AVAILABILITY 

CFD NO LONGER LOCALIZED, NOW IN WIDE USE IN SEVERAL 
DISCIPLINE AREAS 

GROWING EMPHASIS ON MULTIDISCIPLINARY SYSTEM 
SIMULATIONS 


SESSION II 


CENTER OVERVIEWS 
(Continued) 


Chairman: 

Paul Kutler 

Chief, Fluid Dynamics Division 
NASA Ames Research Center 
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Marshall Space Plight Center CPD Overview 
L. A. Schutzenhofer , NASA/MSPC 


Computational Fluid Dynamics (CPD) activities at Marshall Space 
Plight Center (MSPC) have been focused on hardware specific and 
research applications with strong emphasis upon benchmark 
validation. The purpose of this overview is to provide insight 
into the MSPC CPD related goals, objectives, current hardware 
related CPD activities, propulsion CPD research efforts and 
validation program, future near-term CFD hardware related 
programs, and CPD expectations. The current hardware programs 
where CFD has been successfully applied are the Space Shuttle 
Main Engines (SSME), Alternate Turbopump Development (ATD), and 
Aeroassist Plight Experiment (AFE) . For the- future near-term 
CFD hardware related activities, plans are being developed that 
address the implementation of CPD into the early design stages 
of the Space Transportation Main Engine (STME), Space 
Transportation Booster Engine (STBE), and the Environmental 
Control And Life Support System (ECLSS) for the Space Station. 
Finally, CFD expectations in the design environment will be 
delineated . 
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MARSHALL SPACE FLIGHT CENTER CFD 

OVERVIEW 


• COMPUTATIONAL FLUID DYNAMICS BRANCH ACTIVITIES 

- OBJECTIVES 

- INTERACTION 

- APPROACH TOWARD SOLUTIONS OF COMPLEX FLOWS 

• MSFC HARDWARE RELATED ACTIVITIES 

- INHOUSE 

- ROCKETDYNE - SSME 

- PRATT AND WHITNEY -ATD 

• NASA EARTH-TO-ORBIT PROPULSION R&T PROGRAM 

- PROGRAM DEFINITION 

- WORK ELEMENT SUMMARY; CFD EMPHASIS 
“ EXPERIMENTAL APPARATUS 

- CONSORTIUM 

• NEW NEAR TERM CFD ACTIVITIES 

- ADVANCED LAUNCH SYSTEM, ALS 

“ ENVIRONMENTAL CONTROL AND LIFE SUPPORT SYSTEM, ECLSS 

• CFD EXPECTATIONS 
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COMPUTATIONAL FLUID DYNAMICS BRANCH ACTIVITIES 

OBJECTIVES 


• SUPPORT PROGRAM OFFICES 

- "QUICK TURNAROUND" APPLICATIONS 

x 

-INTERACT WITH HARDWARE CONTRACTORS IN DEVELOPMENT OF DESIGN ENVIRONMENTS 
-PROVIDE "SMART BUYER" CAPABILITY FOR LONG-TERM APPLICATIONS 

- DEVELOP SUBSYSTEMS CFD MODELS 

- FOCUS MSFC CFD ACTIVITIES/PROVIDE CENTERWIDE CFD SUPPORT 

• FOCUS DEVELOPMENT OF CFD METHODOLOGY 

- INTERACT WITH ARC, LeRC, LaRC, AND OTHER RESEARCH ORGANIZATIONS TO FOCUS 
TECHNOLOGY DEVELOPMENT TOWARDS MSFC HARDWARE RELATED PROBLEMS 

-DEVELOP REQUIREMENTS FOR CFD CODE VERIFICATION 

-VERIFY CODES THROUGH BENCHMARK COMPARISONS 

-ADVANCE CFD TECHNOLOGY FOR APPLICATIONS 

• DEVELOP ADVANCED HARDWARE TECHNOLOGY CONCEPTS 

-TURBINE STAGE 
-IMPELLER STAGE 
-NOZZLES, PREBURNERS, ETC. 
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COMPUTATIONAL FLUID DYNAMICS BRANCH ACTIVITIES 

ORGANIZATIONAL INTERACTIONS AND UNIQUE 
CFD RESOURCES CAPABILITY 



NASA 

HQ 


DOD 


l 


R&T 

OFF. 


I 


COMPUTATIONAL 
FLUID DYNAMICS 
BRANCH ED32 


I 


PROJECT 

OFF. 


I 


PD 

DIRECT. 



S&E 

DESIGN 

ORG. 


HARDWARE 

CONT. 


SUBSYSTEMS MODELS FOR APPLICATIONS 

- “QUICK TURNAROUND” APPLICATIONS 

- “SMART BUYER” CAPABILITY 

- INTERACT WITH HARDWARE CONTRACTORS 

BENCHMARKED CODES AND SUPPORT 

- FOCUS NATIONAL CFD ACTIVITIES 

- FOCUS MSFC CFD ACTIVITIES 

- REQUIREMENTS FOR CODE VERIFICATION 

- BENCHMARK COMPARISONS 

- TECHNOLOGY TRANSFER IN MSFC 

- ADVANCED CFD TECHNOLOGY 

ADVANCED DESIGN TOOLS 
-AFE 

- TURBINE/IMPELLER STAGES 

- BORE/JOINT FLOW MODELS 
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COMPUTATIONAL FLUID DYNAMICS BRANCH ACTIVITIES 


BASIC RESEARCH 

ALGORITHM DEVELOPMENT 
NUMERICAL METHODS 
GRID GENERATION 

ADAPTIVE GRIDS 
CODE DEVELOPMENT 

FLOW PROCESS MODELING 

ADVANCED COMPUTER SYSTEMS 

ARC, LeRC, LaRC ^ 


CFD CROSS FERTILIZATION 

APPLIED CFD 

— LIBRARY OF CODES 

- BENCHMARK VALIDATION 

- DESIGN CODES FROM 
VARIOUS SOURCES 

— DEVELOP DESIGN CODES 

- DEVELOP CRITERIA TO 
ASSESS CODES 

— EVALUATE ADVANCED HARDWARE 
TECHNOLOGY CONCEPTS 

^ MSFC ^ 


ENGINEERING ANALYSIS 

- PARAMETRIC STUDIES 

- HARDWARE OPTIMIZATION 

- PERFORMANCE 
OPTIMIZATION 

- ANOMALY INVESTIGATION 

- SYSTEM DESIGN 

- DESIGN ASSESSMENT 

MSFC 


5 - 1010 - 8-16 



COMPUTATIONAL FLUID DYNAMICS BRANCH ACTIVITIES 


APPROACH TOWARD SOLUTIONS OF COMPLEX 

FLOWS 


• DEVELOP DATA BASE 

- LITERATURE SEARCH 

- COLLATE RELEVANT EXPERIMENTAL RESULTS 

- IDENTIFY SIGNIFICANT PARAMETERS, SCALING LAWS 

- DEFINE REQUIREMENTS FOR BENCHMARK EXPERIMENTS 

• PERFORM FUNDAMENTAL ENGINEERING ANALYSIS 

- SIMPLIFIED 1 -D OR 2-D ANALYSES 

- ELEMENTARY SENSITIVITY ANALYSIS 

• PERFORM CFD CALCULATIONS 

- EXERCISE BENCHMARKED STATE-OF-THE-ART CODES 

- IMPLEMENT STATE-OF-THE-ART FLOW PROCESS MODELS 



PREDICTION OF SECONDARY FLOW IN CURVED DUCTS OF 

SQUARE CROSS-SECTION 


OBJECTIVE: 


APPROACH: 


COMPUTER 

RESOURCES 

REQUIRED: 

IMPACT: 


TO VERIFY NAVIER-STOKES CODES FOR THE ACCURATE 
PREDICTION OF SECONDARY FLOW IN CURVED DUCTS 


APPLICATIONS OF A 3-D INCOMPRESSIBLE NAVIER-STOKES 
CODE (INS3D) TO FLOW IN A 90° AND 180° BEND AND A 22.5° 
S-DUCT OF SQUARE CROSS-SECTION (25000, 50000, 100000, 
AND 175000 GRID POINTS); COMPARE PREDICTIONS TO LDV 
MEASUREMENTS 


2 TO 6 MW STORAGE ON A CRAY X-MP; 3/4 HOUR RUN TIME 
ON CRAY X-MP FOR THE 1 0 5 GRID POINT CASES 


VERIFICATION OF INCOMPRESSIBLE NAVIER-STOKES CODE 
FOR THE PREDICTION OF SECONDARY FLOWS IN COMPLEX 
INTERNAL FLOW FIELDS 
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TURBOPUMP DESIGN EQUATIONS 
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SSME ► • HPFTP TURBINE BLADES 

• SINGLE CRYSTAL HOLLOW CORE TURBINE BLADES 

• TURBINE DISK CAVITIES 

► • LOX PUMP BEARING INLET CAVITY 

• LOX PUMP BEARINGS 

• FUEL PREBURNER 

• LOX PREBURNER 

• LOX MANIFOLD TEE (4000 Hz) 

• HOT GAS MANIFOLD/MANIFOLD STRUTS 

• PUMP COOLANT FLOW PATHS 

• NOZZLE/M CC MISMATCH 

► • HPOTP NOZZLE PLUG TRAJECTORIES 

• TRANSIENT BEHAVIOR OF FUEL PREBURNER MANIFOLD 

• UTRC HPFTP COOLANT FLOW EXPERIMENT 

• BEARING DEFLECTOMETER (TTBE) 

. ATD • TURBINE INLET TEMP. REDISTRIBUTION 

* • TURBINE TEMP. PROFILE REDISTRIBUTION 

• ROTOR-STATOR INTERACTION 

• TURNAROUND DUCT AND HOT GAS MANIFOLD 

• BEARING ANALYSIS 

• LOX PUMP INLET SCROLL 

• FUEL PUMP INTERSTAGE CROSSOVER DUCTS 

• FUEL PUMP INLET SCROLL 

• LOX PUMP DISCHARGE VOLUTE 

• SEALS 


SRB • BORE FLOW 

—CANTED NOZZLE 
—BROKEN INHIBITOR 

► • FIELD JOINT 

—FLOW AND THERMAL TRANSIENT 
—PRESSURIZATION TRANSIENT 

► • NOZZLE-TO-CASE JOINT 

—FLOW AND THERMAL TRANSIENT 
—PRESSURIZATION TRANSIENT 


ACTIVITIES 


:tivitie< 


INHQUSE 

X 


X 

X 


X 

X 

X 

X 

X 

X 

X 


X 


X 


CONT. 

X 


X 


X 

X 

X 

X 

X 

X 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 



X X 

X X 
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MSFC HARDWARE 


AFE 


U\ 


• AEROTHERMAL ENVIRONMENTS 
► — DSMC 
— NS 


SPACE STATION • CONTAMINATION 

• ECLSS 


ADVANCED PROGRAM • EXTERNAL TANK GAMMA RAY IMAGING 
DEVELOPMENT TELESCOPE 


ACTIVITIES 


1NHQUSE CONT. 


X 

X X 


X 

X X 


X 
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MSFC HARDWARE RELATED ACTIVITIES 
ROCKETDYNE APPLICATION OF CFD TO SSME 


TURBOMACHINERY 


HPFTP IMPELLER CAVITY 
HPOTP P/B BEARING DISCH CAVITY 
HPOTP TURBINE END BRG DISCH CAVITY 
HPOTP 2ND STG TURBINE NOZZLE 
HPFTP 1ST STG TURBINE DISK CAVITY 
HPOTP TURBINE END BRG FLOOD COOL 
HPOTP R/S INTERACTION 
HPFTP TURBINE DISK CAVITIES 
HPFTP 2ND STG TURBINE DISK CAVITY/DIVERTER 
SC/HC 2ND STG TURBINE DISK CAVITY DIVERTER 
^ ROUGH SURFACE SEAL FLOW 
* SC/HC 1 ST STG TURBINE 
HPFTP 1ST STG TURBINE 
SC/HC 1ST STG TURBINE 
HPFTP 1ST STG TURBINE 
SC/HC 1ST STG TURBINE R/S 
HPFTP 1ST STG TURBINE R/S 
HPOTP 1ST STG TURBINE R/S 
LPOTP 4TH STG TURBINE R/S 
HPOTP IMPELLER-DIFFUSER UNSTEADY 
HPOTP BRG FLOWFIELD 
HPFTP 1 ST STG TURBINE W/STRUTS 
HPOTP IMPELLER 

HPOTP 1ST STAGE TURBINE W/STRUTS 
HPOTP P/B DIFFUSER 
T/P HARDWARE DEVIATION SENSITIVITIES 
TO REDUCE MR’S 


CODE 

NONTURBOMACHINERY 

CODE 

STEP 

HOT GAS MANIFOLD 

INS3D 

STEP 

FUEL-SIDE 


STEP 

2-DUCT TURBULENT 


REACT 2D 

3-DUCT TURBULENT 


STEP 

OXIDIZER-SIDE 


STEP 

AXISYMMETRIC 


ROTOR 

3D 


REACT 2D 
STEP 

COMBINED HGM 

INS3D 

STEP 

MAIN INJECTOR 

REACT3D 

REACT 2D 

FLUCTUATING PRESSURE & 


REACT 2D 

DYNAMIC LOADING 


REACT 2D 

REACT 3D 

NOZZLE 

USA 

REACT 3D 

MCC/NOZZLE MISMATCH 

ROTOR 


USA 

ROTOR 

NOZZLE TRANSIENT AND 


ROTOR 

OPERATION 


ROTOR 


USA 

ROTOR 

REACT 3D 

4KHZ RESONANCE 

STEP 

REACT 3D 

REACT 3D 

TEST BED LOX FLOWMETER 


REACT 3D 
REACT 2D/3D 
REACT 2D/3D 
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MSFC HARDWARE RELATED ACTIVITIES 
PRATT AND WHITNEY APPLICATION OF CFD TO ATP 


• MOTIVATION 

- DESIGN VERIFICATION 

• CFD ANALYSIS OF CRITICAL FLOWPATHS IN SSME TURBOPUMPS 

• IDENTIFY POTENTIAL FLOWFIELD NONUNIFORMITIES, REGIONS OF 
SEPARATED FLOWS 

• PROVIDE DETAILED FLOW DATA TO MECHANICAL DESIGN GROUPS FOR 
ADDITIONAL STRUCTURAL, THERMAL ANALYSES 

— ANALYTICAL SUPPORT 

• WHERE NECESSARY, PERFORM FUNDAMENTAL CFD RESEARCH TO 
SUPPORT THE DESIGN VERIFICATION PROCESS 

• INCORPORATE THESE IMPROVEMENTS INTO THE DESIGN DECKS 
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MSFC HARDWARE RELATED ACTIVITIES 
PRATT AND WHITNEY APPLICATION OF CFD TO ATP 


• ACCOMPLISHMENTS 

- ANALYZED COMPLETE HOT GAS FLOW PATH 

- PROVIDED TURBINE INLET CONDITIONS 

- PROVIDED STRUT PRESSURE LOADS FOR MECH DESIGN 

- TAD CALCULATIONS RESULTED IN SUBSTANTIAL REDUCTION IN 
INSTRUMENTATION 

- CFD MODELS OF TAD-HGM, LOX PUMP INLET, PREBURNER AND FUEL PUMP 
CROSSOVER DUCTS READY FOR HOT TEST SUPPORT 

- TURBULENCE MODEL SELECTED FOR COMPLEX DUCT FLOWS 

- ARICC SUBSTANTIALLY UPGRADED 

- CAD TO CFD CAPABILITY IMPLEMENTED 

-2D INVISCID ROTOR-STATOR INTERACTION CAPABILITY DEVELOPED 
AND DEMONSTRATED 
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NASA EARTH-TO-ORBIT PROPULSION R&T PROGRAM 

PROGRAM DEFINITION 


•TECHNOLOGY ACQUISITION PHASE 

- SEEKS IMPROVED UNDERSTANDING OF THE BASIC CHEMICAL AND PHYSICAL PROCESSES 
OF PROPULSION 

- DEVELOPS ANALYSIS METHODS, DESIGN MODELS, AND CODES USING ANALYTICAL 
TECHNIQUES SUPPORTED BY EMPIRICAL LABORATORY DATA AS REQUIRED 

- RESULTS ARE OBTAINED THROUGH TEN DISCIPLINE WORKING GROUPS 


• BEARINGS 

• STRUCTURAL DYNAMICS 

• TURBOMACHINERY V 

• FATIGUE/FRACTURE/LIFE 

• IGNITION/COMBUSTION V 


• FLUID & GAS DYNAMICS V 

• INSTRUMENTATION 

• CONTROLS 

• MANUFACTURING/PRODUCIBILITY/INSPECTION 

• MATERIALS 
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NASA EARTH-TO-ORBIT PROPULSION R&T PROGRAM 

PROGRAM DEFINITION 



• LARGE SCALE SUBSYSTEM TECHNOLOGY VALIDATION 

- VALIDATES TECHNOLOGY EMANATING FROM THE ACQUISITION 
PHASE AT THE LARGE SCALE COMPONENT OR SUBSYSTEM LEVEL 

“ THREE CATEGORIES OF EFFORT 

• LARGE SCALE COMBUSTORS V 

• LARGE SCALE TURBOMACHINERY V 

• CONTROLS AND HEALTH MONITORING 

• TECHNOLOGY TEST BED VALIDATION 

- VALIDATES TECHNOLOGY EMANATING FROM THE ACQUISITION 
PHASE AT THE ENGINE SYSTEM LEVEL 

- THREE CATEGORIES OF EFFORT 

• COMBUSTORS V 

• TURBOMACHINERY V 

• CONTROLS AND HEALTH MONITORING 
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NASA EARTH-TO-ORBIT PROPULSION R&T PROGRAM 
WORK ELEMENT SUMMARY 

• TECHNOLOGY ACQUISITION 

G2 TURBINE DRIVE COMBUSTOR DESIGN 

G29 TTB COMBUSTION MODELS 

G31 COMBUSTION CODE ENHANCEMENTS 

G32 COMBUSTION STABILITY CODE 

G33 TURBULENCE MODELS FOR COMB. ANALYSIS 

G39 ERE PREDICTION METHODS 

H6 FLUID STRUCTURE INTERACTION 

HI 6 VERIFICATION OF INTERNAL FLOW ANALYSIS IN 3D GEOMETRIES 

HI 9 EVALUATION CRITERIA FOR INTERNAL FLOW CFD NUMERICAL MODELING 

H22 ADAPTIVE COMPUTATIONAL METHOD FOR HIGH REYNOLDS NUMBER INTERNAL 

FLOWS IN ADVANCED PROPULSION SYSTEMS 
H23 DEVELOPMENT OF CONVERGENCE ACCELERATION TECHNIQUES FOR 

ALGORITHMS APPLIED TO COMPLEX 3D INTERNAL FLOWS 
H35 ADVANCED INS3D CFD CODE 

H36 CFD CONSORTIUM 

• LARGE SCALE SUBSYSTEM TECHNOLOGY VALIDATION 

LSVT1 EXP. VER. OF CFD TURB. STAGE DESIGN 

LSVT4 HI PRESS TURBOMACHINERY SYS. VALIDATION 
LSVT5 3D TURBOPUMP FLOWFIELD 

LSVT6 EXP. VER. OF IMPELLER STAGE DESIGN 

LSVT1 0 MEASUREMENTS IN MULTI ELEMENT INJECTOR 

LSVT1 2 CFD TURNAROUND DUCT DESIGN VALIDATION 


NASA EARTH-TO-ORBIT PROPULSION R&T PROGRAM 
WORK ELEMENT SUMMARY (CONTINUED) 


•TECHNOLOGY TEST BED VALIDATION 

-TBVC4 INJECTOR DIAGNOSTICS 

-TBVC1 IMPROVED HPOTP PREBURNER PUMP 

-TBVT2 ENHANCED ROTOR CODES 

-TBVT3 IMPROVED BEARING COOLANT PATH 

-TBVT5 WATER FLOW MODELS 

-TBVT8 HGM FUEL SIDE ANALYSIS 

-TBVT9 CFD DATA REDUCTION HARDWARE 

-TBVT1 0 HPOTP JET COOLANT RING 

-TBVT1 3 PREBURNER DOME FILLING FLOW ANALYSIS 

-TBVT24 TURBINE STAGE CFD ANALYSIS AND DATA BASE FOR 

UNSTEADY AERO/HEAT TRANSFER 

-TBVT25 DEV. OF UNSTEADY AERO HEAT/TRANSFER EXPERIMENTS 

DATA BASE FOR AXIAL TURBINE STAGES 
-TBVT26 ADVANCED AXIAL TURBINE STAGE DESIGN METHODS 

-TBVT27 ADVANCED IMPELLER DESIGN METHODS 

-TBVT28 CFD ANALYSIS OF BSMT 

-TBVT29 UTRC ROTOR/STATOR HEAT/TRANSFER 



NASA EARTH-TO-ORBIT PROPULSION 
R&T PROGRAM 

CONSORTIUM OBJECTIVES 


• FOCUS CFD APPLICATIONS IN PROPULSION 

• TECHNOLOGY ACQUISITION PHASE 

•DIRECT BASELINE PROGRAM TOWARDS IMPROVED ACCURACY, 

STABILITY, AND EFFICIENCY 

• LARGE SCALE SUBSYSTEM TECHNOLOGY VALIDATION 

•STIMULATE CFD VALIDATION TOWARDS PROPULSION FLOWS 

• DIRECT APPLICATIONS CODES TOWARD DESIGN TOOLS AND 
ADVANCED HARDWARE TECHNOLOGY CONCEPTS 

• IDENTIFY NATIONAL CFD PROPULSION REQUIREMENTS 

• STIMULATE A FORUM FOR GOVERNMENT, INDUSTRY, AND UNIVERSITY 
INTERACTIONS 

• ENCOURAGE INDUSTRY TO PARTICIPATE IN CFD DEVELOPMENT WITH IRAD FUNDS 

• PROVIDE SYNERGISM IN THE CFD COMMUNITY 

• PROVIDE PEER REVIEW OF CFD PROGRAMS 
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NASA EARTH-TO-ORBIT PROPULSION 
R&T PROGRAM 

CONSORTIUM TASKS 


• DEVELOP A PLAN TO APPLY CFD TO CURRENT AND FUTURE PROPULSION SYSTEMS 

• IDENTIFY AND RANK CRITICAL FLOW PROBLEMS RELATED TO 
PROPULSION SYSTEMS 

• IDENTIFY NATIONAL CFD RELATED RESOURCES 

• DEFINE HIGH PERFORMANCE COMPUTING REQUIREMENTS TO 
ACCOMPLISH CFD FOR PROPULSION APPLICATIONS 

• DIRECT CFD TECHNOLOGY DEVELOPMENT TO PROPULSION APPLICATIONS 

• ASSESS AND VALIDATE CFD APPLICATIONS IN PROPULSION SYSTEMS 

• DEVELOP EVALUATION CRITERIA 

• DEFINE AND IMPLEMENT BENCHMARK VALIDATION 

• DEFINE AND IMPLEMENT VALIDATION TESTS 

• DIRECT THE APPLICATION OF CFD DESIGN TOOLS TOWARDS ADVANCED 
HARDWARE TECHNOLOGY CONCEPTS 

• ACCELERATE THE TRANSFER OF CFD TECHNOLOGY FROM UNIVERSITIES AND 
RESEARCH CENTERS TO INDUSTRY AND HARDWARE DEVELOPMENT CENTERS 
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NASA EARTH-TO-ORBIT PROPULSION R&T PROGRAM 

CONSORTIUM ORGANIZATIONAL STRUCTURE 


00 
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NASA EARTH-TO-ORBIT PROPULSION R&T PROGRAM 

CONSORTIUM TEAMING 


FOCUS DEVELOPMENT OF CFD METHODOLOGY 

AND 

DEVELOP ADVANCED HARDWARE TECHNOLOGY CONCEPTS 



ACTIVITIES SUPPORTED BY ETO RESOURCES 
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CY88 


VALIDATION OF TURBINE STAGE 
DESIGN TOOLS 


PERFORMANCE 


CY89 


CY90 


CY91 


CY92 


CY93 


I I i 

H35 ADVANCED ROTOR 1 & ROTOR 2 CODES 


DEVELOP CODE AT ARC 


I I 

Ele MULTISTAGE TURB. FLOW SIM 





1 30 VISCOUS CODE 1 

r i 

H4m TURB. MODELING 

[ i i 

FOR PROTEUS TURBO. APP. 

| | ANALYTICAL MOUELINC m 



H4a&b 1 



o 



I 


SPECIAL 
TURBINE STAGE 
CFD ANALYSIS 




MSFC INHOUSE TURBINE STAGE BENCHMARK CFD ANALYSES 


I I I I 

LSVT 5 3D TURBOPUMP FLOW FIELD ANALYSIS 


ANALYSIS I CODE ENHANCEMENT 


It I t I t I 

TBVT 24. 25 TURBINE STAGES CFD UNSTEADY FLOW P.B 


AIRFLOW AT TBD CONTRACTOR 





CODE MODIFICATION 


I I I I 

TBVT 26 ADVANCED AXIAL TURB. DES. METHODOLOGY 


ADVANCED CONCEPTS 


PREL. CONCEPT 




!■ 


3D STATOR/ROTOR BLAD 
FLOW EXP. 


LSVT4 


CONCEPT 

DEFINITION 



ACCEPT. 

TEST 


LSVTI EXP. VER OF CFD TURB. STAGE DES. METHOD 


TEST RIG 
DESIGN & FAB 


TESTING 
AT MSFC EATTE 


LSVT 23 IMP AERO H/T PREDICTION 


DEV. D.B. VIA EXP. AT CALSPAN 




o o 


H36 (SUPPORT FROM H19, TBVT 7 & 9) l— F 


CFD CONSORTIUM-VALIDATION OF CFD ACTI VITIES-NATION WIDE 
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NASA EARTH-TO-ORBIT PROPULSION R&T PROGRAMS 
DELIVERABLE PRODUCTS/MILESTONES 
TURBINE STAGES 

• THREE-DIMENSIONAL MULTISTAGE CFD CODES 

• RESEARCH CODE 

- 3-D MULTISTAGE CFD CODE TO PREDICT STEADY AND 1/90 

UNSTEADY FLOW FIELD CHARACTERISTICS, PERFORMANCE, 

LOADS AND HEAT TRANSFER 

• PRODUCTION CODE 

- MODIFICATION TO RESEARCH CODE TO ENHANCE ENGINEERING 8/89 
APPLICATION 

• IMPROVED EFFICIENCIES 8/89 

• STREAMLINE PRE/POST PROCESSING ETC. 

• UNSTEADY THREE-DIMENSIONAL DATA BASE FOR MULTISTAGE TURBINE 

• INITIAL UNSTEADY AERO DATA BASE 6/89 

• ENHANCED UNSTEADY AERO DATA BASE 6/90 

• HEAT TRANSFER DATA BASE 1/90 

• IMPROVED FLOW PROCESS MODELING 

• TURBULENCE MODEL FOR PROTEUS 6/90 

• TURBULENCE MODEL FOR AXIAL TURBOMACHINERY 6/89 

• ADVANCED CONCEPTS AND DEMONSTRATION OF DESIGN TOOLS 

• PRELIMINARY CONCEPT DEFINITION 9/89 

• TESTS 1/89, 6/89, 10/89 

• ADVANCED CONCEPT DEFINITION 1/91 

• FINAL RIG TEST VERIFICATION 1/92 

• HOT FIRE TEST (TTVF) 4/94 


5 - 4588 - 8-53 



NEW NEAR TERM CFD ACTIVITIES 
STB E/STM E DESIGN APPROACH FOR LOW COST 


THE ALS FAMILY 



• DESIGN PRIORITIES 

- RELIABILITY 
-COST 

- PERFORMANCE/WEIGHT 

- COMMONALITY 


• DESIGN BENEFITS 

- REDUCED INTERNAL ENVIRONMENTS 

- ROBUSTNESS 

-REDUCED DEVELOPMENT TIME 

- REDUCED INVENTORIES/INTERCHANGEABILITY 


5 - 1607 - 9-99 



NEW NEAR TERM CFD ACTIVITIES 
STME BASELINE DESIGN REQUIREMENTS 



• GAS GENERATOR CYCLE, SERIES TURBINE DRIVE, 

• LOX/LH2 PROPELLANTS 
CHAMBER PRESSURE = 2250 psi 

• FIXED THRUST OF 580K (VAC) 

• DUAL THRUST: 580K AND 435K (VAC) 

• RELIABILITY = .99, 90% CONFIDENCE LEVEL (DEMONSTRATED) 

• EXPENDABLE OR REUSABLE (15 CYCLES) 

• GIMBAL CAPABILITY FOR TVC, +/-6 DEGREES 

• FIXED NOZZLE, AR = 62:1 

• USABLE IN SINGLE OR MULTI-ENGINE ARRANGEMENT 

• HIGH RELIABILITY, LOW COST 


5 - 1606 - 9-99 



NEW NEAR TERM CFD ACITVITIES 


CFD ACTIVITIES TO SUPPORT STB E/STM E DESIGN 


• THRUST CHAMBER 

-INJECTOR 

-MAIN COMBUSTION CHAMBER 
-NOZZLE 

-COOLING CHANNELS 

• GAS GENERATOR 
: “INJECTOR 

- COMBUSTION CHAMBER 

• PUMPS 

- INLET FLANGE 

- VOLUTE/INDUCER/IMPELLER 

- DIFFUSER/CROSSOVER DUCTS 

- DISCHARGE COLLETOR/DUCTS 
-BEARINGS 

-SEALS 


• TURBINE 

- INLET FLANGE 

- INLET MANIFOLD 

- ROTOR-STATOR INTERACTION 

- MULTISTAGE ANALYSIS 
-AIRFAILS/GUIDE VANES 

- TURBINE EXHAUST — TURN AROUND DUCT 

• SYSTEM ANALYSIS 

- DUCTS 
-MANIFOLDS 
-VALVES 
-CAVITIES 

• ENGINE AIRFRAME INTERACTION 

- PLUME 

-AEROHEATING/LOADS 


5 - 1601 - 9-99 



OF POOR QUALITY 


NEW NEAR TERM CFD ACTIVITIES 
ENVIRONMENTAL CONTROL AND LIFE SUPPORT SYSTEM (ECLSSI 


SPACE STATION CONFIGURATION 


RESPIRABLE ATMOSPHERE AND WATER REQUIREMENTS 



PARAMETER 

UNITS 

OPERATIONAL 

DEGRADED 

CO a PARTIAL PRESSURE 

N/m 2 

400 MAX 

1013 MAX 


(mmHg) 

(3.0 MAX) 

(7.6 MAX) 

0 2 PARTIAL PRESSURE 

N/m 2 x 10® 

19.5-23.1 

16.5-23.7 


(psia) 

(2.83 - 3.35) 

(2.4-3.45) 

TOTAL PRESSURE 

NAn 2 x 10® 

99.9-102.7 

99.9-102.7 


(P8la) 

(14.5-14.9) 

(14.5-14.9) 

TEMPERATURE 

°K 

291.5 - 299.8 

291.5 - 299.8 



(65 - 80) 

(65 - 80) 

DEWPOINT 

•K 

277.6 - 288.8 

273.9 - 294.3 


<T) 

(40-60) 

(35 -70) 

VENTILATION 

msec 

.076 - .203 

.051 - .508 


(ft/mln) 

(15-40) 

(10-100) 

TRACE CONTAMINANTS 

mg/m 3 

TBD 

TBD 


AIRLOCK 


PARTICULATES 


PARTICLES/ 3,530,000 0.5 MICRO- 
NS METERS 

TBD > 1 50 M CROMETERS 


MICRO-ORGANSISMS CFU/m® 


1000 1000 


EMERGENCY 

1600 MAX 
(12 MAX) 

15.8-23.7 

(2.3-3.45) 

99.9-102.7 

(14.5-14.9) 

288.8-302.6 
(60 - 85) 

273.9 - 294.3 
(35 -70) 

.051-1.016 
(10 - 200 ) 

TBD 


1000 


5-1604-9-99 










NEW NEAR TERM CFD ACITVITIES 


CFD ACTIVITIES TO SUPPORT ECLSS 

^ •» , , 

.. T . 

• GENERIC BASELINE CFD MODELS 

- PLANAR, ONE MODULE (NS) 

-3D, ONE MODULE (NS) 

-3D, INNER LOCK DUCTS (NS) 

- PLANAR INTERMODULE (NS) 

-3D, INTERMODULE (NS) 

•, FLOW CONTROL DESIGN PARAMETRIC OPTIMIZATION 

- INTERNAL CONFIGURATION VARIATIONS 
i , “ VENTILATION CONTROL 

- INTRAMODULE VENTILATION (FANS) 

- CONTAMINATION TRANSPORT 
' ~ CO, /FLOW MANAGEMENT 

-BODY FORCE EFFECTS 

• BENCHMARK COMPARISONS 

- PARAMETRIC DESIGN OF EXPERIMENTS 
-CODE VALIDATION 



CFD EXPECTATIONS 


• DIRECT HARDWARE DESIGN UTILIZING CFD 

— PROVIDE INITIAL IMPACT IN DESIGN 

— PERFORM DESIGN OPTIMIZATION STUDIES 

— DEVELOP ADVANCED HARDWARE TECHNOLOGY CONCEPTS 

• ESTABLISH EVALUATION CRITERIA FOR CODES AND CLASSES OF PROBLEMS 

• BENCHMARKED/VALIDATED CODES 

— LAMINAR FLOWS 

— TURBULENT FLOWS u j, p 

— ACOUSTIC PROBLEMS 

— CERTAIN CLASS OF UNSTEADY PROBLEMS 

• USER FRIENDLY CODES 

— B.S. LEVEL ENGINEER 2-3 YRS EXPERIENCE 

— GUIDELINES FOR CLASSES OF PROBLEMS 

— CAD/CAM/CAE; GEOMETRY GRID GENERATION 

— GENERALIZED BOUNDARY CONDITIONS 

— ALGORITHM/GRID OPTIMIZATION FOR SOLUTION EFFICIENCY 

— ARTIFICIAL INTELLIGENCE/EXPERT SYSTEMS 

— MODULAR CODES 

• FLOW ADAPTIVE GRIDS FOR CURRENT CLASS OF PROBLEMS 

• MULTIPLE SCALE AND/OR ZONAL TURBULENCE MODELS, MULTIPHASE, MULTISPECIES, 
COMBUSTION FLOW PROCESS ENGINEERING MODELS EVOLVED FROM EXPERIMENTS 
AND CFD ANALYSIS 


5 - 2005 - 8-43 



V 91-10844 


Johnson Space Center CFD Overview 


C. P. Li 

Advanced Programs Office 
Johnson Space Center, Houston, Tx 77058 


Recent applications and development of CFD technology have focused on flow 
problems that are critically important to the operation and design of space flight 
vehicles. The main effort is spent on the Space Shuttle in order to provide an un- 
derstanding of the cryogenic fluid in the duct connecting the External Tank and 
the Main Engines, the subsonic flow surrounding the Orbiter during crew egress 
maneuvers, the transonic aerodynamic forces on the the Orbiter fuselage and wing, 
the high angle-of-attack abort flight, and the aerodynamic heating during entry. 
To provide in-depth analyses for such diverse problems within a timely schedule, 
matured panel codes and a state-of-the-art incompressible turbulent flow code were 
adapted. Collaboration with Ames Research Center has resulted in a Shuttle ascent 
aerodynamic code; and a viscous chemical nonequilibrium code is being developed 
for predicting Orbiter real-gas aerodynamics and finite-catalytic heating. The re- 
maining activities are devoted to the prediction of the flow environment around the 
Aeroassist Flight Experiment vehicle at hypersonic speeds and high altitudes. A 
thermochemical nonequilibrium Navier-Stokes code has been developed on the ba- 
sis of two- temperature and 11-species models for solving both the shock layer and 
near wake. After validating the code against wind-tunnel aerodynamic, pressure 
and heating data, the code is being used to supplement the ground test facilities 
in predicting a more realistic flight environment. CFD technology is being relied 
upon by other programs as well in the consideration of candidate configurations. A 
biconic cone entering the Martian atmosphere at moderate angles of attack will be 
analyzed for its stability and heating distribution for the proposed mission. Capa- 
bilities of simulating the, low and medium lift-to-drag vehicles flowfiekL flying back 
from the Space Station have been demonstrated and will be enhanced to include 
winglets. The development of hypersonic CFD technology at JSC will continously 
emphasize the modeling of radiation and ablation in continuum flow regime, suffi- 
cient realism of geometry, and efficiency of computational methods. 
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CFD SUPPORT FOR VARIOUS PROGRAMS 


I 


• Shuttle (MY6.5, OSF) 

• Orbiter (MY2.5, OSF) 

• Aeroassist Flight Experiment (MY3.5, OAST) 

• Mars Rover Sample Return Mission (MY1, SE) 

• Crew Rescue Escape Vehicle (MY0.25, SS) 

• High-Energy Aerobraking (MY0.5, OAST) 



Simulation Codes and Computers 

• PAN AIR, VS AERO, QUADPLAN (potential flow, panel method) 

• INS3D (incompressible Navier-Stokes code) 

• F3D (compressible flow NS code) 

• EAGLE, SVTGD3D (grid generation codes) 

• U3D (upwind finite-volume implicit method) 

• NOSIP, AFTB (shock-fitting NS and Parabolized NS codes) 

• VRFNS (shock-fitting, chemically reactive NS codes) 

• VRFLO (shock-fitting, thermochemical nonequilibrium NS 
code) 

• VAX8650s, SCS40 or CX200 (code development) 

• Cray XMP at MSFC and CRAY 2 at NAS (engineering ap- 
plication) 


• Class IV computer (CY89) 



OUTLINE OF PRESENTATION 


• Cryogenic Duct Flow Simulation by Kandula and Pearce 

• Subsonic Orbiter Flow Computation by Slotnick 

• Shuttle Debris Analyses by Gomez, Labbe, Martin 

• Supersonic Orbiter Flow Computation and Comparison by 
Wey and Ma 

• Hypersonic Orbiter Viscous and Reactive Flow Simulation by 
Li 

• Biconic Cone Flow by Stuart 

• AFE results by Gomez, McGary, Tam and Li 

• CERV results 



ANALYSIS OF ET/ORBITER DISCONNECT VALVES 

• Objective: 

- to predict and correlate the hydrodynamic stability of the 
flappers and pressure drop with test data 

• Approach: 

- Adapted and modfied INS3D, SVTGD3D and INGRID 

- Compared with water test data 




<•* 0-0 0.5 1.0 1.5 2.0 2.5 J.t 


GRID IN THE DUCT WITH FLAPPER VALVES 



A PORTION OF X-Z ¥IEW OF THE GRID AT THE PLANE OF SYMMETRY 




Crew Egress Aerodynamic Analysis 


o Objective: Numerically simulate the external flowfield surrounding 
the various flight test vehicles so that an accurate assessment of 
the astronaut exit trajectories may be determined. 

o Methodology: Use production panel code methods to compute aero 
characteristics. 

• Fast, efficient CFD tool 

• Appropriate for subsonic unobstructed external flow 
o Simulations: 

• Space Shuttle Orbiter (VSAERO/PANAIR) 

• Convair C240 (QUADPAN) - Tractor Rocket Concept 

• Lockheed C-141B (QUADPAN) - Pole Concept 


Jeff P. Slotnick / March 8, 1989 
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Johnson Space Center - Houston, Texas 



SPACE SHUTTLE LAUNCH 
VEHICLE (SSLV) CFD ANALYSIS 
STATUS & PLANS 


Advanced Programs Office 


Steven G. Labbe 


3/3/89 


NUMERICAL SIMULATION OF SSLV ASCENT AERO ENVIRONMENT 

- RECENT ADVANCES IN CFD ENABLE FRESH LOOK AT SPACE SHUTTLE ASCENT AERODYNAMICS 

- JOINT EFFORT BETWEEN ARC TEAM (J.L. STEGER, P.G. BUNING) & JSC TEAM (F. W. MARTIN) 

TEAM OBJECTIVES 

- ARC: Develop the Technology to Numerically Simulate Complex Launch Vehicle Geometry 

- JSC: Employ CFD Technology to Gain Insight & Understanding of the Ascent Aero Loads Environment 

CFD ANALYSIS 

- Complex SSLV Ascent Configuration Modelled Using "CHIMERA " Composite Grid Discretization Approach 

- Overset Body-Conforming Grids Of Major Geometry Components - Communication by PEGASUS Code 

- F3D_Implicit Approx. Factored Finite- Difference Procedure - Solution of 3-D Thin-Layer NS Equs. 

CFD ANALYSIS PROGRESS HAS BEEN STEADY & PROMISING 

- EFFORTS CONCENTRATED ON FULLY UPDATING SSLV CONFIGURATION GEOMETRY DETAILS 

Highest Priority Given to Attach Hardware & LOX/Fuel Feedlines 

SRB/IEA Attach Ring and SRB Plume Simulations Incorporated -- Check Out Proceeding 

- ABILITY TO ACCURATELY DETERMIN WING LOADS-* ► MODEL SSLV GEOMETRY DETAILS 







AERO DATABASE COMPARISONS 


MID-FUSELAGE 


WING UPPER SURFACE 
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STS-27 Ascent Debris Trajectory Simulation 


Mach 

Alpha 



MSA— 1 (SRB ablator) released <§ 0° — 90° on nose cone 


ORBITER AERODYNAMICS AT HIGH INCIDENCE ANGLES 


• Objective: 

- To predict aerodynamics and flight characteristics for an- 
gles of attack up to 90 deg 

• Approach: 

- Validate the U3D code with wind-tunnel data for angles 
of attack lower than 50 deg 

- Apply the multi-zonal U3D for higher angles-of- attack flow 



ORIGINAL PAGE IS 
OF POOR QUALITY 


ORBITER GRID AND MACH 5.2 COMPARISON 



CFD & WIND TUNNEL COMPARISONS 


a - 007 
4 - 1000 



MACH 7.32 FLOW RESULTS AND GRID 



1 5 

£r o 



o 


MACH NUMBER 


CONTOUR LEVEL 

1 

0.164 

2 

0.361 

3 

0.558 

4 

0.755 

5 

0.953 

6 

1.150 

7 

1.347 

6 

1.544 

0 

1.742 

10 

1.939 

11 

2.136 

12 

2.333 

13 

2.530 

14 

2.728 

15 

2.925 

16 

3.122 

17 

3.319 

18 

3.517 A 

19 

3.714 1 

20 

3.911 

21 

4.108 

22 

4.306 

23 

4.503 

24 

4.700 

25 

4897 

26 

5.095 

27 

5.292 

28 

5489 

29 

5.686 

30 

5.884 

31 

6.081 

32 

6.278 

33 

6475 

34 

6.672 

35 

6.870 



WEY 1 / 13 / 8 ® 


Orbiter Alone 

Surface Mach No 

Moo 

7.32 ' 

a 

iroo 

Grid 

7fl x 35 x 51 " j 



ORIGINAL PAGE IS 
OF POOR QUALITY 


NUMERICAL INVESTIGATION OF CONTINGENCY ABORT 


o 

00 



T. C. Way/LESC 
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ORBITER ENTRY FLOW SIMULATION 


• Objective: 

- To assess aerodynamic and heating issues for each flight 

- To validate the VRFNS code using flight data 

• Approach: 

-Use a chemical nonequilibrium model and shock-fitting 
technique 

- Use a decouple technique for species 

- Compare the perfect-gas model results with wind-tunnel 
data 



110 


INVISCID CHEMICAL NONEQUILIBRIUM FLOW 

Orbiter canopy 



BOW SHOCK SURROUNDING A SWEPT- WING VEHICLE 

AT MACH 22 and AOA 40 WITH EQUILIBRIUM CHEMISTRY 



O.o- 10.0 20.0 0.0 0.0 10.0 20.0 30.0 10.0 


PROJECTED VIEWS OF TEMPERATURE CONTOURS 



Comparison of shock shape with data 


MRSR aerocapture 
CFD 


Mach 5.997 
Alpha 27 deg 
Air 

Navier-Stokes Solution 
45 x 1 9 x 29 grid 
1 200 iterations 
CFL 2 


Wind tunnel 
shock locations 



Flow 



Mach Number Contours 


12.84°/7° Biconic 
Mach 6.0 Alpha 27° 



Grid Pressure Temperature 


115 


AFE FLOWFIELD SIMULATION 


• Objectives: 

- To predict aerodynamic loads and heating distributions 
for flight conditions 

To define flowfield environment using appropriate physical 
models 

• Approaches: 

- Developed the VRFLO code for typical aerobraking vehi- 
cles 

- Calibrated the code using wind-tunnel data 

Compared the two-temp and 11-species model results with 
RAMC data 

- Performed sensitivity studies of different models 
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PHYSICAL MODELING FOR HYPERSONIC VEHICLES 




Comparison of Shock Shapes and Subsonic Regions 

for an AFE Configuration 



Ideal Air 
Mach 9.81 



Ideal CF^ Equilibrium Air 

Mach 6.29 Mach 32.0 



250a 300 


COMPARISON OF FLIGHT AND WIND-TUNNEL PRESSURES 



(a) Front view (b ) Eear view 

Comparison of flight and wind-tunnel wall 
pressure distributions (the left side corresponds to the 
flight case). 


CONTOUR PLOTS OF AFE FLIGHT FLOWFIELD 

inviscid chemical nonequilibrium flow in the pitchplane 



(b) Isotherm T/T t » 


ORIGINAL PAGE IS 

OF POOR QUALITY 


WALL PRESSURES COMPARISON AND HEATING 

LOW L/D CERV 
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TECHNOLOGY DEVELOPMEMT 


• A multi-zonal Navier-Stokes code to predict the Orbiter flow 
at incidence up to 90 deg 

• A thermochemical nonequilibrium PNS code for high L/D 
configurations 

• An Euler code with adaptive, unstructured grid for complex 
geometry 

• Radiation and ablation modeling and code implementation 
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NASA'S CFD VALIDATION PROGRAM 
by 

Dale R. Satran 
Program Manager 
Aerodynamics Division 
NASA Headquarters 

With computational fluid dynamics (CFD) becoming a productive research 
and design tool, the requirement to validate CFD codes has grown 
significantly. NASA has emphasized CFD validation activities since 1986 
when a separate work element was formed to fund experimental activities 
related to validation. NASA's CFD and CFD validation programs are 
closely coordinated to ensure that experimental data bases are available as 
soon as possible for validating codes. In response to industry and academic 
requirements, four levels of experimental research have been defined as 
part of CFD validation with NASA's Aeronautics Advisory Committee 
(AAC) support although only the fourth level actually has the detailed 
information necessary for validating codes. 

Critical flow physics especially turbulence modeling are key to improved 
CFD codes. NASA has focused additional resources on transition and 
turbulence physics to meet these requirements. With improved turbulence 
models, CFD codes will be more accurate, robust, and efficient. However, 
with the level of detailed information available from CFD codes, highly 
accurate and detailed experiments are required to capture the critical 
information for validating codes. Advanced instrumentation especially 
non-intrusive instrumentation is required to acquire this information in 
validation experiments. The CFD validation program is being coordinated 
and managed to address these critical activities. A list of experiments 
which are currently being supported at least partially has been included 
with this presentation. 


PRECEDING PAGE BLANK NOT FILMED 
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CFD CODE VALIDATION DEFINITION 


• DETAILED SURFACE-AND-FLOW-FIELD COMPARISIONS WITH EXPERIMENTAL 
DATA TO VERIFY THE CODE’S ABILITY T O ACCURATELY MODEL THE 
CRITICAL PHYSICS OF THE FLOW. VALIDATION CAN OCCUR ONLY WHEN THE 
ACCURACY AND LIMITATIONS OF THE EXPERIMENTAL DATA ARE KNOWN AND 
THROUGHLY UNDERSTOOD AND WHEN THE ACCURACY AND LIMITATIONS OF 
THE CODE’S NUMERICAL ALGORITHMS, GRID-DENSITY EFFECTS, AND PHYSICAL 
BASIS ARE EQUALLY KNOWN AND UNDERSTOOD OVER A RANGE OF SPECIFIED 
PARAMETERS. 
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CFD VALIDATION CATEGORIES 


USSSSL 


CATEGORIES OF CFD-RELATED EXPERIMENTATION 

A. EXPERIMENTS DESIGNED TO UNDERSTAND FLOW PHYSICS 

B. EXPERIMENTS DESIGNED TO DEVELOP PHYSICAL MODELS FOR CFD CODES 

C. EXPERIMENTS DESIGNED TO CALIBRATE CFD CODES 

D. EXPERIMENTS DESIGNED TO VALIDATE CFD CODES 


ALL FOUR CATEGORIES ARE IMPORTANT AND ARE NECESSARY 


TO BUILD A MATURE CFD CAPABILITY 
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IMPLEMENTATION PLAN 


• ALL EXPERIMENTS HAVE BEEN CLASSIFIED AND DOCUMENTED 

- GOALS 

- LIMITATIONS 

- MODELING 

■ PARTICIPATION 

- LEVEL OF EFFORT 

• SEVERAL KEY EXPERIMENTS INVOLVE MULTIPLE RESEARCH CENTERS 

• CFD VALIDATION WORKSHOP HELD TO IDENTIFY CRITICAL NEEDS 

• COORDINATING BOARD FOR CFD VALIDATION DEVELOPING UPDATED 
DETAILED IMPLEMENTATION PLAN 

• EFFORTS INITIATED TO INVOLVE THE AEROSPACE INDUSTRY AND 
UNIVERSITIES 
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CFD VALIDATION PROGRAM 



EXPERIMENTS HAVE BEEN CLASSIFIED 
INTO MULTIPLE CATEGORIES 


CATEGORY 

AMES 

LANGLEY 

LEWIS 

TOTAL 

A. FLOW PHYSICS 

16 

35 

29 

80 

B. FLOW MODELING 

8 

7 

13 

28 

C. CODE CALIBRATION 

6 

19 

12 

37 

D. CODE VALIDATION 

7 

24 

16 

47 






TOTAL NUMBER 
OF EXPERIMENTS 

27 

45 

29 

101 
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CFD VALIDATION PROGRAM 


• EXPERIMENTS COVER LARGE SPEED RANGE 


- SUBSONIC: 

- TRANSONIC: 

- SUPERSONIC: 

- HYPERSONIC: 


33 EXPERIMENTS 
27 EXPERIMENTS 
23 EXPERIMENTS 
18 EXPERIMENTS 


• EXPERIMENTS FALL INTO SEVERAL VEHICLE CLASSES 


- GENERIC 

- FIGHTER/ATTACK 

- SUBSONIC TRANSPORT 

- ROTORCRAFT 

- ASTOVL 

- PROPULSION SYSTEMS 
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CFD VALIDATION EVENTS 


- ■■■■ ■■■ 

• NEW RTOP ELEMENT, FLOW MODELING AND VERIFICATION, 
CREATED 

• NRC ASEB REVIEW OF CFD ACTIVITIES 

• NASA REVIEW AND DEVELOPMENT OF IMPLEMENTATION PLAN 
FOR CFD VALIDATION 

• AAC AD HOC SUBCOMMITTEE REVIEW OF CFD VALIDATION 

• NASA COORDINATING BOARD FOR CFD VALIDATION FORMED 

• FIRST NASA CFD VALIDATION WORKSHOP AT AMES 

• IMPLEMENTATION PLAN REVISED BY COORDINATING BOARD 

• AGARD CFD VALIDATION CONFERENCE IN LISBON 

• CFD VALIDATION ACTIVITIES AND IMPLEMENTATION PLAN 
REVIEW 

• NASA CFD CONFERENCE AT AMES 

• SECOND NASA CFD VALIDATION WORKSHOP 


FY 1986 

FY 1986 
FEB., 1986 

FY 1987 
JUNE, 1987 
JULY, 1987 
AUG., 1987 
MAY, 1988 
NOV., 1988 

MAR., 1989 
JULY, 1990 



130 


NASA CFD VALIDATION WORKSHOP 



• 103 PERSONS ATTENDED FROM NASA, DOD, INDUSTRY, AND UNIVERSITIES 

• 31 PRESENTATIONS WERE GIVEN ON CFD VALIDATION STATUS 

• 6 WORKING GROUP SESSIONS FOCUSED ON NEAR AND FAR TERM NEEDS 

• NUMEROUS RECOMMENDATIONS 

/ 

- STANDARDIZED TEST CASES FOR CALIBRATION 

- CLOSE COOPERATION BETWEEN CFD DEVELOPERS AND EXPERIMENTALISTS 

- INCREASE FLIGHT-BASED ACTIVITIES 

- DETAILED MEASUREMENTS OF FLOW FIELD AND BOUNDARY CONDITIONS 

- IMPROVED OR NEW NON-INTRUSIVE MEASUREMENT CAPABILITIES 
■ REDUNDANCY IN BOTH MEASUREMENTS AND EXPERIMENTS 
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SUMMARY 



• AAC AD HOC TASK TEAM RECOMMENDATIONS 
IMPLEMENTED 

• NASA PROGRAM EXPANDING TO COVER 
ADDITIONAL AREAS 

• INSTRUMENTATION HAS BEEN ADDED TO 
SEVERAL FACILITIES BY VALIDATION PROGRAM 

• INCREASED EMPHASIS ON COOPERATIVE PROGRAMS 
WITH UNIVERSITIES AND INDUSTRY 


AMES CFD VALIDATION PROGRAM FOR FY 1989 

LEX/DELTA VORTICAL FLOW 
TRANSONIC LOW ASPECT RATIO WING-BODY 
REARWARD FACING STEP 
SSME TURNAROUND DUCT 

SUPERSONIC SHOCK BOUNDARY LAYER INTERACTION 
COMPRESSIBLE PRESSURE-DRIVEN 3-D INTERACTIONS 

2- D TRANSONIC CIRCULATION CONTROL 

3- D SPIN FLOWS 

3-D LOW SPEED WEDGE FLOW WITH SEPARATION 

TRANSONIC SUPERCRITICAL AIRFOIL 

LOW SPEED HIGH ALPHA INVESTIGATION 

CFD VALIDATION FOR WING AERODYNAMICS 

3-D HIGH ASPECT RATIO SEPARATED FLOW 

STOVL AERO/PROPULSION INTERACTION 

THERMO-CHEMICAL NONEQUILIBRIUM FLOWS 

PHOTODIAGNOSTIC INSTRUMENTATION 

UNSTEADY VISCOUS FLOW 

HYPERSONIC SHOCK BOUNDARY INTERACTION 

TURBULENT SHEAR LAYERS 

TURBULENT BOUNDARY LAYERS 

ALL-BODY HYPERSONIC TEST 

HIGH SPEED ROTOR FLOWS 

HYPERSONIC REAL GAS 

SHOCK TUNNEL NOZZLE TESTS 

3.5' HWT NOZZLE TESTS 

COMBUSTION/DETONATION 

FLIGHT/CFD CORRELATION OF F-18 WING PRESSURES AT HIGH ALPHA 
SUPERSONIC VORTEX-SHOCK WAVE INTERACTION 
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LANGLEY CFD VALIDATION PROGRAM FOR FY 1989 

TRANSONIC HIGH ASPECT-RATIO WING 
TRANSONIC LOW ASPECT RATIO WING 
REARWARD FACING STEP IN WATER TUNNEL 
REARWARD FACING STEP IN BART 
DELTA WING VORTEX FLOWS 
SUPERSONIC COAXIAL JET 
TURBULENT MODELING IN SEPARATED FLOWS 
45-DEG SWEEP AIRFOIL 
BARF LDV TEST 

SUPERSONIC BOUNDARY LAYER TRANSITION 

NTF FLAT PLATE TEST 

VORTEX BURST EXPERIMENTS 

HYPERSONIC FLIGHT INSTRUMENTATION 

HYPERSONIC INLET TESTS IN HELIUM 

HYPERSONIC SHOCK-ON-LIP 

HALIS ORBITER EXPERIMENT 

BLUNT BODIES (AOTV/AFE) EXPERIMENT 

HYPERSONIC WINGED SLENDER BODY 

OSCILLATING CANARD/WING UNSTEADY PRESSURES 

VALIDATION OF JET PLUME MODULES 

SUPERSONIC JET PLUME DYNAMICS 

SUPERSONIC HIGH-ALPHA FLOWFIELD 

OFF-AXIS WING-BODY STUDY 

STORE/CAVITY SEPARATION EXPERIMENTS 

WAVERIDER DESIGN PROCEDURE 

5 DEG CONE EXPERIMENT 

75/76-DEG DELTA WINGS 

NTF FOREBODY/MISSLE MODEL 

LEADING EDGE VORTEX FLAP 

X-29 EXPERIMENT IN NTF 

3-D TRANSONIC CAVITY FLOW 

LOW REYNOLDS NUMBER AIRFOIL EXPERIMENTS 

CONFLUENT BOUNDARY LAYER 

GORTLER INSTABILITY ON AIRFOILS 

EXPERIMENTAL INVESTIGATION OF TURBULENCE 

RANGE AND ACCURACY OF THIN FILM ARRAYS 

JUNCTURE FLOW EXPERIMENT 

SWEPT SUPERCRITICAL HLFC AIRFOIL EXPERIMENTS 
TWIN ENGINE AFTERBODY EXPERIMENT 
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LEWIS CFD VALIDATION PROGRAM FOR FY 1989 

3-D SHOCK WAVE/TURBULENT BOUNDARY LAYER INTERACTIONS 

3-D FLOWS IN HIGH SPEED TURBOMACHINERY 

BLADE SURFACE BOUNDARY LAYER 

FUNDAMENTAL SEPARATION BUBBLE RESEARCH 

AIRFOIL (BLADING) FLOW CONTROL 

LEADING EDGE STAGNATION REGION 

BOUNDARY LAYERS IN TRANSITION 

UNSTEADY HEAT TRANSFER IN ROTOR WAKES 

TRANSITION DUCT - AERO & HEAT TRANSFER 

VORTEX GENERATORS 

SHEAR LAYER EXCITATION - JET MIXING 

SHEAR LAYER EXCITATION - SLOT RESONATOR 

MULTI-PHASE FLOWS 

MULTI-PHASE FLOW AND FLUID SPRAY STUDY 
LOW TEMPERATURE HEAT TRANSFER 
FUEL SWIRLER CHARACTERIZATION 

COMBUSTION CHARACTERISTICS OF HYDROCARBON FLAMES 
KINETIC STUDY OF H2/02 SYSTEM 
FLOW INTERACTION EXPERIMENT 
HOT GAS INGESTION 

COHERENT STRUCTURES IN SUPERSONIC SHEAR LAYER 
AERO CHARACTERISTICS OF AIRFOIL WITH ICE ACCRETION 
TURBOMACHINERY BLADE ROW INTERACTIONS 
SUPERSONIC THROUGH-FLOW CASCADE RESEARCH 
CENTRIFUGAL COMPRESSOR FLOW RESEARCH 
SUPERSONIC THROUGH-FLOW FAN RESEARCH 
HIGH REYNOLDS NUMBER (HEAT TRANSFER) 

DETAILED AERO OF ADVANCED TURBOPROPS 
FUEL RICH CATALYTIC COMBUSTION 
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UNDERSTANDING TRANSITION AND TURBULENCE 
THROUGH DIRECT SIMULATIONS 


P. R. Spalart & J. J. Kim 

A. R. C., C. F. D. Branch, Turbulence Physics Section 


Direct simulations consist in solving the full Navier-Stokes equations, 
without any turbulence model, and describing all the detailed features of the 
flow. Usually the flows are three-dimensional and time-dependent and contain 
both coarse and fine structures, which makes the numerical task very chal- 
lenging in terms of both the algorithm and the computational effort. Most 
of the work until now has involved spectral methods, which are highly accu- 
rate but not very flexible in terms of geometry or complex equations. For 
that reason, future work will also rely on high-order finite-difference or other 
methods. 

Direct simulations complement experimental work, and both contribute 
to the theory and the empirical knowledge of turbulence. Once such a sim- 
ulation has been shown to be accurate the flow field is completely known, in 
three dimensions and time, including the pressure, the vorticity and any other 
quantity. On the other hand, most simulations to date solved the incompress- 
ible equations in rather simple geometries, and direct simulations will always 
be limited to moderate Reynolds numbers. Extensive simulations have been 
conducted in homogeneous turbulence, channel flows, boundary layers, and 
mixing layers. Much effort is devoted to addressing flows with compressibil- 
ity and chemical reactions, and to new geometries such as a backward-facing 
step. 
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UNDERSTANDING TRANSITION AND TURBULENCE 

I 

THROUGH DIRECT SIMULATIONS 

P. R. Spalart & J. J. Kim 

A. R. C., C. F. D. Branch, Turbulence Physics Section 

00 

• Secondary instability in channel flow. 

• Stability and relaminarization of the swept attachment-line flow. 

• Improvement of the pressure term in turbulence models. 

• Coherent structures in homogeneous and wall-bounded shear flows. 

• Lyapunov exponents of a turbulent channel flow. 

• Growth of compressible mixing layers. 
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Secondary instability in channel flow 

• Direct simulation is applied to “natural” transition; only small 
(« 10 -10 ) random disturbances are introduced. 

• By the primary instability, a 2D TS wave grows (wavenumbers 
k x = 1, k z = 0) (x streamwise). The final 3D breakdown may 
start as a “K” type (fundamental, k x — 1, aligned A vortices) or as 
an “H” type (subharmonic, k x = 1/2, staggered A’s). Theory and 
experiments don’t fully agree in a channel. 

• We find that span wise nonuniformities ( k x = 0, k z ^ 0) accumulate 
energy much larger than the initial disturbances (by a factor Re) 
and transfer some to the “K” modes before the Herbert secondary 
instability. 

• This explains the difficulties in observing the “H” type breakdown 
in experiments. 
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Stability and relaminarization of the swept attachment-line flow 

• Turbulence has been found to propagate from the fuselage along 
the attachment line and render Laminar Flow Control impossible. 
This is “leading-edge contamination” (Gregory, 1960). 

• Direct simulations are conducted in this region. Curvature, com- 
pressibility, and span wise variations are neglected. 

• Previous linear-stability results (Gortler Sz Hammerlin, 1955, 

Hall, Malik Sz Poll, 1984) are confirmed and generalized. 

• The unswept flow is found to be linearly and nonlinearly stable. 

• The instability and relaminarization boundaries are computed in 
(K, R) space, where K is the suction parameter and R the Reynolds 
number based on strain rate and sweep velocity component. Suction 
has much less effect on relaminarization than on linear stability. 
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Critical Re at Leading Edge 
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Pressure terms in a Reynolds-stress turbulence model 


• The complete budgets of the Reynolds-stress tensor and the 
dissipation tensor have been extracted from a direct simulation of 
turbulent channel flow. 


• This allows one to test a turbulence model term by term, instead 
of as a whole (e.g., using the Cf). 


• The published estimates for the terms that are difficult to measure 
experimentally (pressure, dissipation) were often very inaccurate, es- 
pecially near a wall. 


• As a result, the widespread k — e and Reynolds-stress models incur 
large errors near walls. 
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Pressure 



Kkm kk 27. Pressure 
computed from the 
<m | iaal ion (4h)- 
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Coherent structures in shear flows 


• Streaks have long been observed near the wall in boundary layers 
(y + < 20), and horseshoes away from the wall. These observations 
were confirmed by simulation results. 

• Why the difference? 

• Near the wall, the nondimensional shear rate S* = Sq 2 /e takes 
values much larger than in the homogeneous flows that had been 
studied ( q 2 : turbulent energy; e: dissipation rate). 

• When S* is given such values in a homogeneous shear flow, streaks 
appear, indicating that the presence of a wall is not essential to form 
them. 
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Profile of S* in Boundary Layer 



0 
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The dimension of strange attractors in turbulent shear flows 

• Does a turbulent solution live on a “strange attractor”? (i.e., does 
its state always return to the same region, but never settle down, 
with almost identical initial conditions rapidly moving apart?) 

• If so, chaos (nonlinear dynamics) theory may provide a new frame- 
work to unify known turbulent phenomenology and to understand 
the dynamics rather than the statistics of turbulence. 

• Using numerical simulation of low Reynolds number turbulent 
channel flow we have confirmed the existence of a strange attrac- 
tor by measuring its Lyapunov exponent spectrum and calculating 
its dimension. 

• At a Reynolds number R r (= 8u T /v) of 80 the dimension of the 
attractor is ~ 360. This is the first measurement of the intrinsic 
complexity of a fully developed turbulent flow. 

• This result shows that shear flow turbulence cannot be considered 
to result from the interaction of a “few” degrees of freedom. 
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Direct simulation of compressible (free shear) flows 


• Goal: study compressibility effects on turbulent flows, in particular 
on their structure, global evolution, and aerodynamic noise. 

• Tool: D. N. S. using high-order compact finite differences (close 
to spectral resolution). Compressible ideal-gas Navier- Stokes equa- 
tions. No artificial viscosity or filtering. 

• Example: Spatially evolving mixing layers. 2D, forced at inflow 
witth 1% amplitude. R = Pi{U\ — £72)<Wmi ~ 200 — 400. 

• Experiments have shown that: a) compressibility reduces the growth 
rate; b) it scales with the convective Mach number 
M c = (U\ — U 2 )/(a 1 + <22) (71 = 72). 


• Simulations have: a) reproduced this effect; b) validated M c \ 
c) provided a physical explanation. 
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Direct Simulation of Compressible Turbulence 

T.A. Zang and G. Erlebacher 
NASA Langley Research Center 

M.Y. Hussaini 

Institute for Computer Applications in Science and Engineering 
NASA Langley Research Center 


The physics of turbulence remains one of the most challenging problems in fluid dynam- 
ics. Although more than a century of effort has been devoted to it, a lot of fundamental 
issues are still unresolved. This is particularly so in the case of turbulence in high speed 
flows because of the increased number of possible mode interactions due to compressibility 
effects. For example, the cubic non-linearities in the momentum equations allow the vor- 
ticity, acoustic and entropy modes to interact with each other. The dynamics are further 
complicated by the possible existence of non-stationary shocks and/or eddy shocklets. 

In this paper, several direct simulations of 3-D homogeneous, compressible turbulence 
are presented with emphasis on the differences with incompressible turbulent simulations. 
A fully spectral collocation algorithm, periodic in all directions coupled with a 3rd order 
Runga-Kutta time discretization scheme is sufficient to produce well-resolved flows at Tay- 
lor Reynolds numbers below 40 on grids of 128x128x128. A Helmholtz decomposition of 
velocity is useful to differentiate between the purely compressible effects and those effects 
solely due to vorticity production. In the context of homogeneous flows, this decompo- 
sition is unique. Time-dependent energy and dissipation spectra of the compressible and 
solenoidal velocity components indicate the presence of localized small scale structures. 
These structures are strongly a function of the initial conditions. We concentrate on a 
regime characterized by very small fluctuating Mach numbers Ma (on the order of 0.03) 
and density and temperature fluctuations much greater than Ma 2 . This leads to a state in 
which more than 70% of the kinetic energy is contained in the so-called compressible com- 
ponent of the velocity. Furthermore, these conditions lead to the formation of curved weak 
shocks (or shocklets) which travel at approximately the sound speed across the physical 
domain. Various terms in the vorticity and divergence of velocity production equations are 
plotted versus time to gain some understanding of how small scales are actually formed. 
Possible links with Burger turbulence are examined. 

To visualize better the dynamics of the flow, new graphic visualization techniques have 
been developed. The 3-D structure of the shocks are visualized with the help of volume 
rendering algorithms developed in house. A combination of stereographic projection and 
animation greatly increase the number of visual cues necessary to properly interpret the 
complex flow. The presence or absence of shocks is automatically detected by monitoring 
of the minimum and maximum divergence of the velocity field over the physical domain. 
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NAVIER-STOKES EQUATIONS 

Non-Dimensionalizations 

Velocity, C/ 0 , Temperature T 0 , density po, pressure pUfi. 

Viscosity and Prandl number are constant. 

It + V-(H = 0 


Dv 

P-RI = - 


1 — ♦ 

XrVP + — V-f 


1 Ml 


= _ 7 pv-t?+ 7 V-(/cVr) + 7 ^ 7 f : I. 

Z>£ e p PrRe e p Re 

e p P = pT 

where the stress tensors and dissipation function are defined by 

|=i(Vtr+Vt? r ) ) 

2 i 


t = 2pe- -p{V-v)I, 

?: ?= £(W+ Vtf) : (Vi?+ Vtf) - ^p(V-v) 2 . 
L o 
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OBJECTIVE 


• Understand if, when, and why shocks occur in homogeneous turbulence, shear 
layers and mixing layers in the supersonic and hypersonic regimes. 

• Develop mixing enhancement strategies in supersonic shear layers (e.g. scram- 
jets) 

• The study of initial condition effects on the time-dependent characteristics of 
supersonic “isotropic” turbulence is the first step toward the aforementioned 
objectives. 
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PHYSICS 

• Solve time-dependent, 3-D full Navier-Stokes equations 

• Code set up for variable viscosity, but set to a constant 

• Pr = 0.7 (air) 

• ideal gas 

• no chemistry 

• initial energy and thermodynamics autocorrelation spectrums identical 
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INITIAL CONDITIONS (1) 


• Generate random v,6p,6T subject to 


V-v s = 0 
Vxv c = o 


• Impose autocorrelation spectrum on v,6p,6T 


E(k) = fcV* 2/2 *o 


• Specify rms levels for 6p , ST 

• Calculate 


P = 1 + Sp 
T = 1 + ST 



INITIAL CONDITIONS (2) 


Specify fluctuating Mach number M a , and compute 


Ml = 


«o 


7 RT 0 


according to the approximate formula 


M a = Mq u 


rms 


• Specify initial level of compressibility 

fu 2 c dV 

x ~;(u! + u 2 )dV 



NUMERICAL METHOD 

• Fourier Collocation in both directions (periodic box) 

• Conservative scheme in absence of time discretization errors 

• Isotropic truncation in Fourier space at every iteration 

• Splitting algorithm 

— 3rd order Runga-Kutta method in time on the explicit ter ms (1st step) 

- Acoustic terms treated implicitly (2nd step) 





K.E> dissipation 


turb3d/run 1 32 
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omega**2 


5.30 

5.15 

5.01 

4.86 

4.71 

4.57 

4.42 

4.27 

4.13 

3.98 

3.84 
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CONTOUR LEUELS 
-3. 10000 
-3.00000 
-2.90000 
-2.80000 
-2.70000 
-2.60000 
-2.50000 
-2.400DD 
-2.30000 
- 2.20000 
-2. 10000 
- 2.00000 
-1 . 90000 
-1 . 80000 
-1.70000 
-1.60000 
-1.50000 
-1.40000 
-1 . 30000 
- 1 . 20000 
- 1 . 10000 
- ' . oooo r 


-c, 60000 
-o, 500 OC 

-0.40000 
-0.30000 
- 0.20000 
- 0. 10000 
0.00000 
0. 10000 
0.20000 
0.30000 


OIUERGENCE OF UEL0C1TY 


turb3d/run132, t=3.3, it=600 



1 . 000 fIRCH 

0.00 DEG RLPHR 

-10 Re 

5.5 TIME 

65x65x166 GRID 
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CONCLUSIONS 

• Weak shocks can be generated by an initially random 2-D velocity and ther- 
modynamic field. 

• These shocks propagate at the speed of sound across the domain 

• The presence of these shocks is reflected in the structure of the compressible 
energy spectrum 

• Although the flow is isotropic at the onset, the compressible component of 
velocity quickly becomes anisotropic (as evidenced by the shock structure) 

• Sophisticated visualization techniques are necessary to capture the essence of 
the dynamics in 3 dimensions 


Future Work 

• Mixing layers with shock interactions 

• Three-dimensional isotropic turbulence (in progress) 

• Turbulence modeling Testing 
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Non Linear Evolution of a Second Mode Wave in Supersonic Boundary Layers 

G. Erlebacher 

NASA Langley Research Center 
M.Y. Hussaini 

Institute for Computer Applications in Science and Engineering 
NASA Langley Research Center 


Recent advances in supersonic and hypersonic aerospace technology have led to a renewed 
interest in the stability and transition to turbulence of high speed flows. The last 30 years have 
intermittently witnessed some vigorous attempts to understand some of the fundamental routes 
to transition for incompressible flows. While a fairly comprehensive picture of the initial stages 
leading to the breakdown of an incompressible laminar boundary layer has emerged (mostly under 
controlled conditions) the non-linear effects responsible for transition at high speeds are still very 
much a mystery. However, the current nonlinear incompressible theories, numerical simulations and 
experiments will, hopefully, serve as a guide in gaining a better understanding of the m *r h*ni amp 
present in the supersonic and hypersonic regimes. 

Compressible linear stability theory diverges from incompressible linear stability theory in sev- 
eral ways. Incompressible inviscid instabilities, linked with the existence of an inflection point in 
the mean streamwise velocity profile are replaced by a correlation between inviscid compressible 
instabilities and a generalized inflection point, which brings into play the mean density profile. 
Furthermore, as the Mach number increases, the growth rate of the 3-D modes begin to overtake 
those of the 2-D modes. Beyond Mach 2.2, multiple unstable modes (at fixed Reynolds number 
and frequency) can coexist. The higher modes are inviscid in nature and have different behaviors 
with regard to wall cooling. Not only are they more unstable than their first mode counterparts 
(viscous in nature) above a certain Mach number, but they are destabilized by wall cooling, which 
is detrimental for high altitude hypersonic aircraft. 

It is of vital importance that the nonlinear nature of these second mode instabilities be under- 
stood, and that their role in the context of transition be elucidated. The objective of the work is 
to understand the possible equilibrium state of a second mode wave, before initiating a study of 
3-D wave interactions. 

Two years ago, a spectral code was developed to perform direct simulations of subsonic and 
supersonic flows over flat plates. In this paper, we present several direct simulations of one 2-D 
second mode perturbation wave, superimposed upon a prescribed mean flow. Periodicity is assumed 
in the streamwise direction (Fourier) and the variables are expanded in Chebyshev series in the 
direction normal to the plate. The code is fully explicit and is time advanced with a 3rd order 
Runga-Kutta scheme. The second mode wave (Rf = 8000), interacts with itself to generate higher 
streamwise harmonics. Physical parameters are chosen to maximize the linear growth rate at the 
prescribed Reynolds number. Initial results indicate that the nonlinear processes begin in the 
critical layer region and are the result of the cubic interactions in the momentum equations, rather 
than due to the higher streamwise harmonics. Analysis of the various terms in the momentum 
equations combined with numerical experiments in which various modes are artificially suppressed, 
lead to the conclusion that asymptotic methods will produce the saturated state in one or two 
orders of magnitude less computer time than that required by the direct numerical mm nlat.mnff 
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CURRENT STATUS OF FLAT PLATE 
STABILITY AND TRANSITION 


• Theory 

— Inviscid Stability Theory (Lees and Lin, 1946) 

— Linear Parallel Theory (Mack 1965) 

— Linear Non-Parallel Theory (Nayfeh and El-Hady, 1980) 

— Non-Linear Theories 

* Wave-Interactions: NONE 

* Secondary Instability: (El-Hady 1988, Ng 1988) 

• Experiment 

— Stability of Supersonic Boundary Layers (Laufer and Vrebalovich 1960) 

• Numerical 

— Linear Stability of Ideal Gases (Mack 1975, Malik 1982, Macaraeg Si Streett 1988, 
Erlebacher Si Herbert 1988, Ng 1988) 

— Linear Stability of real gases (equilibrium air) (Malik 1988, Macaraeg Si Streett 1988) 

— Non-Linear Stability of (Erlebacher and Hussaini 1987) 



PHYSICAL PARAMETERS 


Non-Dimensionaiization 


length: S* 

u,p,T: free-stream values 

P- pul 

• boundary layer flow 

• ideal gas (air) 

• Pr = 0.70, Re = 8000 

• M=4.5 

• a = 2.25, xj) = 0° 

• — linear frequency: 2.05 => period = 3.07 

- linear growth: .0215 

- 20 periods in time gives amplification of 3.76 

- e 9 amplification in 136 periods 
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DISCRETIZATION 


• Fully explicit 

• 3 rd order low-storage Runga-Kutta in time 

• Fourier collocation in two periodic directions (stream and span) 

• Chebyshev collocation in vertical direction 

• Zero normal stress boundary-conditions in the free-stream 

• Continuity equation is imposed at the wall and in the far-field 

• Zero temperature perturbations at the wall 



MODAL ANALYSIS 


Let 


u oc e *"* 
u 

— J=S —COJ 

u 

w,- = Real(—) 

u 

oj r = Imag(~) 


OO 


max 


w (^s/)= £ u n (y,t)d 

n=—oo 

After substitution into the x-momentum equation, 
ii = linear + cubic + viscous terms 


where 

Linear = lcxuq uipo + • 
Cubic = taui u-i pi + • 



growth rate 


0.30e— 01 
0.27a— 01 
0.24e— 01 
0.21a— 01 
0.18a— 01 
0.15a— 01 
0.1 2e— 01 
0.90a— 02 
0.60e— 02 
0.30e— 02 
0.00 









2nd mode 



run 19 M=4.5 









2nd mode run19 M=4.5 




0.35e— 01 
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^ O.ISa-01 
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rho ampl (10) 


O.IOo+OO 


0.90a— 01 
0.80a— 01 
0.70a— 01 
0.60a— 01 
0.50a— 01 
0.40a— 01 
0.30a— 01 
0.20a— 01 
0.1 0a— 01 
0.00 





2nd mode run 19 M=4-.5 



0.35e— 02 


0.30e— 02 


0.2Se— 02 

o 

CM 

^ 0.20e— 02 

Q. 

£ 0. 1 5e— 02 
J 
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0.00 


O. 
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vort prod dilatation 


0.50a— 01 
0.40a— 01 
0.30a— 01 
0.20a— 01 
0.1 0a— 01 
0.56a— 08 
— .1 0a— 01 
—.20a— 01 
—.30a— 01 
—.40a— 01 
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vort prod barocl torque 
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0 . 1 2 o — 01 


0.80e— 02 


0 . 40 o — 02 


—. 130—08 


—. 400—02 


—. 800—02 


-.12o—01 


—.I 60 —OI 
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COST OF DIRECT SIMULATION 

• The number of operations for a single iteration is approximately twice that of 
an incompressible direct simulation 

• Acoustic effects must be treated implicitly for low Mach number simulations. 
At high Mach numbers, time step restrictions are a function of the Reynolds 
number 

Steep gradients of v! and T' in the critical layer increase the initial resolution 

Compressibility effects weaken the secondary instability, further increasing the 
required computer time to capture the initial stages of breakdown 

Cost can decrease with use of static or dynamics adaptive grids, or multido- 
main decompositions in the vertical direction 

. VERY EXPENSIVE !!! 

25 periods on 8 x 2 x 65 grid: 10 CPU hrs on Cray II at 
100 Mflops average yields factor 4 amplification 




CONCLUSIONS 

• Compressible boundary-layer code is suitable for the study of second modes, 
albeit is still very expensive due to the slow growth and high spatial and 
temporal frequencies of the instabilities 

• Non-linearities (not due to higher streamwise harmonics) induces strong growth 
rate departures from the linear values in the critical layer region. The actual 
cause is still unknown. 

• Density perturbations probably play (as expected) a fundamental role in the 
development of the non-linear saturated state. Further information awaits 
more detailed considerations of intermodal energy transfers. 
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NUMERICAL SIMULATION OF NONLINEAR DEVELOPMENT OF INSTABILITY WAVES 

Reda R. Mankbadi 

Institute for Computational Mechanics in Propulsion 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135, U.S.A. 

and 

Cairo University 
Cairo, Egypt 


SUMMARY 

The present work is concerned with the nonlinear interactions of high 
amplitude instability waves in turbulent jets. In plane shear layers Riley 
and Metcalf (1980) and Monkewitz (1987) have shown that these interactions are 
dependent, among other parameters, on the phase-difference between the two 
instability waves. Therefore, in the present work we consider the nonlinear 
development of both the amplitudes and the phase of the instability waves. The 
development of these waves are also coupled with the development of the mean 
flow and the background turbulence. In formulating this model it is assumed 
that each of the flow components can be characterized by conservation equations 
supplemented by closure models. Results for the interactions between the two 
instabi 1 ity waves under high-amplitude forcing at fundamental and subharmonic 
frequencies are presented here. Qualitative agreements are found between the 
present predictions and available experimental data. 


CONSERVATION EQUATIONS 

Each flow component is split in the form: 

U.(x,r,t) = U.(x,r) + u.(x,r,t) + u!(x,r,t) 

U is the time-averaged mean flow velocity which is taken to be given by the 
two-stage hyperbolic tangent profile. 0 is the periodic component which is 
split into two frequency-components in the form: 

i 4>, (x)-iw,t ~ i<tu(x)-ico-t 

u^ = A^" u. i'(.r,0) e + A 2 u 2 .(r,e) e + c.c. 

u is the radial shape taken as the eigen-function solution of the local ly- 
parallel linear stability equation for each frequency-component. A and <|> 
are the amplitude and phase to be obtained from the nonlinear interaction equa- 
tions, and 0 is the momentum thickness u' is the turbulence component which 
is related to the turbulence' energy T through an assumed Gaussian profile. 


THE NONLINEAR INTERACTION EQUATIONS 

Time-averaging and phase-averaging techniques are applied to the full 
unsteady Navier-Stokes equations to derive the governing equations for each 
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flow component. These equations are manipulated to obtain nonlinear equations 
for the amplitudes A], A 2 , phases 4 >i , $ 2 . momentum thickness 0 and the tur- 
bulence energy T: 


Mean flow 


Turbulence 


1 ^ma d0_ = _ I T _ x a* I A? 

2 d0 dX mt com 1 2u 2 


BX [ I ta T ] = l ~ J + I “^ A ’ T + I, “ < - A ’ T ' - 


l mt + ut 1 2tot 2 £ 

uwcomponent 

h [ : »a A l] ■ U A 1 - Lt A ? T * 'a-Vl C0S<2 *» • *2“ ■ * C> 

d<t> 

I T-Sf = irS +1* + A, I, sin(2<|> - 4u - + o) 

1 ua dx g) wm 2 2u (a 2u o 

2u-component 

5X [wH ■ I 2»™*2 - W* 1 - ““*« ‘ +2 “ ' 

.2 

▼ - i - / _ _ 

" *2» ‘ *0 


d(b« A, 

*2ua ”d)T = . irS 2u + *2wm “ A^ l Zwa Sin(2< * > u " *2u ^o + a 


The integrals I appearing in the above equations are functions of 0, fre- 
quencies, and the closure assumptions. S is the Strouhal number defined as # 
wd/(2irU) . The solution of the above system of equations is subject to the ini- 
tial conditions at x = 0: 9 q. T 0 , Ai 0 , A 2q ^ and $< 3 . 


RESULTS AND DISCUSSIONS 

The calculated fundamental and subharmonic components at Strouhal numbers 
0.3 and 0.6 are shown in figure 1 for several initial phase angles. The ini- 
tial momentum thickness is 0.026 R, initial turbulence energy 
the initial energies of the fundamental and subharmonic are 
the initial instability velocity components at the jet centerline are 1.2 and 
0.6 percent, respectively. Figure 1 shows that the fundamental is not sensi- 
tive to the phase-difference as much as the subharmonic does. Bradley and Ng s 
(1989) measured integral spectral amplitude shows similar features. The funda- 
mental is less dependent than the subharmonic on the phase angle 
subharmonic amplification occurs at 4>o = ar, d fni n i mum 
fication occurs at 4> 0 = 0 °, same as the present results in figure Kb). 

The calculated centerline phase-averaged velocities are shown in figure 2 
in comparison with the data of Arbey and Ffowcs-Wi 1 1 iams (1984). The Strouhal 


levels is 0.0001 
taken such that 


Maximum 
subharmonic's ampli- 
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numbers are 0.3 and 0.6, the initial turbulence energy level is 0.00001, ini- 
tial momentum thickness is 0.012. The initial centerline velocity of the 
S = 0.3 component is taken 1.5 percent and that of the S = 0.6 is 0.38 of the 
S = 0.3 component as in the experiment. At S = 0.3, figure 2(a) shows that 
calculate peak occurs further down stream as compared to the measured ones. 
However, the calculated peak has the same level as the measured one. The peak 
increases when 4> 0 is changed from 0 to 180°, as the present computations also 
predict. The calculated phase averaged velocities at S = 0.6 shown in fig- 
ure 2(b) has the same features as th'e measured one; same level of amplification 
and same dependency on <J> 0 . The measured component increases again after it 
decays which is probably due to Its interaction with other frequency compo- 
nents. This mechanism is not accounted for here. 

“The dependency of the subharmonic amplification on the initial phase angle 
2s shown in figure 3 for Strodjial numbers 0.2 and 0.4. The initial levels are 
uf 0 = 7 percent, u s0 = 0.5 percent. The peak of the subharmonic at $0 = 270' 
is three times higher than its peak at <p 0 « 90°. The corresponding momentum 
thickness is shown in figure 4, compared to the unexcited momentum thickness. 
The figure shows that the momentum thickness is only weakly dependent on the 
phase angle. This indicates that the direct role of the subharmonic in turbu- 
lent jets in control 1 ing- the mixing is less pronounced as compared to its role 
in controlling the mixing in Laminar jets. However, the subharmonic can still 
have a strong role in the mixing process through enhancing the background tur- 
bulence which in turn increases the mixing. 

If both the fundamental and subharmonic's initial levels are high, the 
dependency ori the phase angle is less pronounced as figure 5 indicates. The 
Strouhal numbers are 0.3 and 0.6 and the initial levels are Uf = u s = 3 per- 
cent. At high initial levels, large energy levels are drained from the mean 
flow and therefore the fundamental-subharmonic energy exchanges are relatively 
smaller and consequently less pronounced. 

The effect of the forcing level at a fundamental frequency of S = 0.4'on 
the subharmonic's amplification is shown in figure 6. The figure shows that 
the peak of the subharmonic increases with increasing the forcing level. How- 
ever, a saturation condition occurs around a forcing level of 10 percent. 

Higher forcing levels result in no further ‘increase of the subharmonic's peak 
over 20 percent. 
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Figure 2 


COMPARISON BETWEEN CALCULATED PHASE-AVERAGED CENTERLINE 
VELOCITIES AND ARBEY & FFOWCS-WILLIAMS’ (1984) DATA 




188 


Figure 3 


DEPENDENCY OF SUBHARMONIC’S AMPLIFICATION ON THE INITIAL 
PHASE-ANGLE AT STROUHAL NUMBERS OF 0.4 AND 0.2 
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Figure 4 


DEVELOPMENT OF THE MOMENTUM THICKNESS UNDER 
TWO-FREQUENCY EXCITATION AND FOR THE UNEXCITED CASF 
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Figure 5 


DEPENDENCY OF THE SUBHARMONIC’S PEAK ON THE PHASE-ANGLE AT 
HIGH LEVELS BOTH THE FUNDAMENTAL AND THE SUBHARMONIC 
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Figure 6 


DEPENDENCY OF THE SUBHARMONIC’S PEAK ON THE 
FORCING LEVEL. AT STROUHAL NUMBERS OF 0.2 AND 0.4 
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More Accurate Predictions with Transonic Navier-Stokes 
Methods Through Improved Turbulence Modeling 

by 

Dennis A. Johnson 

Experimental Fluid Dynamics Branch 
NASA Ames Research Center 


Because the aerodynamic characteristics of aircraft in the transonic 
regime are so sensitive to viscous effects, the selection of the turbulence model 
for a transonic prediction method is no less important than the selection of the 
numerical algorithm. Yet, the usual practice in transonic airfoil, Reynolds- 
Averaged, Navier-Stokes codes has been to employ "equilibrium" algebraic tur- 
bulence models. Satisfactory results are obtained with these turbulence models 
for weak interaction cases (i.e., cases where the upper surface shock wave is too 
weak to have a major effect on the turbulent boundary layer). Such is not the 
situation for cases where the shock wave is sufficiently strong to cause 
separation. The danger in using these "equilibrium" turbulence models for 
airfoil design is that they can result in unduely optimistic projections of aircraft 
performance at off-design conditions. 

Significant improvements in predictive accuracies for off-design con- 
ditions are achievable through better turbulence modeling; and, without neces- 
sarily adding any significant complication to the numerics. One well established 
fact about turbulence is it is slow to respond to changes in the mean strain field. 
With the "equilibrium" algebraic turbulence models no attempt is made to model 
this characteristic and as a consequence these turbulence models exaggerate the 
turbulent boundary layer’s ability to produce turbulent Reynolds shear stresses 
in regions of adverse pressure gradient. As a consequence, too little momentum 
loss within the boundary layer is predicted in the region of the shock wave and 
along the aft part of the airfoil where the surface pressure undergoes further 
increases. 

Recently, a "nonequilibrium" algebraic turbulence model was formulated 
which attempts to capture this important characteristic of turbulence. This 
"nonequilibrium" algebraic model employs an ordinary differential equation to 
model the slow response of the turbulence to changes in local flow conditions. In 
its originial form, there was some question as to whether this "nonequilibrium" 
model performed as well as the "equilibrium" models for weak interaction cases. 
However, this turbulence model has since been further improved wherein it now 
appears that this turbulence model performs at least as well as the "equilibrium" 
models for weak interaction cases and for strong interaction cases represents a 
very significant improvement. 

The performance of this turbulence model relative to popular 
"equilibrium” models is illustrated for three airfoil test cases of the 1987 AIAA 
Viscous Transonic Airfoil Workshop, Reno, Nevada. A form of this 
"nonequilibrium" turbulence model is currently being applied to wing flows for 
which similar improvements in predictive accuray are being realized. 
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TRANSONIC FLOW TURBULENCE MODELING 
AIRFOIL AND WING FLOWS 


WHY BE CONCERNED ABOUT THE TURBULENCE MODEL? 

• NUMERICAL PREDICTIONS CAN BE EXTREMELY 
SENSITIVE TO TURBULENCE MODEL IN TRANSONIC REGIME 


• WIDELY USED EQUILIBRIUM ALGEBRAIC TURBULENCE 
MODELS OVERPREDICT AIRCRAFT PERFORMANCE 



NONEQUILIBRIUM JOHNSON-KING 
ALGEBRAIC TURBULENCE MODEL 


SLOW RESPONSE OF TURBULENT SHEAR STRESS TO 
CHANGES IN STRAIN FIELD MODELED 

ORDINARY DIFFERENTIAL EQUATION FOR MAXIMUM 
REYNOLDS SHEAR STRESS ESTABLISHES EDDY VISCOSITY 
(SIMPLIFIED "REYNOLDS STRESS" MODEL?) 

DUE TO RECENT IMPROVEMENTS, MODEL PERFORMS BETTER 
FOR ATTACHED FLOW CASES (SURPRISINGLY, ORIGINAL 
MODEL WORKED BETTER FOR STRONGLY SEPARATED CASES) 
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SURFACE PRESSURE COEFFICIENT AXISYMMETRIC BUMP 




CLOSURE MODEL 

PRESENT 

JONES-LAUNDER 

CEBECI-SMITH 

Q EXPERIMENT 
S SEPARATION 

R REATTACHMENT 






198 


MAXIMUM REYNOLDS SHEAR STRESS DISTRIBUTIONS 




IMPROVEMENTS IN JOHNSON-KING MODEL 


NEW INNER EDDY VISCOSITY EXPRESSION 

• BETTER SATISFIES "LAW OF THE WALL" FOR ATTACHED 
ADVERSE PRESSURE GRADIENT CONDITIONS (ORIGINAL 
MODEL OVERPREDICTED C f DOWNSTREAM OF SHOCK WAVE) 


• INCLUDES VELOCITY SCALE WHICH ACCOUNTS FOR 
COMPRESSIBILITY 
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IMPROVEMENTS IN JOHNSON - KING TURBULENCE MODEL FOR 
WEAK SHOCK WAVE/BOUNDARY LAYER INTERACTIONS 

RAE 2822, M exp = 0.729, cv exp = 3.19° 



C L C D 
0.803 0.0168 

0.787 0.0159 

0.828 0.0170 

0.825 0.0170 
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ATTACHED AIRFOIL CASE 

RAE 2822, M exp = 0.73, c* exp = 3.19° 



SEPARATED AIRFOIL CASE 

RAE 2822, M a = 0.75, = 3.19' 


BALDWIN-LOMAX "EQUILIBRIUM" MODEL 

□ EXPERIMENT 


iSoooa 


0.837 (0.743) 
0.0312 (0.0242) 



X/C 


.8 
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SEPARATED AIRFOIL CASE 

NACA 0012, M exp = 0.8, c* exp = 2.86° 
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SUMMARY 


SIGNIFICANT IMPROVEMENTS IN PREDICTIVE ACCURACY WITH 
NONEQUILIBRIUM ALGEBRAIC TURBULENCE MODEL 

PRESENT APPROACH IS TO MODEL MOST RELEVANT PHYSICS 
WHILE MAINTAINING MATHEMATICAL SIMPLICITY 

NEGLIGIBLE INCREASE IN COMPUTATIONAL EFFORT COMPARED 
TO EQUILIBRIUM ALGEBRAIC TURBULENCE MODELS 



SESSION IV 


CFD CODES 


Chairman: 

Jerry C. South, Jr. 

Head, Analytical Methods Branch 
NASA Langley Research Center 



precedes page blank ^ot filmed 


RECENT ADVANCES IN RUNGE-KUTTA SCHEMES FOR SOLVING 
3-D NAVIER-STOKES EQUATIONS 

By 

Veer N. Vatsa 

NASA Langley Research Center 
Bruce W. Wedan, Ridha Abid 


Abstract 

A thin-layer Navier-Stokes has been developed for solving high Reynolds number, turbulent flows past aircraft components under tran- 
►o sonic flow conditions. The computer code has been validated through data comparisons for flow past isolated wings, wing-body configu- 
rations, prolate spheroids and wings mounted inside wind-tunnels. The basic code employs an explicit Runge-Kutta time-stepping scheme 
to obtain steady state solution to the unsteady governing equations. Significant gain in the efficiency of the code has been obtained by 
implementing a multigrid acceleration technique to achieve steady-state solutions. The improved efficiency of the code has made it feasi- 
ble to conduct grid-refinement and turbulence model studies in reasonable amount of computer time. The non-equilibrium turbulence 
model of Johnson and King has been extended to three-dimensional flows and excellent agreement with pressure data has been obtained 
for transonic separated flow over a transport type of wing. 
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OBJECTIVES 

• Develop an efficient Navier-Stokes code for high 
Reynolds number, transonic, separated flows 

• Assess the effect of grid refinement 

— on solution accuracy 

— on convergence properties 

• Improve turbulence model for separated flows by 
including non-equilibrium effects 

• Validate the code via data comparisons 
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GOVERNING EQUATIONS 

• Reynolds-averaged Navier-Stokes equations 

• Thin-layer approximation 

• Equations are written in conservation law form 

• Turbulence models 

Equilibrium model : Baldwin-Lomax model 
Non-equilibrium model : Johnson-King model 
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BOUNDARY CONDITIONS 

• Solid surface : viscous 

no slip and zero injection 

zero normal pressure gradient 

specified temperature or adiabatic condition 

• Solid surface : inviscid 

zero flux across surface 
extrapolate surface pressure 

• Inflow/outflow : free-air 

Riemann invariants for farfield 
extrapolate all variables at downstream 

• Inflow/outflow : in-tunnel simulations 

Riemann invariants at inflow 

specify initial profile for viscous sidewall 

specify pressure and extrapolate other 
variables at downstream 



NUMERICAL ALGORITHM 


Based on Runge-Kutta schemes of Jameson 
and co-workers : 5-stage scheme with 3 
dissipation evaluations 

Finite volume, central-difference scheme 

Non-isotropic artificial dissipation added for 
stability 

Variable coefficient, grid aspect-ratio dependent, 
implicit residual smoothing for increasing stability 
bound 

Multigrid acceleration technique 

— Full multigrid ( FMG ) strategy 

— V-cycle ( saw-tooth ) 

— Viscous fluxes evaluated only on fine mesh 



LOG(RESIDUAL) 


Effect of grid refinement on convergence 
history for ONERA M6 wing 

( = 0.84, a = 6.06° ) 

Baldwin-Lomax Turbulence Model 


193x65x33 grid 289x65x49 grid 



WORK WORK 


SUP/SUPF 



213 


Effect of grid refinement on pressure 
distributions for ONERA M6 wing 

( IVL = 0.84, a = 6.06° ) 

Baldwin-Lomax Turbulence Model 

O □ Experimental data 193x65x33 grid results 289x65x49 grid results 




x/c 
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Effect of grid refinement on pressure 
distributions for ONERA M6 wing 

( IVL = 0.84, a = 6.06° ) 
Baldwin-Lomax Turbulence Model 


O □ Experimental data 193x65x33 grid results 289x65x49 grid results 
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Convergence history for ONERA M6 wing 
with Johnson-King model, 289x65x49 grid 

( M„ = 0.84, a = 6.06° ) 



0 400 800 1 200 


WORK 
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Effect of turbulence model on 

pressure contours for ONERA M6 wing 

( = 0.84, a = 6.06° ) 

Upper surface, 289x65x49 grid 

Johnson-King Baldwin-Lomax 
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Effect of Turbulence Model 
ONERA M6 Wing 

M = 0.84, AOA = 6.06 deg, Re = 11.7 million 


Baldwin — Lomax 


Johnson 4 — King 
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Effect of turbulence model on pressure 

distributions for ONERA M6 wing 

( M„ = 0.84, a = 6.06° ) 

289x65x49 grid computations 


O □ Experimental data Baldwin-Lomax results Johnson-King results 



x/c x/c 
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Effect of turbulence model on pressure 
distributions for ONERA M6 wing 

( = 0.84, a = 6.06° ) 

289x65x49 grid computations 


O □ Experimental data Baldwin-Lomax results Johnson-King results 





CONCLUDING REMARKS 


Significant gains in efficiency are achieved 
through multigrid acceleration technique 

Grid-convergence studies feasible due to 
improved efficiency 

Baldwin-Lomax model gives good solutions for 
attached flows, but is found inadequate for 
separated flows 

Johnson-King model results in improved 
comparison with data for separated flows 

Block-structured grids must be employed for 
more efficient use of mesh points and for 
computing more complex configurations 
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COMPUTATIONS OP THREE-DIMENSIONAL STEADY 
AND UNSTEADY VISCOUS INCOMPRESSIBLE FLOWS 


Dochan Kwak, Stuart E. Rogers 
NASA Ames Research Center 

Seokkwan Yoon, Moshe Rosenfeld, and Leon Chang 
MCAT Institute 


The INS3D family of computational fluid dynamics computer codes is presented. These codes 
are used to as tools in developing and assessing algorithms for solving the incompressible Navier- 
Stokes equations for steady-state and unsteady flow problems. This work involves applying the 
codes to real-world problems involving complex three-dimensional geometries. The algorithms 
utilized include the method of pseudo compressibility and a fractional step method. Several ap- 
proaches axe used with the method of pseudocompressibility including both central and upwind 
differencing, several types of artificial dissipation schemes, approximate factorization, and an im- 
plicit line-relaxation scheme. These codes have been validated using a wide range of problems 
including flow over a backward-facing step, driven cavity flow, flow through various type of ducts, 
and steady and unsteady flow over a circular cylinder. Many diverse flow applications have been 
solved using these codes including parts of the Space Shuttle Main Engine, problems in naval 
hydrodynamics, low-speed aerodynamics, and biomedical fluid flows. The presentation details 
several of these including the flow through a Space Shuttle Main Engine inducer, vortex shedding 
behind a circular cylinder, and flow through an artificial heart. 
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OUTLINE 


© Objective and Approach 

© Summary of Flow Codes 

• INS3D Family of codes 

• CENS3D 

© Applications and Results 

• Space Shuttle Main Engine (SSME) components 

• Artificial Heart Flow 

© Summary and Future Work 
© Movie 

• Circular cylinder vortex shedding 

• Artificial heart flow 
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OBJECTIVE AND APPROACH 


© Objective 

• To develop CFD capability for simulating steady and unsteady 
viscous incompressible flows (Incompressible Navier- Stokes) 

© Approach 

• Develop and assess algorithms, and implement in codes 

• Develop / implement physical models for engineering analysis 
(turbulence, cavitation, porous medium, etc.) 

• Apply the codes to real-world problems 
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SUMMARY OF INS3D 


0 Governing equations 

• Incompressible Navier- Stokes equations in generalized 3-D co- 
ordinates for steady-state solutions 

• Pseudocompressibility approach 

0 Numerical scheme 

• Finite difference, central differencing plus artificial dissipation 

• Approximate Factorization 

• Single or multiple zones 

© Turbulence Models 

• Algebraic models 

• k — e model 

© Applications 

• Numerous SSME related simulations 

© Status 

• Distributed to numerous users across the nation 

• Available through COSMIC 
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EXTENSIONS TO INS3D 


0 INS 3D family of research codes used to study various approaches to 
solving the INS equations in generalized 3-D coordinates 


Pseudo compressibility Approach 

© INS3D-UP (Steady-State and Time-accurate calculations) 

• Upwind flux- difference splitting of uniformly high order 

• Line-relaxation implicit scheme 

• Characteristic boundary conditions 


© INS3D-LU (Steady-State and Rotating reference frame) 

• Finite-volume method 

• Spectral radius or flux- difference split based dissipation 

• LU-SGS Implicit Scheme 

• Non-reflecting boundary conditions 

• Completely vectorized 
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EXTENSIONS TO INS3D, continued 


Fractional Step Approach 

0 INS3D-FS (Time- accurate problems) 

• Finite- volume method on a staggered mesh 

• Accurate treatment of geometry 

=> Exact discrete mass conservation 

• Two-step fractional step method : 

Solve momentum equations in time (AF scheme) 
Correct for pressure and velocity (Poisson equation) 



SUMMARY OF CENS3D CODE 


© Governing Equations 

• Compressible Euler and Navier- Stokes equations and species 
transport equations in generalized 3-D coordinates 

© Numerical Methods 

• Fully-coupled and implicit thermal-chemical nonequilibrium 
finite-rate-chemistry 

• Finite volume / flux-limited TVD, optional high-order flux dif- 
ference split upwind scheme 

• LU-SGS implicit scheme 

© Applications 

• SSME preburner, main combustor and nozzle 

© Status 

• Research code 
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VALIDATION CASES 


© INS3D 

• Internal flow : Channel, Backward- facing step, Rectangular 
duct, Turn-around duct 

• External flow : Circular cylinder (steady-state), Ogive cylinder 

• Juncture flow : Cylinder-plate, Wing-plate, Cavity 

© INS3D-UP , r . 

• Internal flow : Driven cavity, Backward-facing step, Square 

duct with a 90° bend _ 

• External flow : Oscillating plate, Circular cylinder (steady and 

vortex-shedding flows) 

© INS3D-FS 

• Internal flow : Driven cavity 

• External flow : Impulsively started circular cylinder, vortex 
shedding from a circular cylinder 
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SUMMARY OF APPLICATIONS 


0 Space Shuttle Main Engine (SSME) 

(NASA/MSFC, Rocketdyne) 

• Hot Gas Manifold, Main Injector 

• Bearing 

• Impeller /Inducer 

• Preburner 

0 Artificial Heart / Biofluid Mechanics 

(NASA Tech Utilization, Penn State Univ and Stanford Univ) 

© Low Speed Aerodynamics 

• High lift device 

• External flow over Automobiles and Trucks 

© Naval Hydrodynamics (Submarine) 

(DARPA, ONR and David Taylor Research Center) 
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Particle Traces 


for SSME Inducer 





Surface Pressure for SSME Inducer 

(INS3D-LU) 
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FLOW OVER A CIRCULAR CYLINDER AT RE = 105 


Computation 





Experiment 
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ARTIFICIAL HEART 


0 Current artificial devices have problems stemming from fluid dy- 
namic phenomena 

• High shear stress damages the red blood cells and arterial walls 

• Stagnation and secondary flow regions lead to clotting 

• Desire short residence time in artificial environment 

• Large pressure losses cause heart to work harder 

0 Apply CFD technology to analyzing blood flow through artificial 
hearts and to suggest improved design 

• Develop moving boundary capability 

• Apply time accurate flow solvers to Penn State Artificial Heart 

• Develop simple non-Newtonian fluid model 
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Penn State Artificial Heart — INS3D-UP Code 
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CONCLUDING REMARKS 


0 Incompressible and low speed flow simulation codes have been de- 
veloped (INS3D-xx, CENS3D). 

0 Results of computer simulations have made significant impact on 
analysis and redesign of the SSME power head. 

© These codes are being extended to analyze other important real 
world problems. 

0 Future work includes further enhancement of these codes and im- 
provement in physical modeling. 
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SAGE 


A 2-D Self-Adaptive Grid Evolution Code and its Application in 

Computational Fluid Dynamics 

by 

Carol B. Davies 1 , Ethiraj Venkatapathy 2 and George S. Deiwert 3 


Abstract 


SAGE is a user-friendly, highly efficient, 2-dimensional self-adaptive grid code based on Nakahashi 
& Deiwert's variational principles method. Grid points are redistributed into, regions of high flowfield 
gradients while maintaining smoothness and orthogonality of the grid. CPU efficiency is obtained 
by splitting the adaption into 2 directions and applying one-sided torsion control, thus producing a 
1 -D elliptic system that can be solved as a set or tridiagonal equations. 


1 Sterling Software, Palo Alto, Ca. 

2 Eloret Institute, Sunnyvale, Ca. 

3 Chief, Aerothermodynics Branch, NASA Ames Research Center, Moffett Field, Ca. 
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Outline of Presentation 


Brief Description of Self-Adaptive Method (based on 
Nakahashi & Deiwert's variational principles scheme.) 

SAGE code features, including input parameters. 

Application of SAGE to various flow problems: 

• Hypersonic blunt body 

• Supersonic shock impingement 

• Hypersonic inlet with cowl 

• Supersonic plume flow 

• Supersonic inlet 
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Self-Adaptive Grid Method of Nakahashi & 

Deiwert 

Objective: to redistribute points into regions of high flow gradients (utilizing minimization principles) 
while maintaining smoothness and orthogonality between adapted lines. The technique is 
analogous to connecting each node by tension and torsion springs and locating the equilibrium 
position of the resulting system. 

initial grid 



marching direction 
first adaption 
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Adaption Equation 

Splitting the adaption into 2 directions and applying one sided torsion control reduces the problem 
to the following 1-D elliptic system which can be solved as a tridiagonal system of equations: 

(OiASi + i+cOi.iASi-Ci0i=O 

where: 

tension spring constants, coj =l+Afj B 

fj=f(gradient of Qi) {error measure} 

self-adaptive mesh size controlled by A and B: 

A= As max/ As min ~1, B enforces computed ASmin=A.?/?jift 

torsion spring constant, Ci =f(coj,A,*,cell aspect ratio) 

and 6i=f(C7, smoothness, orthogonality) 
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SAGE Input Parameters 


Control parameters : 

X, Ct, Asmin and Asmax 

Adaption parameters : 

Limits of adaption domain: 

Stepping direction. 

Adaption variable (or combination or user defined). 
Boundary spacing controls. 

Addition of grid points. 

Order of interpolation and smoothing. 



SAGE Code Features 


Fast and efficient - -0.005 secs/grid point on a VMS/VAX 
e.g. 30x30 grid takes 5 secs 

Small size: 1 900 lines of code (40% comments) 

Structured code and detailed user guide available - 
easy to customise 

Few user input parameters necessary but many options 
for great flexiblility 

Can be applied to zonal, patched and multiple grids 
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Application 1: Hypersonic Blunt Body 

Flow features: Simple 1 -directional problem, shock shape aligned with grid. 



Initial grid (32x32) 


Initial flow solution 
(density contours) 


Adapted grid 
(i step): Af=.001 
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Application 2: Cowl Lip with Shock Impingement 

Flow features: blunt body shocks, impinging shock and shear layers 




Initial Grid(91x46) 


Initial density solution 


Adapted grid 

h As m j n =.25, Cp.7, A^.005 
j: As m j n =.25, orthe=false 
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Application 3: Hypersonic Inlet with Cowl 

Flow features: Cowl blunt body shock, Mach stem & reflected shocks 


4 ^ 

'-J 
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Application 3 continued: Adapted Grid 

Example of zonal adaption 


i-direction - full domain adaption 
j-direction - 3 zonal adaptions: 


1 ) cowl region 

2) lower inlet region 

3) upper inlet region 




Application 4: Plume Flow 

Flow features: Outer shear layer, barrel shock, Mach 
disc, reflected shock, triple-point shear layer 



Initial grid and 
Mach contours 




Adapted grid and 
solution (after 3 
iterations) 



Comparison with 
shadowgraph 
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Application 5: Supersonic Inlet 

Flow features: Corner shock, reflected shock and 

expansion fan 

Initial grid (101x81) and density solution contours 
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Application 5 continued: Adapted 
grid (marching in j) and Solution 

Input parameters: As m j n =.25, As max =2.5, X=.0005 




251 


Application 5 continued: Adapted 
Grid (marching in i) and Solution 

Input parameters: jstep=false, As m j n =.25, X^.001 
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Concluding Remarks 


SAGE is a new 2-D self-adaptive grid code that is 
user-friendly, flexible and efficient. 

Appropriate for a variety of CFD applications. 

Use of the SAGE code will efficiently improve the 
flowfield solution 
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Time Dependent Viscous Incompressible Navier-Stokes Equations 


By 

John W. Goodrich 

Computational Fluid Dynamics Branch 
Internal Fluid Mechanics Division 
NASA Lewis Research Center 
Cleveland, Ohio 44135 
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TIME DEPENDENT INCOMPRESSIBLE NAVIER-STOKES EQUATIONS 


au _ , . 1 

— + V • (uu) — — A u = — Vp + F, for x in 0, and t > 0, 

V • u = 0, for x in H, and t > 0. 

* fl is a bounded open region in R 2 , with boundary 3D. 

* Initial and boundary conditions must be supplied. 

* F is the volume force per unit mass, assumed to be 0. 

=* > The continuity equation is not given in a time evolution form. 

The pressure gradient couples the continuity equation to the momentum equations. 
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STREAMFUNCTION EQUATIONS FOR UNSTEADY INCOMPRESSIBLE FLOW 


with 


d A tp 1 2 dtp dtp dtp dtp 

~dT ~ Te A * + a? ^ ■ for * in n ' and ‘ > °; 

/ .x , , x dtp 

«(*.*) = 0^, and v(x,t) = - — , for x in fl, and t > 0. 


* For Cl in R 2 . 

* Initial and boundary conditions must be supplied. 

=* Vorticity and pressure do not enter into the streamfunction formulation. 
^ The velocity solution is always divergence free, and so incompressible. 
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THE STREAMFUNCTION ALGORITHM FOR UNSTEADY INCOMPRESSIBLE FLOW 


La(z" + 1 ) - -^-Bi(z" + 1 ) 
v ' 2Re v 1 

= La( 2 n ) + ^ Bi(n _3At 


<5* ^,(z n )La(z")^ - 6 y (s x { z”)La(z n ) 


+ Y [«. («. (2- “ 1 M*“- 1 )) - «„ («, (*' ' - 1 )La(2" - : ‘ )) 


with 


u n . = 


2Ay 


«y + 1 “ *"y-i). and <y = ~ 2Ax ( 2 *" + 1,J “ 


* La and Bi are central difference approximations to the Laplace and Biharmonic operators. 

* 6 X and 6 y are conventional centered difference operators. 


In R 2 there is one unknown { z ? y } per grid cell instead of three. 

The velocity components and streamfunction are all defined at each grid point. 
The discrete solution is exactly incompressible, 6 X (u?V) + 6 y (uTV) = 0. 
Stability limit is Courant number <1. 
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A MULTIGRID SOLVER FOR THE LINEAR IMPLICIT EQUATIONS 
La(z" +1 ) - ^Bi(z ri + 1 ) = Source Term{z n ,z n ~ 1 ) 


=>• Use a multigrid solver for the implicit equations at each time step. 

* The Biharmonic operator is factored as two Laplacians. 

On a 256 by 256 fine grid, 7 grid levels are used in 6.8 MBytes storage. 

* Point Gauss-Seidel smoothing, linear restriction and prolongation. 

* A V-cycle with 3 iterations per grid level while coarsening, none while refining. 

* 10 to 15 iteration cycles reduce residuals to less than 5.0 x 10 -11 . 
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STREAM FUNCTION CONTOURS STREAM FUNCTION CONTOURS - NORMALIZED VECTOR 'PLOTS 

R e = 5k, 128*128 grid, 1=491.80625 Re=5k. 128*128 grid, t=491. 80625 
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STREAM FUNCTION SURFACE 
Re = 5k, 96*192 grid, t=4000 
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KINETIC ENERGY SURFACE 
Re=5k, 96*192 grid, t=4000 
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STREAMFUNCTION CHANGE PER TIME STEP 
Relative LI norm for the change 
Re=5k, 96*192 grid, 3800<=t<=4 100 
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STREAM FUNCTION CONTOURS - NORMALIZED VECTOR PLOTS 
Re = 5k, 96*192 grid, t = 4100.25 
0.0<=x< = 1.0 and 0.85<=y< = 1 .35 
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STREAM FUNCTION CONTOURS - NORMALIZED VECTOR PLOTS 
Re = 5k, 96*192 grid, t=4101.25 
0.0<=x< = 1 .0 and 0.85<=y< = 1 .35 
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STREAM FUNCTION CONTOURS - NORMALIZED VECTOR PLOTS 
Re = 5k, 96*192 grid, t=4101.50 
0.0<=x< = 1.0 and 0.85<=y< = 1 .35 
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STREAM FUNCTION CONTOURS - NORMALIZED VECTOR PLOTS 
Re = 5k, 96*192 grid, t=4102.50 
0.0<=x< = 1 .0 and 0.85<=y< = 1 .35 


Os 

SO 


U22I 


1 .099 


0 .977 


0.854 

0 






s*' 


}f^\\\\^\\NNN\\\\\\\\V 
^\\\\\\\\\\\\\\\\\\\\\\\*: 
VJA 

w \ WWW *W 
W\*>* ^ 

Mi i wwsss vw w w w 

M v *■ v ^ **wv ^ wC^ 

*/ yl >> N s\\sss'^ * “ — 

i << i/ 4-^^^^^^^ r: 

l ^,<*T^.ww\w 
M t W WNW*>^V*. 

l <\eV\\.\ \ \ \\\W^ 


rrm 

i-4t'*'4r'*' tty. 


ft'***'* k' / * *. 


risx'x'* tftc 


//////l////J//a/ 
////////////////.U/^ 
/s/z/s/j ////// //// * 
•szzzzzz/zzz//7Zz///«' 
ZZ/ZZZ zZ //////// z// /v 
z z z z z z z Z& z z z / /% z / / /V*. 

ZrfZZZZz/ ^ ? ' ' 


^w wkk 


*~4r~4c'4r'4r' *<*'*-*'#' 4^. 


■M* 




„ 'z Z z z z z/z z z/z z zz 
r zzzz zz/z Z Z/ ZZ z^{ }'f-i * f i , 
z z z z z zVz z z’/z*/ z*\ t tt'tff f f } 
■*s*s*ss*yi//u/74</ { ttttfttf ft 
. 'z^^u'zzz^zzzyzZzz ft ft fit f til 

ls,/„f r t' n t [ 

vzzzzSz/z z )»,*/• fisx? / ?? r ttl 
**/.*. ft**???? ft 

* *^z***s s’ s s s s s 

\j$ 

■vv' $Sl littltt f '^' , -*-C'> v * 

^ 5J N«ijl|l/lte||! 

555555535^>* UNI j 


r**zzje Z * 
ZrzzZz* 


\\\\\\\\<^ t I | 

► i * 


U/*T 


liZZ** 
azzzz^^ 
i azzz^^^ 


11/ Z ZZZk'*'*'*'* 

a/ zz * ** **'«'}**'+■ 


y iii i /ZZZ***'*’«'«'*'*~* rA 

u am/ tst/y*'*'*'*'*'*'^ 

lam ///</!/**’*'*'*'*'«'*'*'*■ 
la am / zzz zz**'**'*'*'*^** 

\ aim /• ///**!/ ****'«'*'*'*'*'* 
WiiUM / ///*!/ */*/«*'*'*'**'* 

1 1 1 1 a /z////zzzz *'*’*'*<'*'*’ 

■ f i/ / / / 


oo 


hB#* 


JTaVV * V zzzzzuzvw y * 

0.25 


k wwwww\\wwo ^ ; 

. ^W\WWW\W\WW£/ l 

vwwww\w w \w wwSc / 1 

i 

.\\ww\\\\\\\*Xl / i „ 

* W W W W W W S A \ \W ^ + 

‘ \ss\ \ ^ \ \ ^ * \ \ \ *. \ V i u 

— ™ , > ; 

Kx \ ^ \ \ ^ \ \ \ \ \ \ \ ^a \ \ \ \ \x ; 

w<V, \ \ N \ N S\\ \ *\ 's \ \ \ \ l 

\\ \ \ \ \ \ \ \ \\ \ \ \ \ \* 

,. i 5iiSJi55555555y4S^5 3 } 

v \ \ \ 




\ \ ^ h v \ \ \ \ \ \ 1 VA \ \ U } 


0.50 


0.75 


1 .00 


X 


270 


SUMMARY: A NEW ALGORITHM 

=> has one unknown per grid cell in two space dimensions; 

=>• requires storage that increases linearly with the number of grid points; 

=>• CPU time per time step increases linearly with the number of grid points; 
==>■ is second order accurate in both time and space; 

==>• stability limit is Courant number < 1; 

=> is robust with respect to Reynolds number. 


SUMMARY: A NEW PERIODIC FLOW SOLUTION 
=r- is exactly periodic; 

==>■ does not use a time dependent forcing term; 

==> has no periodic or artifical throughflow boundary conditions; 
=> is probably driven by the wall jet descending from the lid; 

==> is evidence of a Hopf bifurcation; 

=$■ may lead to period doubling bifurcations and a chaotic flow. 


..4 
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CFD for Applications to Aircraft Aeroelasticity 

Guru P. Guruswamy 
Applied Computational Fluids Branch 
NASA Ames Research Center 
Moffett Field, California 

Abstract 

Strong interactions of structures and fluids are common in many engineering 
environments. Such interactions can give rise to physically important phenomena 
such as those occurring for aircraft due to aeroelasticity. Aeroelasticity can signif- 
icantly influence the safe performance of aircraft. At present exact methods are 
available for making aeroelastic computations when flows are in either the linear 
subsonic or supersonic range. However, for complex flows containing shock waves, 
vortices and flow separations, computational methods are still under development. 
Several phenomena that can be dangerous and limit the performance of an aircraft 
occur due to the interaction of these complex flows with flexible aircraft components 
such as wings. For example, aircraft with highly swept wings experience vortex in- 
duced aeroelastic oscillations. Correct understanding of these complex aeroelastic 
phenomena requires direct coupling of fluids and structural equations. This paper 
provides a summary of the development of such coupled methods and its applica- 
tions to aeroelasticity since about 1978 to present. A part of the paper discusses 
the successful use of the transonic small perturbation theory (TSP) coupled with 
structures. This served as a major stepping stone for the current stage of aeroe- 
lasticity using CFD. The need for the use of more exact Euler /Navier-Stokes(ENS) 
equations for aeroelastic problems is explained. The current development of un- 
steady aerodynamic and aeroelastic procedures based on the ENS equations are 
discussed. The paper illustrates aeroelastic results computed using both TSP and 
ENS equations. 
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HISTORY OF CFD APPLICATIONS TO AEROELASTICITY 

BASED ON UNSTEADY TIME ACCURATE METHODS 



NAVIER 

STOKES 

1988 

? 

? 
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? 
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MAJOR ISSUES FOR ADVANCED CFD METHODS 


Computational speed 

- Aeroelastic computations require two orders more compu- 
tational time than steady computations 

Time accuracy 

- An essential requirement for accurate aeroelastic computa- 
tions 

Grid and its unsteady movement v 

- Time accuracy between zones 

Validity of turbulence models for unsteady and separated flows 

Robustness of solution methods 

- Other issues like artificial viscosity, upwinding, etc. 
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APPROACH FOR COMPUTER SIMULATION 


• GOVERNING EQUATIONS 

- Aerodynamics : 3-D Euler/Navier- Stokes equations (ENS) 
and transonic small perturbation equation (TSP) 

- Aeroelastic : Modal equations of motion 

• ALGORITHM 

- Aerodynamics : time accurate finite difference methods 
based on alternate direction implicit schemes 

- Aeroelastic : Simultaneous time integration method 


• Note 

- ENS computations are made using aeroelastic adaptive dy- 
namic grids 

- For TSP computations, aerodynamic and structural prop- 
erties of the fuselage tip stores, and control surfaces are 
modeled 



COUPLED AEROELASTIC EQUATIONS OF MOTION 

• Deformed shape is a sum of modal coordinates 

• Equations are solved by simultaneous integration technique 

• Equations of Motion 

- Assuming displacement vector {d} = [</>]{#} where [</>] is 
the modal matrix and {qr} is the generalized displacement 
vector, the aeroelastic equation of motion is 

[M]{q} + [<?]{«} + {K]{q} = {F} 

[M], [G], and [K] mass, damping and stiffness matrices 

{F} = (|)pU^[</>] T [A]{AC p } is the aerodynamic force vector 

[A] is the diagonal area matrix of the aerodynamic control points. 



OF POOR QUALITY 


MODES OF A RECTANGULAR WING 




SOME APPLICATIONS OF TSP THEORY 


Transonic flutter boundaries of transport and fighter wings 

Aeroelasticity of a variable sweep wing (B-l wing) 

Aeroelasticity of wings with tip stores 

Aeroelasticity of wings with active control surfaces 

Aeroelasticity of full span wing-body configurations (Symmetric 
and Asymmetric modes) 
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TYPICAL RESULTS FROM TSP CODE ATRAN3S 
Full-Span Wing-Body Aeroelasticity 


ANTISYMMETRIC MODES 
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DEVELOPMENT OF ENSAERO 


• PURPOSE 

- To develop an aeroelastic code to solve Euler/Navier Stokes 
equations coupled with structural equations of motion for 
full aircraft 

• CHARACTERISTICS 

- Solves either Euler or Navier Stokes equations 

- Models structure by either modal or finite element equa- 
tions 

- Includes aeroelastic configuration adaptive grid scheme 

- Modular to adopt different finite difference schemes 

- Transportable to different computer configurations 

• ENSAERO- version 2.0 

- Solves Euler/Navier Stokes equations with modal structural 
equations of motion for wings 
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CONFIGURATION ADAPTIVE DYNAMIC GRID 

• Grids are generated by an algebraic method 

• Grids conform to the wing surface defined by displacements {d} 

• Grids are generated every time-step of integration 

• Time metrics are computed every time step 

= UrCy 

'H t = 3 'tVx Ur'Hy %tVz 

Ct ~ Cx UrCy %tCz 

J -1 = + x^Zr, + X v y i z i - x^y^Zr, - x n y ^ - x ( y v zt 

• Note - Present technique can be used for both structured and 
unstructured grids 



VORTEX DOMINATED UNSTEADY PRESSURES 

Navier-Stokes Computations 

Rectangular Wing in ramp motion , AR = 4.0, NACA0015 




282 ' 


COMPARISON OF UNSTEADY PRESSURES 
BETWEEN RIGID AND FLEXIBLE WINGS 
(Navier-Stokes Computations) 
Rectangular Wing in ramp motion , AR = 4.0, NACA0015 


GRID 151x20x40 
Moo = 0.50, A = 0.30, 
R e = 60000.0 
a = 20°, r = 0.94 

Rigid Wing 

Flexible Wing 




X/C 
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CONCLUDING REMARKS 

• During last decade TSP applications have progressed from air- 
foils to almost full aircraft 

- Computational speed has increased by a factor of about 100 

- Robust codes such as ATRAN3S are now available 

- Applied for advanced applications such as active controls 

• Euler/Navier Stokes (ENS) equations are currently being used 
for aeroelastic problems of wings 
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FUTURE DIRECTIONS 


• Improve time accurate Euler /Navier Stokes(ENS) algorithms 

• Extend unsteady ENS algorithms for full aircraft configurations 

• Couple advanced CFD methods with advanced CSM methods 

i 

• Conduct research in unsteady aerodynamics and aeroelasticity 
of full aircraft at high angles of attack 

• Maintain TSP codes for immediate industrial use 
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FUTURE DIRECTIONS (continued) 
Algorithm Development 


• Typical results from a new upwind scheme that will be imple- 
mented in ENSAERO 


Vortical flow on a 75° delta wing 
Steady pressures at M = 2.80, a = 16.0° 



0.0 0.1 0.2 0.3 
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FUTURE DIRECTIONS (continued) 


• Typical steady results from Transonic Navier Stokes (TNS) code 

• Unsteady algorithm will be implemented in full aircraft TNS 
code 
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APPLICATION OF UNSTRUCTURED GRID METHODS TO 
STEADY AND UNSTEADY AERODYNAMIC PROBLEMS 


John T. Batina 

Unsteady Aerodynamics Branch 
NASA Langley Research Center 
Hampton, Virginia 23665-5225 


Abstract 


The presentation summarizes recent work in the Unsteady Aerodynamics Branch at NASA 
Langley Research Center on developing unstructured grid methods for application to steady and 
unsteady aerodynamic problems. The CAP-TSD transonic aeroelasticity code, which is based on 
the transonic small-disturbance (TSD) theory, is described first to provide background 
information to put the present work in context. The CAP-TSD code is the most fully-developed 
code for aeroelastic analysis of complete aircraft configurations at the TSD equation level and 
has been widely accepted throughout the U.S. aerospace industry. Currently, aeroelastic 
analysis capabilities are being developed at NASA Langley for the Euler and Navier-Stokes 
equations based on both structured and unstructured grids. The purpose of the presentation is to 
describe the development of unstructured grid methods which have several advantages when 
compared to methods which make use of structured grids. Unstructured grids, for example, 
easily allow the treatment of complex geometries, allow for general mesh movement for 
realistic motions and structural deformations of complete aircraft configurations which is 
important for aeroelastic analysis, and enable adaptive mesh refinement to more accurately 
resolve the physics of the flow. The presentation is therefore organized in three parts 
including: (1) steady Euler calculations for a supersonic fighter configuration to demonstrate 
the complex geometry capability; (2) unsteady Euler calculations for the supersonic fighter 
undergoing harmonic oscillations in a complete-vehicle bending mode to demonstrate the general 
mesh movement capability; and (3) vortex-dominated conical-flow calculations for highly- 
swept delta wings to demonstrate the adaptive mesh refinement capability. The basic solution 
algorithm is a multi-stage Runge-Kutta time-stepping scheme with a finite-volume spatial 
discretization based on an unstructured grid of triangles in 2D or tetrahedra in 3D. The moving 
mesh capability is a general procedure which models each edge of each triangle (2D) or 
tetrahedra (3D) with a spring. The resulting static equilibrium equations which result from a 
summation of forces are then used to move the mesh to allow it to continuously conform to the 
instantaneous position or shape of the aircraft. The adaptive mesh refinement procedure 
enriches the unstructured mesh locally to more accurately resolve the vortical flow features. 
These capabilities are described in detail along with representative results which demonstrate 
several advantages of unstructured grid methods. The presentation further discusses the 
applicability of the unstructured grid methodology to steady and unsteady aerodynamic problems 
and suggests directions for future work. 
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PRESENTATION OVERVIEW 


Background on CAP-TSD transonic aeroelasticity code 
Unstructured grid methods 

- steady and unsteady Euler calculations for a supersonic 
fighter configuration 

- vortex-dominated conical-flow calculations including 
adaptive mesh refinement 


Concluding remarks 



CAP-TSD: COMPUTATIONAL AEROELASTICITY PROGRAM - 

TRANSONIC SMALL DISTURBANCE 


• Based on time-accurate approximate factorization algorithm 

• Complete aircraft modeling involving arbitrary combinations 
of lifting surfaces and bodies 

• Static and dynamic aeroelastic analysis 

• Aircraft trim capability 

• Longitudinal short-period response 

• Entropy and vorticity effects included to treat cases with 
strong shock waves 
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CAP-TSD INSTANTANEOUS PRESSURES ON F-16C AIRCRAFT 

DUE TO RIGID PITCHING MOTION 


• M „ = 0.9, <x 0 = 2.38°, k = 0.1 
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ADVANTAGES OF UNSTRUCTURED GRID METHODOLOGY 


• Allows the treatment of complex geometries 

• Allows general mesh movement for realistic motions and 
structural deformations of complete aircraft 

• Enables adaptive mesh refinement to more accurately 
resolve the physics of the flow 
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DESCRIPTION OF EULER SOLUTION ALGORITHMS 


• Four-stage Runge-Kutta time-stepping scheme 

• Finite-volume spatial discretization on unstructured grids 
of triangles in 2D or tetrahedra in 3D 

• Adaptive blend of harmonic and biharmonic operators for 
artificial dissipation 

• Enthalpy damping, local time-stepping, and implicit residual 
smoothing to accelerate convergence to steady state 

• Dynamic mesh algorithm employed for unsteady applications 
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SURFACE TRIANGULATION FOR LANGLEY FIGHTER 


• Total grid has 13,832 nodes and 70,125 tetrahedra 
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STEADY PRESSURE CONTOURS ON LANGLEY FIGHTER 


• Moo = 2.0 and <x 0 = 0° 
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PRECEDING PAGE BLANK NOT FILMED 


OVERVIEW OF DYNAMIC MESH ALGORITHM 

• Each edge of each triangle is modeled using a spring 

• Spring stiffness is inversely proportional to the length of 
the edge 



airfoil 


• Points on outer boundary of grid are fixed 

• Locations of points on inner boundary of grid are specified 



300 


OVERVIEW OF DYNAMIC MESH ALGORITHM 


• Displacement of interior nodes determined by solving the 
static equilibrium equations 

• Equations solved using predictor-corrector procedure 


- new locations of nodes predicted by extrapolation 




5 


= 2 




2 8 


n 





- locations corrected by several Jacobi iterations 




ASSUMED BENDING MODE FOR LANGLEY FIGHTER 
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INSTANTANEOUS PRESSURE CONTOURS ON LANGLEY FIGHTER 
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DESCRIPTION OF CONICAL EULER/NAVIER-STOKES ALGORITHM 


• Solves the conical Euler/Navier-Stokes equations 

• Multi-stage Runge-Kutta time-stepping with finite-volume 
spatial discretization on unstructured grid of triangles 

oj 

o 

• Scheme is a zonal method 

• Mesh enrichment capability enables automatic refinement in 
regions of high flow gradients 

• Presently limited to laminar flow 
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EFFECTS OF MESH ENRICHMENT ON CONICAL EULER 
VORTICAL FLOW SOLUTION 

• 75° swept flat plate delta wing at M~ = 1.7, a = 12°, p = 8° 

• Meshes • Total pressure loss contours 





2555 nodes 



3735 nodes 



CROSS FLOW VELOCITY VECTORS FROM CONICAL 
NAVIER-STOKES SOLUTION 

• 70° swept elliptic cone delta wing; thickness ratio 14:1 

• = 2.0, a = 10°, and Re = 5 x 10 5 
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CONCLUDING REMARKS 


• CAP-TSD code developed for transonic aeroelastic analysis 
at the transonic small-disturbance equation level 

• Unstructured grid methods for the Euler and Navier-Stokes 
equations under development for steady and unsteady 
aerodynamic applications 

- Steady Euler calculations for a supersonic fighter 
demonstrated ability to treat complex geometries 

- Unsteady Euler calculations for a supersonic fighter 
demonstrated general dynamic mesh capability 

- Vortex-dominated conical-flow calculations 
demonstrated adaptive mesh refinement capability 
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FIGHTER AIRCRAFT 


Chairman: 

Terry L. Holst 

Chief, Applied Computational Fluids Branch 
NASA Ames Research Center 
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Grid Generation and Inviscid Flow Computation 
About Aircraft Geometries 

Robert E. Smith 

Analysis and Computation Division 
NASA Langley Research Center 

Abstract 


Grid generation and Euler flow about fighter aircraft are described. A 
fighter aircraft geometry is specified by an area ruled fuselage with an inter- 
nal duct, cranked delta wing of strake/wing combinations, canard and/or 
horizontal tail surfaces, and vertical tail surfaces. The initial step before 
grid generation and flow computation is the determination of a suitable 
grid topology. The external grid topology that has been applied is called a 
dual-block topology which is a patched C 1 continuous multiple-block system 
where, inner blocks cover the highly-swept part of a cranked wing or strake, 
rearward inner-part of the wing, and tail components. Outer-blocks cover 
the remainder of the fuselage, outer-part of the wing, canards and extended 
to the fax field boundaries. The grid generation is based on transfinite in- 
terpolation with Lagrangian blending functions. This procedure has been 
applied to the Langley experimental fighter configuration and a modified 
F-18 configuration. Supersonic flow between Mach 1.3 and 2.5 and angles 
of attack between 0° and 10° have been computed with associated Euler 
solvers based on the finite- volume approach. When coupling geometric de- 
tails such as boundary layer diverter regions, duct regions with inlets and 
outlets, or slots with the general external grid, imposing C 1 continuity can 
be extremely tedious. The approach taken here is to patch blocks together 
at common interfaces where there is no grid continuity, but enforce conser- 
vation in the finite-volume solution. The key to this technique is how to 
obtain the information required for a conservative interface. We have used 
the Ramshaw technique which automates the computation of propor- 
tional areas of two overlapping grids on a planar surface and is suitable for 
coding. We have generated internal duct grids for the Langley experimen- 
tal fighter configuration independent of the external grid topology, with a 
conservative interface at the inlet and outlet. 
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FEATURES 

• MULTIPLE-BLOCK STRUCTURED GRIDS 

• FINITE-VOLUME EULER SOLVERS 

• CONSERVATIVE INTERFACE BETWEEN GRID BLOCKS 




CONTENTS 


• FIGHTER CONFIGURATIONS 

• BOUNDARY GRIDS >1, 

• VOLUME GRID GENERATION 

• CONSERVATIVE INTERFACING 

• EULER FLOW SOFTWARE AND SOLUTIONS 

• VIDEO DISPLAY 

• CONCLUSIONS AND COMMENTS 


FIGHTER CONFIGURATIONS 



EXPERIMENTAL FIGHTER 


MODIFIED F-18 
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BOUNDARY GRIDS 


• CONFIGURATION SURFACE 

1. Component Cross Sections 

2. Interpolation, Fitting and Smoothing 

3. Grid Point Distributions 

• INTERMEDIATE BOUNDARY SURFACES 

1. Algebraic Functions 

2. Embedded Grid Point Distributions 


• FAR FIELD BOUNDARY SURFACES 


BOUNDARY-SURFACE GRIDS 



Original Definition 


Grid-Surface Definition 



317 


VOLUME GRID GENERATION 


• GRID TOPOLOGY 

* Block Location and Interfaces 
•TRANSFINITE INTERPOLATION 

^Various Interpolation Techniques 

• GRID SPACING CONTROL 

* Exponential Functions 
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CONSERVATIVE GRID INTERFACES 


Rai’s Approach 


• Ramshaw Redistribution 
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INVISCID-COMPRESSIBLE FLOW 


•FINITE-VOLUME METHODOLOGY 

•CUSTOM EULER SOLVER 

‘CONTINUITY, MOMENTUM AND ENERGY 
*3 rd ORDER RUNGE-KUTTA 

•GENERAL MULTI-BLOCK EULER SOLVER 

‘CONTINUITY, MOMENTUM AND CONST. ENTHALPY 
*4 fft ORDER RUNGE-KUTTA 



SUMMARY OF SOLUTIONS 

Mach Number = 2 
Coefficient of pressure 
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| EXPERIMENTAL FIGHTER CONFIGURATION 
CONTOURS OF PRESSURE COEFFICIENT 

Mach number = 2 Angle of attack = 0° 




PRECEDING PAGE BLANK NOT FILMED 325 



326 


CONCLUSIONS 


ALGEBRAIC GRID GENERATION FEASIBLE 


• CONSERVATIVE BLOCK INTERFACES POSSIBLE 


FINITE- VOLUME EULER SOLVERS USED 


• FLOW CHARACTERISTICS IN GOOD AGREEMENT 


VORTEX COMPUTATION NEEDS MORE STUDY 
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A ZONAL NAVIER-STOKES METHODOLOGY FOR FLOW 
SIMULATION ABOUT A COMPLETE AIRCRAFT 

Jolen Flores 


Abstract: The thin-layer, Reynolds-averaged, Navier-Stokes equations are 
used to simulate the transonic viscous flow about the complete F-16A fighter air- 
craft. These computations demonstrate how computational fluid dynamics (CFD) 
can be used to simulate turbulent viscous flow about realistic aircraft geometries. 
A zonal grid approach is used to provide adequate viscous grid clustering on all 
aircraft surfaces. Zonal grids extend inside the F-16A inlet and up to the com- 
pressor face while power on conditions are modeled by employing a zonal grid 
extending from the exhaust nozzle to the far field. Computations are compared 
with existing experimental data and are in fair agreement. Computations for the 
F-16A in side slip are also presented. 


Applied Computational Fluids Branch 
NASA Ames Research Center 
Moffett Field, CA. 94035 

Presentation at NASA OAST CFD CONFERENCE 
March 7-9, 1989 
NASA Ames Research Center 
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CANOPY 


OBJECTIVE 

FUSELAGE 

> To numerically simulate viscous 
transonic flow about realistic 
aircraft configurations. 


STRAKE 



INLET 


EXHAUST PLUME 


• MOTIVATION 

> Technical demonstration of state-of-the-art CFD research 

> To provide bench-mark calculations (validated by measurements) 

> Reveal areas requiring future research emphasis 

> Catalyst for future cooperative efforts between NAS, 

aerospace industry and academia 
industrial use for prediction of integrated aircraft performance 

• APPROACH 

> Thin-layer Navier- Stokes equations 

> Zonal grid approach 
t> Modular program 
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F-16 





OF POOR QUALITY 


INLET /DIVERTER PRESSURE COEFFICIENT COMPARISONS 

Moo = 0.9, a = 6.0°, Re c = 4.5X10 6 




PRESSURE CONTOURS ON UPPER SURFACE OF F- 16 A 
(Mach=0.9, Alpha=6.0°, Beta=0.0°, Re c =4.5x10 6 ) 
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CROSS-FLOW INFLUENCE ON UNRESTRICTED PARTICLE TRACES 
(Mach=0.9, Alpho=6.0°, Beta=5.0 o , Re c =4.5x10 6 ) 
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SUMMARY 


• Benchmark Navier-Stokes simulation of a complete aircraft including sideslip 

g = 0.0° 

>Good comparison with Cp, Cp, and C d 
>Successful implementation of internal inlet grids 

U) 

>Successful simulation of power-on conditions 
oConvergence in 5000 iterations/ 25 hours of cpu 

P = 5.0° 


c>Pressure contours/particle traces indicate proper physical trends 
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NUMERICAL SIMULATION OF F 18 FUSELAGE FOREBODY 
FLOWS AT HIGH ANGLES OF ATTACK 


Lewis B. Schiff* 

Russell M. Cu mmi ngs t 
Reese L. Sorenson* 
Yehia M. Rizk^ 

NASA Ames Research Center 
Moffett Field, CA 94035 


Abstract 

As part of the NASA High Alpha Technology Program, fine-grid Navier-Stokes solu- 
tions have been obtained for flow over the fuselage forebody and wing leading edge exten- 
sion of the F/A-18 High Alpha Research Vehicle at large incidence. The resulting flows 
are complex, and exhibit crossflow separation from the sides of the forebody and from the 
leading edge extension. A well-defined vortex pattern is observed in the leeward-side flow. 
Results obtained for laminar flow show good agreement with flow visualizations obtained in 
ground-based experiments. Further, turbulent flows computed at high-Reynolds-number 
flight-test conditions (M^ = 0.2, a = 30°, and Rec = 11.52 x 10 6 ) show good agreement 
with surface and off-surface visualizations obtained in flight. 


* Research Scientist, Applied Computational Fluids Branch. 

f National Research Council Research Associate. Associate Professor, on leave from 
California Polytechnic State University, Aeronautical Engineering Department. 

J Member of the Professional Staff, Sterling Federal Systems, Inc., Palo Alto, CA. 
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OBJECTIVE 


• DEVELOP FLIGHT-VALIDATED DESIGN METHODS THAT 
ACCURATELY PREDICT THE AERODYNAMICS OF AIRCRAFT 
MANEUVERING AT LARGE ANGLES OF ATTACK 

APPROACH 

• UTILIZE A THREE-DIMENSIONAL NAVIER-STOKES CODE, WITH 
SUITABLE GRIDS AND AN EDDY-VISCOSITY TURBULENCE 
MODEL, TO COMPUTE HIGH-ALPHA FLOWS OVER THE F-18 
FUSELAGE FOREBODY AND LEX 

• VALIDATE THE NUMERICAL RESULTS BY COMPARISON WITH 
FLIGHT-TEST DATA OBTAINED ON THE NASA F-18 HIGH 
ALPHA RESEARCH VEHICLE (HARV) 
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GOVERNING EQUATIONS 

dQ dF dG dH _ J _5S 

57 + ai - + + ac _ fle 9C 

• THIN-LAYER NAVIER-STOKES EQUATIONS 

• CURVILINEAR, BODY-CONFORMING COORDINATES 

• HIGH REYNOLDS NUMBER FLOWS 

• LAMINAR VISCOSITY FROM SUTHERLAND’S LAW 

• ALGEBRAIC EDDY-VISCOSITY MODEL CORRECTED FOR 
CROSSFLOW SEPARATION 
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NUMERICAL METHOD 


{ 


I+h 


^(A+) + S Q C - ±- e 6 c M 






R.H.S. 


• TWO-FACTORED ALGORITHM (F 3 D) 

• FIRST OR SECOND-ORDER ACCURACY IN TIME 

• SECOND-ORDER SPATIAL ACCURACY 

• FLUX-VECTOR SPLITTING AND UPWIND DIFFERENCING IN 
£ (STREAMWISE) DIRECTION 

• CENTRAL DIFFERENCING IN THE 77 (CIRCUMFERENTIAL) 
AND ( (RADIAL) DIRECTIONS 

• COMBINATION OF SECOND AND FOURTH-ORDER SMOOTHING 
USED IN THE 7 7 AND C DIRECTIONS 

• SMOOTHING TERMS SCALED BY q/q,*, 

• SINGLE-BLOCK AND TWO-BLOCK GRIDS USED 



349 


F-18 FOREBODY TWO-BLOCK GRID 
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ONE-BLOCK GRID: F-18 FOREBODY CLOSE-UP 




TWO-BLOCK GRID 


5m:' 


F-18 FOREBODY CLOSE UP 
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TANGENT OGIVE-CYLINDER SINGLE-BLOCK GRID 



NASA-AMES 
HIGH ALPHA GROUP 



COMPUTED SURFACE OIL FLOW 

= 0 . 2 , a = 20 ° 

Re n = 5.0 x 10 6 (TURBULENT) 
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HELICITY 

Moo = 0.2, a = 20° 

Re D = 5.0 x 10® (TURBULENT) 


SURFACE FLOW PATTERN 

Meo = 0.2, * = 30’ 

Re c = 11,540,000 (TURBULENT) 
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FLIGHT SURFACE FLOW VISUALIZATION 

QUARTER VIEW, a = 30° 



SECONDARY VORTEX 
SEPARATION LINE 


PRIMARY VORTEX 
SEPARATION LINE 


SURFACE 

STREAMLINE 


dL 
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HELICITY DENSITY 

Moo = 0.2, * = 30° 

Re c = 11,540,000 (TURBULENT) 
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Wingtip Photograph of F-18 

a - 20.8° and /? = +1.15° 



LEX vortices visualized using smoke 
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SUMMARY REMARKS 


• NAVIER-STOKES COMPUTATIONS FOR HIGH-ALPHA SEPARATED 
TURBULENT FLOW ABOUT THE F-18 (HARV) FUSELAGE FOREBODY 
AND LEX SHOW GOOD AGREEMENT WITH FLIGHT-TEST DATA 

• ONLY MINOR DIFFERENCES BETWEEN SINGLE-BLOCK AND 
TWO-BLOCK RESULTS 

• EFFECTS OF INCREASING INCIDENCE CONSISTENT WITH 
EXPERIMENT 

• CFD RESULTS HAVE GIVEN NEW INSIGHT INTO HIGH-ALPHA 
FLOW STRUCTURE 

• COMPUTATION-TO-FLIGHT PREDICTIONS OF FULL F-18 
CONFIGURATIONS ARE NEXT STEP 

• USE OF CFD AS A DESIGN TOOL FOR VORTEX CONTROL 
CONCEPTS IS AT HAND 
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NAVIER-STOKES SOLUTIONS ABOUT THE F/A-18 
FOREBODY-LEX CONFIGURATION 


Farhad Ghaffari 
Vigyan Research Associates 

James M. Luckring, James L. Thomas 
NASA Langley Research Center 

Brent L. Bates 

g Vigyan Research Associates 


Abstract 


Three-dimensional viscous flow computations are presented for the F/A-18 forebody-LEX geome- 
try. Solutions are obtained from an algorithm for the compressible Navier-Stokes equations which in- 
corporates an upwind-biased, flux-difference-splitting approach along with longitudinally- patched grids. 
Results are presented for both laminar and fully turbulent flow assumptions and include correlations 
with wind tunnel as well as flight-test results. A good quantitative agreement for the forebody surface 
pressure distribution is achieved between the turbulent computations and wind tunnel measurements at 
Moo = 0.6. The computed turbulent surface flow patterns on the forebody qualitatively agree well with 
in-flight surface flow patterns obtained on an F/A-18 aircraft at = 0.34. 
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Overview 


• Navier-Stokes Formulation 

o CFL-3D 

• Grid Generation 

o Transfinite interpolation 

• Results 

o Laminar, turbulent flow 
o Comparisons with wind-tunnel experiment 
o Comparisons with flight test 


Summary 
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Grid Generation - Transfinite interpolation 

• H-0 topology 

• Far field 

o Inflow, outflow ^ 1 c 
o Radial w 1.5 c 

• Baseline grid 

o Block 1 : 31 x 65 x 27 
o Block 2 : 65 x 65 x 31 
o Approximately 185,000 points 
o 2 for wind-tunnel conditions 

o y~*~ & 8 for flight conditions 

• Refined grid 

o Doubled number of radial points 
o Normal surface spacing « 0.25 X baseline 
o & 3 for flight conditions 
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Computed Results 


Wind tunnel conditions 
o M 0 o = 0.6, R d = 0.8 x 10 6 , a 
o Laminar, turbulent flow 
o Comparison with experiment 

Flight conditions 
o M 0 o = 0.34, R d = 13.5 x 10 6 , 
o Turbulent flow 

Comparison with experiment 


o 
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LEX Upper Surface Flow -■ Laminar 

1^=0. 6, R_=0.8x10 6 , Alpha=20° 

0 


u> 



Grid 1 : 27x31x65 
Grid2: 31x65x65 
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LEX Upper Surface Flow — Turbulent 

M =0.6, R_=0.8x10 6 , Alpha = 20° 

00 c 



Grid 1 : 11x31 x65 
Grid2: 17x31 x65 
Grid3: 31 x65x65 


u> 

^4 

U> 


O o 

^ 33 

% a 
Q a 

el 

O 

c: Si 

jI 

3 » 



Grid 1 : 27x3 1 x65 
G rid 2 : 3 1 x65x65 
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Forebody Surface Pressure 

M. = 0.6, R f = 0.8 x 10®, a = 20° 




LEX Surface Pressure 




C13W1IJ tOW MMVl^ SSW*? 


Turbulent Surface Flow — Side View 

M rjo =0.34, R_=1 3.5x1 0 6 , Alpha= 1 9° 

Vmtrnr O 



Grid 1 : 11x31 x65 
Grid2: 17x31x65 
Grid3: 31x65x65 
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Summary of Results 


• Significant differences between laminar and turbulent solutions 
o Forebody 
o LEX upper surface 
o Body-LEX juncture on lower surface 

u> 

00 

u> 

• Turbulent solutions provide good correlation with experiment 

o Surface C p comparison with wind tunnel data 
o Surface flow comparison with flight test data 


• Convergence achieved with practical resource utilization 
o & 185,000 points 
o & 2400 cycles 
o ^2 hours of Cray-2 time 
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ATTRIBUTES OF MaGGiE COPE 



CURVE C 


Search Circle 


P 


FRINGE 

TARGET 


VERTEX 


- SUBDOMAIN GRIDS GENERATED SEPARATELY AND INDEPENDENTLY. 

- DISTANCE BETWEEN CELL CENTERS USED FOR INTERPOLATIONS. 

- REGIONS OF A FINE LEVEL GRID COMMON TO OTHERS REMOVED 

(HOLES). 

- REGIONS OF A COARSE LEVEL GRID OVERLAPPING SOLID 

SURFACES REMOVED (ILLEGAL ZONES). 

- EDGES OF HOLES OR ILLEGAL ZONES ARE INTERGRID BOUNDARIES. 

- 3-D VECTOR OPERATIONS TO LOCATE 

CELLS IN A HOLE OR AN ILLEGAL ZONE. 

CELLS OF OUTER BOUNDARY OF OVERLAP REGIONS 
(FRINGE CELL). 

IMMEDIATE-NEIGHBOR CELLS OF A HOLE OR AN 
ILLEGAL-ZONE CELL (FRINGE CELL ). 

- INFORMATION HEXAHEDRONS AROUND FRINGE CELLS FORMED. 


- CELLS IN A HOLE BUT NOT IN AN ILLEGAL ZONE UPDATED 


FROM 


FINE LEVEL 
EQUAL COARSE LEVEL 


OF OTHER GRIDS 


-COEFFICIENTS FOR TRILINEAR INTERPOLATIONS. 


-VECTORIZED DATA STRUCTURES AS PREPROCESSOR FOR VUMXZ3 


-DEVELOPED FROM CHIMERA. 

-CHARACTERISTIC INTERGRID BOUNDARY CONDITIONS. 
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ATTRIBUTES OF VUMXZ3 CODE 

- DEVELOPED FROM CFL3D CODE, HENCE INCLUDES ALL CFL3D ATTRIBUTES. 

-ALL BUT CROSS-DERIVATIVE VISCOUS TERMS. 

- BALDWIN-LOMAX TURBULENCE MODEL MODIFIED FOR VORTICAL FLOWS, MULTIPLE-WALLS, AND 

TURBULENT MEMORY RELAXATION FOR WAKES AND SHEAR LAYERS. 

- BLOCK OR DIAGONAL INVERSIONS ACCOMODATE HOLES AND ILLEGAL ZONES FOR OVERLAPPED GRIDS. 

- AVOIDS INFORMATION POLLUTION NEAR HOLES OR ILLEGAL ZONES FOR HIGHER ORDER SCHEMES. 

- NULLIFIES WEIGHT OF CONTRIBUTIONS FROM ILLEGAL ZONES DURING 3-D MULTIGRID PROLONGATION. 

- INTERGRID INFORMATION EXCHANGED BY TRILINEAR INTERPOLATIONS 

COUPLED WITH CHARACTERISTIC BOUNDARY CONDITIONS. 

- COMBINED ZONAL AND OVERLAPPED EMBEDDING. 

- DEMONSTRATIVE CASES : 

SUPERSONIC FLOW PAST A BLUNT-NOSE CYLINDER (LVD=6.7) 


AT 32 -ANGLE OF ATTACK ON 

SINGLE C-OGRID. 

C-0 GRID EMBEDDED IN CARTESIAN GRID. 
C-0 GRID PATCHED TO H-0 GRID. 

SUPERSONIC FLOW PAST AN OGIVE-NOSE CYLINDER (UD=1 8) NEAR 


FLAT PLATE. 
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EFFECTS OF SHAPE AND NCDENCE ON FLOW PAST A CYLINDRICAL SECTION 

Mi = 1.6 Rex = 2x10” / ft 
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EMBEDDED GRIDS FOR BLUNT-NOSE-CYLINDER WITH STING 


C-0 GRID : 73x65x57 
CARTESIAN GRID : 81x73x73 



overlap of c-o with cartesian 
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MACH CONTOURS OF FLOW (M=1.6 , RE=2E6/FT) PAST OGIVE-NOSE CYLINDER (L/D=6.7 , alpha=32) 


flflCH NUMBER 


CONTOUR LEUELS 






0*25000 

0.27500 

0.32500 
0.3SOOO 
0.37500 
0. -40000 
0.42500 
0.45000 
0.47500 
0.50000 
0.52500 
0.55000 
0.57500 
0.60000 
0.62500 
0.65000 
0.67500 



1.600 

32.00 DEG 
5. 00x10**5 

1.0 

40x63x63 

68x65x45 


" 18 CH 

| 

Re 
TIME 
GRID 
6R ID 
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CROSSFLOW DENSITY CONTOURS OF FLOW PAST A BLUNT-NOSE-CYLINDER 

Mx = 1.6 Rex = 2x10* / ft 




IE6/FT ) t| =0 LEESIDE, n = 180 WINDSIDE 



X/D 


OF POOR QUALITY 


PRESSURE CONTOURS ON THE SYMMETRY PLANES OF BLUNT-NOSE-CYLINDER 


Mx = 1.6 Rex = 2x10" / ft a = 32° 



(a) composite grid 
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(c) Cp on leeside 
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0.379SO 

0.39600 


0.41250 
U . -»2U(Hl 


DENSITY CONTOURS OF FLOW (M=2.86 , RE=2E6/FT) PAST OGIVE-NOSE CYLINDER 


(L/D=18) NEAR A FLAT PLATE 


JNTOUR LEUELS 


0. 13200 
0. 1 4850 
0.16S00 
0. 18150 
0. 19800 
0.21450 
0.23100 
0.24756 


2.860 
0.00 DE6 
2. 02x10*#5 

to 

83x57x53 

73x65x17 


*1RCH 
BLPHft 
Re 
TME 
GRID 1 
GRID 2 



FLOW OVER OGIVE-NOSE-CYLINDER NEAR FLAT PLATE 
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Pressure Coefficient Contours on the surface of the flat plate 



Skin Friction Pattens near the flat plate surface 



(a) two cells above the surface 



leading edge 
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ATTRIBUTES OF EMCAV3 


- 3-D REYNOLDS-AVERAGED NAVIER-STOKES EQUATIONS FOR UNSTEADY COMPRESSIBLE FLOW. 

- EXPLICIT, UNSPLIT, PREDICTOR-CORRECTOR TIME INTEGRATION (MACCORMACK). 

-FINITE DIFFERENCE SPATIAL DISCRETIZATION. 

- SECOND ORDER ACCURATE IN TIME AND SPACE. 

-VECTORIZED DATA STRUCTURE. 

- DEVELOPED FROM SCRM32. 

- ALGEBRAIC GRID GENERATOR. 


4 ^ 
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- BALDWIN-LOMAX TURBULENCE MODEL MODIFIED FOR 


VORTICAL FLOWS. 

MULTIPLE WALLS. 

TURBULENT MEMORY RELAXATION FOR SHEAR LAYER. 


- IMPLICIT MARCHING SOLUTION FOR 2-D COMPRESSIBLE BOUNDARY LAYER PROFILE. 


- FOURIER TIME SERIES ANALYSIS FOR CAVITY ACOUSTICS. 


- DEMONSTRATIVE CASES 


n-Di 

SUBSONIC 

r~o° 

“1 

DEEP 

La-oJ 

TRANSONIC 

FLOWS AT “ 

1 4 c 

,° YAW PAST 

TRANSITIONAL 

SUPERSONIC 

L w 

> J 

SHALLOW 


CAVITIES. 



FREC^tnS PAGE BLANK NOT FILMED 


INSTANTANEOUS STREAMLINES OF FLOW PAST RECTANGULAR CAVITIES 
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M x =0.58 L/D=1 1.7 Re=5.02x10‘/m 
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FREQUENCY SPECTRA OF SOUND PRESSURE LEVEL ON THE CAVITY FLOOR 
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CONCLUSIONS 


-3-D CFD CAPABILITIES FOR UNSTEADY CAVITY FLOWS. 

- TIME-AVERAGED AND TIME SERIES ANALYSES OF CAVITY FLOWS : 

EFFECTS OF LENGTH-TO DEPTH, FLOW REGIMES, YAW ANGLES. 

- 2-D ANALYSES OF FLOWS PAST VARIOUS CYLINDRICAL SECTIONS AT VARIOUS ANGLES OF ATTACK. 

- 3-D ANALYSES OF LAMINAR AND TURBULENT FLOWS PAST A BODY OF REVOLUTION 

AT VARIOUS ANGLES OF ATTACK. 

- 3-D CFD CAPABILITIES DEVELOPED FOR VISCOUS FLOWS ABOUT 

COMPLEX AND/OR MULTICOMPONENT CONFIGURATIONS. 

- COMBINED ADVANTAGES : GEOMETRICALLY CONSERVATIVE, MINIMALLY DISSIPATIVE 

: NUMERICALLY AND COMPUTATIONALLY EFFICIENT. 

: FLEXIBILITY IN BODY GEOMETRIES AND GRID TOPOLOGIES. 

- FLUX CONSERVATION ACROSS INTERGRID BOUNDARIES NEED FURTHER STUDY. 

- ANALYSES OF AERODYNAMIC INTERFERENCE OF A STORE NEAR A FLAT PLATE (EXTERNAL CARRIAGE). 

- WORK-IN-PROGRESS : CFD SOLUTIONS FOR FLOWS PAST AN OGIVE-NOSE-CYLINDER WITH L-SHAPE 

STING NEAR A CAVITY AND CODE VALIDATATION WIND TUNNEL EXPERIMENTS. 
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TranAir: Recent Advances and Applications 

Michael D. Madson 

NASA Ames Research Center, Moffett Field f CA 
March 8, 1989 


Abstract 


TranAir is a computer code which solves the full-potential equation for transonic flow about very general 
and complex configurations. Piecewise flat surface panels are used to describe the surface geometry. This 
paneled definition is then embedded in an unstructured cartesian flow field grid. Finite elements are used 
in the discretization of the flow field grid in a manner which is fully conservative and 2 n<i -order accurate. 
Since geometries may be defined with relative ease, and since the user is not involved in the generation of 
the flow field grid, computational results may be generated rather quickly for a wide range of geometries. 
For transonic cases in the cruise angle-of-attack range, TVanAir has generated results which are in generally 
good agreement with both Euler results and wind tunnel data. A typical transonic case runs in 1-2 CPU 
hours on a Cray X-MP. For subcritical cases, the code runs in 15-30 CPU minutes, even for geometries 
in which several thousand surface panels are used in the definition. This ability to rapidly and accurately 
provide both subsonic and transonic predictions about very complex aircraft configurations gives TVanAir 
the potential of being a very powerful and widely used design tool. 


Acknowledgements 
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dedicated work in the development of the code, and for providing many of the results presented in the charts 
which follow. 
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TRANAIR: RECENT ADVANCES AND 

APPLICATIONS 

MICHAEL D. MADSON 


NASA Ames Research Center 


ACKNOWLEDGEMENTS 

Forrester T. Johnson, Satish S. Samant, David R Young, 
Robin G. Melvin, and John E. Bussoletti 


OUTLINE 

• OBJECTIVE • APPLICATIONS 

• APPROACH • OBSERVATIONS 

• RECENT ADVANCES • POSSIBLE FUTURE DIRECTIONS 


MDM-1/15 



OBJECTIVE 


TO DEVELOP AND VALIDATE A COMPUTATIONAL METHOD WHICH 
ELIMINATES THE USE OF SURFACE-CONFORMING GRIDS IN THE 
ANALYSIS OF COMPLEX AIRCRAFT CONFIGURATIONS IN THE 
TRANSONIC FLOW REGIME. 



APPROACH 


GEOMETRY DEFINED BY 
SURFACE PANELS 

EMBEDDED IN UNIFORM 
CARTESIAN GRID 

LOCAL GRID REFINEMENT 
CAPABILITY BASED ON 

• LOCAL PANEL SIZE 

• USER INPUT 

SOLUTION PROCESS 

• GRID DISCRETIZED W/FINITE 
ELEMENTS, CONSERVATIVE, 

2ND ORDER 

• SET OF NONLINEAR ALGEBRAIC 
EONS SIMULATE FULL-POTENTIAL 

• 1ST ORDER DISSIPATION FOR 
SHOCK CAPTURING 

• SET OF EONS SOLVED BY 
ITERATIVE PROCESS 





415 


RECENT ADVANCES 


• LOCAL GRID REFINEMENT 

• CURRENTLY REQUIRES SUBSTANTIAL USER INPUT 

• SOLUTION ADAPTIVE REFINEMENT NEARLY COMPLETE 


IMPROVEMENTS IN SOLUTION PROCESS 

• DYNAMIC DROP TOLERANCE IN SPARSE MATRIX SOLVER 

• “SHOCK SMEARING” FOR CONVERGENCE IMPROVEMENTS 

• REGIONS OF DIFFERING T T , P T 

• MORE GENERAL BOUNDARY CONDITION OPTIONS 

• SUPERSONIC FREESTREAM) 


HIERARCHICAL MULTIGRID ) 


CURRENT EFFORTS 



APPLICATIONS 


F-16A W/TIP MISSILES, FUEL TANKS 
BOEING 747-200 

OBLIQUE WING RESEARCH AIRCRAFT 
GENERIC FIGHTER 
ONERA M6 

AXISYMMETRIC NACELLE 

4 

ADVANCED TURBOPROP MODEL 



F- 16 A W/TIP MISSILES, FUEL TANKS 
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F-16A W/TIP MISSILES, FUEL TANKS 


Moo = 0.9, a = 4° 



MDM-6/15 



F-16A, TIP MISSILE EFFECT 

Moo = 0.6, a = 4°, rj = 0.95° 
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BOEING 747-200 





GENERIC FIGHTER 


I 


i 
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GENERIC FIGHTER 

x, = 0.8, a = 5°, WING/BODY 


EXPERIMENT 

TRANAIR 


-.9 

-.8 





ONERA M6 

REFINED GRID AT ROOT SECTION 
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ONERA M6 WING 


Moo = 0.84, a = 3.06° 




PLANE OF SYMMETRY 




t) = 0.70 




OBSERVATIONS 


SHOWS PROMISE OF BEING EXCELLENT DESIGN TOOL 

• EXPECTED TO REPLACE PANAIR FOR ANALYSIS OF COMPLEX CONFIGURATIONS 

• LINEAR CASES TAKE 15-30 CPU MINUTES ON CRAY X-MP 

• NONLINEAR CASES TAKE 1-2 CPU HOURS ON X-MP 

• FOR CASES W/LIMITED VORTICAL FLOW, GOOD AGREEMENT W/FL057 

• INITIAL Y-MP RUNS INDICATE 25-40% IMPROVEMENT IN CPU TIME 

AREAS FOR IMPROVEMENT 

• FURTHER MODS COULD REDUCE CPU TIME BY A FACTOR OF 2 

• BOUNDARY LAYER COUPLING FOR MORE ACCURATE SHOCK LOCATION 

• ADDITION OF “WAKE CAPTURING” TERM FOR VORTICAL FLOW TRACKING 



POSSIBLE FUTURE DIRECTIONS 


MODIFIED FULL-POTENTIAL 

• SOURCE TERM ADDED TO FULL-POTENTIAL EQN FOR WAKE TRACKING 

• VISCOUS COUPLING 

EULER EQUATIONS 

• MAINTAIN CURRENT GEOMETRY, GRID FORMAT 

• SOLUTION ADAPTIVE REFINEMENT FOR ACCURATE VORTEX CAPTURING 

NAVIER-STOKES?? 

• PROBABLY REQUIRES “PSEUDO-GRID” FOR BOUNDARY LAYER SOLUTION 



SESSION VI 
ROTORCRAFT 


Chairman: 

W. J. McCroskey 

U.S. Army Aeroflightdynamics Directorate 
NASA Ames Research Center 
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NUMERICAL SIMULATION OF ROTORCRAFT* 


W. J. McCroskey 

U. S. Army Aeroflightdynamics Directorate-AVSCOM 

and 

NASA Ames Research Center, Moffett Field, California 


The objective of this research is to develop and validate accurate, user- 
oriented viscous CFD codes (with inviscid options) for three-dimensional, 
unsteady aerodynamic flows about arbitrary rotorcraft configurations. This 
effort draws heavily from the supercomputer capabilities of the National 
Aerodynamic Simulation project, and it will provide significantly better 
design and analysis tools to the rotorcraft industry. Better vehicles can be 
designed at lower cost, with less expensive testing, and with less risk. 

Unsteady, three-dimensional Euler and Navier-Stokes codes are being 
developed, adapted, and extended to rotor-body combinations. Flow solvers 
are being coupled with zonal grid topologies, including rotating and 
nonrotating blocks. Special grid clustering and wave-fitting techniques have 
been developed to capture low-level radiating acoustic waves. 

Significant progress has been made in computing the propagation of acoustic 
waves due to the interaction of a concentrated vortex and a helicopter airfoil. 
In this study, the need for higher-order schemes was firmly established in 
relatively inexpensive two-dimensional calculations. In three dimensions, 
the number of grid points required to capture the low-level acoustic waves 
becomes very large, so that large supercomputer memory becomes essential. 

Good agreement was obtained between the numerical results obtained with a 
thin-layer Navier-Stokes code and experimental data from a model rotor. In 
addition, several nonrotating configurations that are sometimes proposed to 
simulate rotor blade tips in conventional wind tunnels were examined, and 
the complex flow around the radical tip shape of the world's fastest 
helicopter is under investigation. These studies demonstrate the flexibility 
and power of CFD to gain physical insight, study novel ideas, and examine 
various possibilities that might be difficult or impossible to set up in 
physical experiments. 

As a prelude to studies of rotor-body aerodynamic interactions, a 
preliminary grid topology and moving-interface strategy has been developed. 
A new Euler / Navier-Stokes code using these techniques computes the 
vortical wake directly, rather than modeling it, as in most previous rotorcraft 
studies. Several hover cases were run for conventional and advanced- 
geometry blades. Numerical schemes using multi-zones and/or adaptive 
grids appear to be necessary to simulate the complex vortical flows in rotor 
wakes. 

Although major improvements both in supercomputers and in codes will be 
required, the present trends and rate of progress indicate that practical 
computations of rotor-body combinations will be feasible in the mid-1990's. 


♦This research is performed by the Rotorcraft CFD Group, consisting of 
James Baeder, Ryan Border, Earl Duque, G.R. Srinivasan, and Sharon 
Stanaway, whose contributions are gratefully acknowledged. 
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NUMERICAL SIMULATION OF ROTORCRAFT 


OBJECTIVE: 

• DEVELOP AND VALIDATE CFD CODES FOR 3-D 
VISCOUS FLOWS ABOUT ARBITRARY ELASTIC 
ROTORCRAFT CONFIGURATIONS 

APPROACH: 

• DEVELOP AND VALIDATE EULER AND NAVIER- 
STOKES CODES FOR FUTURE NAS 
SUPERCOMPUTERS 



BACKGROUND: 

• PRESENT DESIGN AND ANALYSIS TOOLS FOR ROTORCRAFT ARE INADEQUATE 

• TRIAL-AND-ERROR TESTING IS EXPENSIVE AND TIME-CONSUMING 

• FOREIGN COMPETITION IS GROWING RAPIDLY 

• CFD TECHNOLOGY FOR ROTORCRAFT LAGS FIXED-WING DEVELOPMENTS BY YEARS, 

BUT FUTURE SUPERCOMPUTERS WILL PERMIT REALISTIC ROTORCRAFT APPLICATIONS 


NUMERICAL SIMULATION OF ROTORCRAFT 

I. Very Difficult Problems 

II. We're Doing Great Work 

Specific Examples 

III. We Have Great Plans 

Complete Aeroelastic Rotor- Body Combinations , etc. 

IV. BUT... 

Hardware 
Software 
Algorithms 
Grids 

TUrbulence Model 

Summary and Conclusions 

Supercode R C 2 2 2 




V. 


ARMY/NASA ROTORCRAFT CFD PROGRAMS 


• INDIVIDUAL COMPONENTS 

1. TRANSONIC AIRFOIL CHARACTERISTICS 

2. BLADE-VORTEX INTERACTIONS 

3. ROTOR TIP-VORTEX FORMATION 

4. 3-D ACOUSTIC PROPAGATION 

i 

• COUPLED FINITE-DIFFERENCE CODES AND WAKE MODELS 


• ISOLATED ROTOR (IMA VI ER STOKES ) 


• ROTOR-BODY COMBINATIONS 







Transonic Airfoil-Vortex Interaction 


Formation of Acoustic Wave — SCI 095 



Computations Graphics : James D. Baeder, Army AFDD / NASA Ames 














BERP ROTOR 

a = 20°, Re = 1.5X10 6 , M = 0.2 



E.P.N. DUQUE, U.S. ARMY AEROFUGHTDYNAMICS DIRECTORATE. AVSCOM 
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THE BRITISH EXPERIMENTAL ROTOR PROGRAM BLADE 

M = 0.2, a = 13°, Re = 1.5 X 10 6 , NONROTATING 


— CALCULATION 
a EXPERIMENT 


i i 1 — 

1 .5 0 

x/c 


.5 

X/C 


TIP SHAPE COMPARISON 

M=0.6, alpha=6 degrees, Inviscid 





MflCH NUMBER 

1.50 

1.25 

1.00 

0.75 

0.50 

0.25 

0.00 




EULER HOVERING ROTOR CALCULATIONS 

WITH AND WITHOUT COMPUTED VORTEX WAKE 


M t = 0.794, 0 C = 8° 


WAKE 


TIP 

VORTEX 




X/C X/C X/C 


COMPUTATIONAL GRIDS FOR ROTORCRAFT 


* 



ROTATING OUTER GRID 

ROTATING INNER GRID 


ROTATING BLADE 


ROTATING INNER GRID 


COMPUTER SPEED AND MEMORY REQUIREMENTS 

FOR 1 HOUR RUNS 
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What Can We Do ? 


1. Accept Longer Run Times 

2. Speed Up the Hardware 

3. Change the Hardware and Software 

Different Architectures, 

Different Operating Systems, 
Different Languages, 

New / Improved Coding 

4. Improve the Algorithms 

Increase Stability -» increase A t 
Reduce Numerical Dissipation 

5. Use Dynamic, Solution — Adaptive Grids 

6. Simplify the Turbulence Model 



SUPERCODE R C 2 2 2 

Two Blades + Body : 

10 6 grid pts/blade + 0.5xl0 6 for body 

200 M'worcLs , 2 Gflops ( 4 psec/g rid pt/time step ) 

Major Improvements Are Required 

20 Times Faster than Today's Unsteady Navier— Stokes Codes 
-> At 10 times larger ( 1° azimuth per time step ) 

— > Flow solver 2 times faster 

Rotating and Nonrotating Grid Zones 

Solution — Adaptive Grids, Minimum, Artificial Dissipation 
Near — Wake Vortex Capturing, Far — Wake Modeling 
Improved Transition Modeling and Separation Prediction 
Flow Solver Coupled with Finite — Element Structural Model 



SPEED, MFLOPS 


ROTORCRAFT CFD PROJECTIONS 


FOR 1 HOUR RUNS 





. i . ' ‘ * • , • , „ ' v. • - : 

SUMMARY AND CONCLUSIONS 

• CFD IS VIABLE AND USEFUL FOR ROTORCRAFT 

• FUTURE DIRECTIONS 

- Detailed Study of BERP and Other Advanced Tips 
— Rotor — Body Interactions 

— Improved Wake Computations 

— Increased Collaboration with Industry ^ , 

— New Tilt — Rotor Initiatives 

• CONCERNS AND LIMITATIONS 1 ' 

— Manpower — Trained in Both CFD and Rotorcraft 

- Fhr — Field Aeroacoustics and Structural Coupling 

- Wake Capturing vs. Wake Modeling ; 

— TUrbulence Models 

— Grids for Complex Bodies in Relative Motion 
— Code Validation, Accuracy, and Reliability 

- Computer Power — CPU and Clock Time 

- Mass Storage, Post - Processing, and Graphical Display 

of 3 - D Time - Dependent Results 

• PRACTICAL ROTOR - BODY COMBINATIONS BY MID - 1990's 



N91 -10864 


CALCULATION OF THE ROTOR INDUCED 
DOWNLOAD ON AIRFOILS 


C. S. Lee 

Sterling Federal Systems, Inc. 
NASA Ames Research Center 


ABSTRACT 

Interactions between the rotors and wing of a rotary wing 
aircraft in hover have a significant detrimental effect on its payload 
performance. The reduction of payload results from the wake of 
lifting rotors impinging on the wing, which is at -90 degrees angle of 
attack in hover. This vertical drag, often referred as download, can 
be as large as 15% of the total rotor thrust in hover. 

The rotor wake is a three-dimensional, unsteady flow with 
concentrated tip vortices. With the rotor tip vortices impinging on 
the upper surface of the wing, the flow over the wing is not only 
three-dimensional and unsteady, but also separated from the leading 
and trailing edges. 

A simplified two-dimensional model was developed to 
demonstrate the stability of the methodology. The flow model 
combines a panel method to represent the rotor and the wing, and a 
vortex method to track the wing wake. A parametric study of the 
download on a 20% thick elliptical airfoil below a rotor disk of 
uniform inflow was performed. Comparisons with experimental data 
are made where the data are available. This approach is now being 
extended to three-dimensional flows. Preliminary results on a wing 
at -90 degrees angle of attack in free stream is presented. 
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SCHEMATIC OF ROTOR/WING FLOW FIELD 



AT WING ROOT 
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TWO-DIMENSIONAL ANALYSIS FORMULATION 


• DOUBLET PANELS ON ROTOR, VORTICITY PANELS ON AIRFOIL, 

AND POINT VORTICES IN WAKE 

* 

• UNSTEADY CALCULATION 

— IMPULSIVELY STARTED FLOW 

— time stepping for following solutions 

• BOUNDARY CONDITION 

--- KNOWN normal velocity distribution on actuator disk 

— ^CONSTANT STREAM FUNCTION ALONG AIRFOIL 

— ZERO TOTAL VORTICITY IN FLOW FIELD 


T 




TWO-DIMENSIONAL ANALYSIS FORMULATION 

CONTINUED 


KUTTA CONDITION 


- ROTOR WAKE STRENGTH 
DETERMINED BY TOTAL 
PRESSURE DIFFERENCE ACROSS 

SLIPSTREAM : Tr= y VrAt 

— AIRFOIL WAKE STRENGTH 
RELATED TO STRENGTH OF 
BOUND VORTICITY AT 

SEPARATION POINT : Ta = ysVsAt 

TOTAL PRESSURE VARIATION 
ACROSS ROTOR SLIPSTREAM 
AND AIRFOIL WAKE 

APr= p Vr y 




NACA 64A223 AIRFOIL (XV-15 WING) IN FREE STREAM 

-90 DEGREE ANGLE OF ATTACK 
25% CHORD FLAP DEFLECTED AT 45 DEGREE 




NACA 64A223 AIRFOIL (XV-15 WING) IN FREE STREAM 

-90 DEGREE ANGLE OF ATTACK 
25% CHORD FLAP DEFLECTED AT 45 DEGREE 



SURFACE PRESSURE DISTRIBUTION 



AIRFOIL/ROTOR INTERACTION : EFFECT OF ROTOR/AIRFOIL SPACING 

ELLIPTICAL AIRFOIL 


& 

U> 





AIRFOIL/rtOTOR INTERACTION : EFFECT OF ROTOR HEIGHT ABOVE GROUND 

ELLIPTICAL AIRFOIL 




THREE-DIMENSIONAL ANALYSIS FORMULATION 


CONSTANT SOURCE AND DOUBLET PANELS ON WING, 
DOUBLET PANELS IN WAKE 

UNSTEADY CALCULATION 

— IMPULSIVELY STARTED FLOW 

— TIME STEPPING FOR FOLLOWING SOLUTIONS 

WAKE CORE SIZE GROWS WITH AGE 
ro ~VF 

BOUNDARY CONDITIONS 

— FLOWTANGENCY 

— VELOCITY POTENTIAL JUMP ACROSS BODY PANEL =0/2 

KUTTA CONDITION 

WAKE STRENGTH, ^lw = JIu - |LlL 



ELLIPTICAL WING IN FREE STREAM 


-90 DEGREE ANGLE OF ATTACK 
ASPECT RATIO =4.0 



ELLIPTICAL WING IN FREE STREAM (CONTINUED) 



VIEW FROM BENEATH 


FUTURE WORK 


ROTOR MODEL 

- ACTUATOR DISK MODEL 

LINEAR DOUBLET PANELS IN STREAMWISE DIRECTION FOR ROTOR WAKE 

— ROTOR BLADE MODEL 

DOUBLET PANELS ON ROTOR BLADE TO INCLUDE THE EFFECT OF BLADE 
TWIST, AND SENSE OF ROTOR ROTATION 


WAKE MODEL 

- AMALGAMATE AND REDISTRIBUTE WAKE PANELS TO REDUCE 
COMPUTATIONAL TIME 

— DISCRETIZE FAR WAKE PANELS TO MODEL OSCILLATING WAKE 



'' 93-10865 


THREE-DIMENSIONAL VISCOUS DRAG PREDICTION 
FOR ROTOR BLADES 


CHING S. CHEN 

NATIONAL RESEARCH COUNCIL 
NASA AMES RESEARCH CENTER 


SUMMARY 

The state-of-the-art in rotor blade drag prediction involves the use 
of two-dimensional airfoil tables to calculate the drag force on the 
blade. One of the most serious problems with the current methods is 
that they cannot be used for airfoils that have yet to be tested. Most 
of the drag prediction methods also do not take the Reynolds number 
or the rotational effects of the blade into account, raising doubts 
about the accuracy of the results. This project addresses these 
problems with the development of an analytical method which 
includes the shape of airfoil, the effects of Reynolds number, and the 
rotational motion of the blade. 
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OBJECTIVES 


PROVIDE AN EFFICIENT ANALYSIS FOR ROTOR 
BLADE DRAG PREDICTION 

ENABLE ACCURATE DESIGN OF NEW, ADVANCED 
ROTOR SYSTEMS WITH IMPROVED PERFORMANCE 



MOTIVATIONS 


ROTOR BLADE DRAG PREDICTION CURRENTLY USES TWO- 
DIMENSIONAL AIRFOIL TABLES 

MAJOR SHORTCOMINGS: 

• CANNOT BE USED FOR AIRFOILS NOT YET TESTED 

• NO REYNOLDS NUMBER OR ROTATIONAL EFFECTS 
INCLUDED 

CURRENT APPROACH INCLUDES: 

• AIRFOIL SHAPE 

• REYNOLDS NUMBER EFFECTS 

• ROTATIONAL MOTION OF BLADE 
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RELATION OF COMPUTATIONAL DOMAIN OF INVISCID ANALYSIS 

AND ROTOR WAKE 



THE COUPLING SCHEME OF VISCOUS, INVISCID, AND FREE-WAKE METHODS 

DIRECT APPROACH 





UNSTEADY, COMPRESSIBLE, THREE-DIMENSIONAL BOUNDARY-LAYER EQUATIONS 

ROTATING REFERENCE FRAME 

X-MOMENTUM EQUATION 

p(u t + uu x + vu y + wu z - 2w Cl - xQ 2 ) = - p x + (|iu y ) y 

Z-MOMENTUM EQUATION 

p(w t + uw x + vw y + ww z + 2u£2 - zQ. ) = - p z + (|xw y ) y 

PERFECT GAS RELATION 

p = pRT 

ENERGY EQUATION 

CpT + \(m 2 + v 2 + w 2 ) = CONST 

CONTINUITY EQUATION 

p t + (pu) x + (py)y + (pw) z =0 



THE COUPLING SCHEME OF VISCOUS, INVISCID, AND FREE-WAKE METHODS 

INVERSE APPROACH 






COMPARISON OF PREDICTED AND MEASURED DRAG COEFFICIENT 
FOR NACA 0012 TWO-DIMENSIONAL AIRFOIL FLOWS 
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COMPARISON OF Cp DISTRIBUTION WITH AND WITHOUT 
BOUNDARY-LAYER INTERACTION 
TIP MACH NUMBER = 0.9, r/R = 0.8453 





COMPARISON OF Cp DISTRIBUTION WITH AND WITHOUT 
BOUNDARY-LAYER INTERACTION 
TIP MACH NUMBER = 0.9, r/R = 0.9759 
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VELOCITY vector plot showing separation region 

AT TRAILING EDGE 

TIP MACH NUMBER = 0.925, r/R = 0.9759, AR = 13.71 
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TORQUE, FT-LBF 


COMPARISON OF PREDICTED AND MEASURED TORQUE 
FOR A NONUFITNG, HOVERING ROTOR 
TWO BLADES, RECTANGULAR NACA 0012 BLADE 
RADIUS = 3.428 FT, AR = 13.71 



FUTURE WORK 


SIMULATE NONLIFTING ROTOR FLOWS WITH 
ADVANCED BLADE PLANFORMS, AIRFOILS, AND 
TWIST DISTRIBUTIONS 

INCORPORATE FREE-WAKE ANALYSIS INTO VISCOUS- 
IN VISCID APPROACH TO SIMULATE LIFTING ROTOR 
FLOWS 

INVESTIGATE EFFECTS OF ROTATIONAL MOTION ON 
DEVELOPMENT OF BOUNDARY-LAYER 


Progress Toward the Development of an 
Airfoil Icing Analysis Capability 


§ 


Mark G. Potapczuk 
Colin S. Bidwell 
NASA Lewis Research Center 
Cleveland, Ohio 

Brian M. Berkowitz 
Sverdrup Technology, Inc. 
Mayfield Heights, Ohio 
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Progess Toward the Development of an Airfoil Icing Analysis Capability 


Af.G . Potapczvk 
C.S. Bidwell 

NASA Lewis Research Center , Cleveland , OAio 
R.J/. Berkowitz 

Sverdrup Technology , Jnc., Middleburg ffts ., OAio 

The NASA-Lewis aircraft icing analysis program is composed of three major sub-programs. 
These sub- programs are ice accretion simulation, performance degradation evaluation, and ice 
protection system evaluation. These topics cover all areas of concern related to the simulation of 
aircraft icing and its consequences. The motivation for these activities is twofold, reduction of time 
and effort required in experimental programs and the ability to provide reliable information for aircraft 
certification in icing, over the complete range of environmental conditions. In addition to the 
analytical activities associated with development of these codes, several experimental programs are 
underway to provide verification information for existing codes. These experimental programs are also 
used to investigate the physical processes associated with ice accretion and removal for improvement of 
present analytical models. The NASA-Lewis icing analysis program is thus striving to provide a full 
range of analytical tools necessary for evaluation of the consequences of icing and of ice protection 
systems. 

Recently, two of these tools were used to produce a computational evaluation of the ice accretion 
process and resulting performance changes for a NACA0012 airfoil. The ice accretion code, LEWICE, 
provided the ice shape geometry at several points in time during the simulated icing encounter. The 
predicted shapes are a function of several environmental input parameters, including airspeed, 
temperature, water droplet size and distribution, liquid water content, and duration of the encounter. 
These ice shape geometries are then used as input for a Navier-Stokes analysis code, ARC2D, which 
calculates the flowfield and determines changes in performance characteristics of the airfoil. Presently, 
there is no direct link between the two codes and all interfacing is done by the user. One of the 
objectives of the icing analysis program is to combine codes such as these into a comprehensive icing 
analysis method. Work in this area is currently underway via a number of grant supported activities. 
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STEPS IN AIRFOIL ICING ANALYSIS METHODOLOGY 



ARC2D 




AERODYNAMIC PERFORMANCE DEGRADATION DUE TO ICING 
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MULTIZONE MODEL OF ICE ACCRETION PROCESS 


ZONE B 

LARGE SURFACE ROUGHNESS 



VARYING ROUGHNESS 
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CONTROL VOLUME ANALYSIS OF ICE ACCRETION PROCESS 


CONVECTION WATER F LOW 
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COLLECTION EFFICIENCY COMPARISONS 
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COMPARISON OF CALCULATED AND MEASURED HEAT 
TRANSFER COEFFICIENTS ON SMOOTH SURFACE 
“ICED” CYLINDER MODEL 



CD-88-38182 



COMPARISON OF ICE SHAPE PREDICTIONS WITH 

AIRFOIL ICING DATA 
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EFFECT OF TURBULENCE MODEL ON NAVIER-STOKES PREDICTIONS 
OF LEADING EDGE PRESSURE DISTRIBUTIONS, a = 0° 




BALDW IN-LOMAX JOHNSON-KING 

MODEL MODEL 



0 .125 .250 


x/c 

k-e model 


CD-88-38 J 8? 



EFFECT OF BOUNDARY LAYER TRANSITION SPECIFICATION 
ON NAVIER-STOKES PREDICTED VELOCITY PROFILES IN 
SEPARATION-REATTACHMENT ZONE, a = 0° 
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COMPARISON OF ICED AIRFOIL CODE PREDICTIONS WITH 
EXPERIMENTAL MEASUREMENTS 
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CONCLUDING REMARKS 


First generation airfoil icing capability exists 

Code validation activities are ongoing 
Droplet trajectories / impingement 
Ice accretion 

Aerodynamic performance 

Supporting analytical/experimental efforts underway 
to improve physical modeling in codes 
Movies/photographs of ice accretion 
Ice surface roughness 

Extension to 3D icing analysis has been initiated 
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The Breakup of Trailing- Line Vortices 


David Jacqmin 

NASA Lewis Research Center 
Cleveland, Ohio 44135 


It is now known that Batchelor’s trailing-line vortex is extremely unstable to 
small amplitude disturbances for swirl numbers in the neighborhood of .83. We 
present results of numerical calculations that show the response of the vortex in 
this range of swirl numbers to finite amplitude, temporal, helical disturbances. 
Phenomena observed include: 1) ejection of axial vorticity and momentum from 
the core resulting in the creation of secondary, separate vortices; 2) a great in- 
tensification of core axial vorticity and a weakening of core momentum; 3) the 
production of azimuthal vorticity in the form of a tightly wrapped spiral wave. The 
second phenomenon eventually stabilizes the vortex, which then smooths and grad- 
ually returns to an axisymmetric state. The calculations are mixed spectral-finite- 
difference, fourth-order accurate, and have been carried out at Reynolds numbers 
of 1000-2000. Some linearized results will also be discussed in an attempt to explain 
the process of vortex intensification. 


PRECEDING PAGE BLANK NOT FILMED 


489 



Applications 


Axisymmetric evolution of vortices, jets, and wakes. 

3D instability of axisymmetric flows (weakly nonlinear). 

Spiral and axisymmetric vortex breakdown - strongly nonlinear instabilities. 
Acoustic excitation of swirling jets. 

Mixing caused by vortices - effects on combustion. 



Numerical Methods 


3D direct numerical simulation using primitive variables. 
Unbounded cylindrical coordinates - use of mapping r = tan A. 
Mixed spectral and finite difference methods 
Azimuthal: spectral (expm0). 

Axial: spectral (exptamz) or 2nd order FD. 

Radial: spectral (sin n 2A cos 2; A) or 4th order FD. 

Fourth order accurate (Runge-Kutta) in time. 



Batchelor’s Trailing-Line Vortices 


Defined by axial velocity = exp(— r 2 ), axial vorticity = q exp(—r 2 ). 

Model for aircraft trailing line vortices, created in laboratory. 

Very unstable in some ranges of q. q cz .85 most unstable. 

From inviscid theory, all azimuthal wavenumbers unstable. 

As n — ♦ oo Ci — *• .4. 

As n — > oo most unstable axial wavenumber — >• .52 n 
From viscous theory, c,- modified by 0(i2~ 1 ). # unstable waves 0(.R 3 / 5 ). 


Numerical Solution 


Quasi-3D calculation can capture all the most unstable modes. 
Navier Stokes equations have helical solutions: 


n=-f-oo 

fn[r t t)expin[Px - 9) 

n=-oo 


Calculations have q = .82, (3 = .52, R s* 1600. 
Calculations spectral in x and 6, fourth order in A. 



Features of Solution 


• Initially, elliptical deformation of vortex. Helical displacement from center line. 

• Kinematic deformation of vortex then produces spiral structure. 

• Axial divergence intensifies axial vorticity, is associated with azimuthal vorticity. 

• “Turbulent” patches form at vortex outer edges. 

• In later stages, viscous diffusion weakens spiral structure. Turbulent patches separate from 

vortex. Vortex core remains intense. 

• Intensified vortex returns to axisyjxmietric state. 

Axial vorticity intensifies. Axial velocity weakens. Vortex stabilizes. 


APPENDIX 


LIST OF ATTENDEES 
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